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Abstract 

I Searches for possible new quantum phases and classifications of quantum phases have been central problems in 

y—{ physics. Yet, they are indeed challenging problems due to the computational difficulties in analyzing quantum many- 

^^ body systems and the lack of a general framework for classifications. While frustration-free Hamiltonians, which 

Cn appear as fixed point Hamiltonians of renormalization group transformations, may serve as representatives of quan- 

?-H tum phases, it is still difficult to analyze and classify quantum phases of arbitrary frustration-free Hamiltonians ex- 

haustively. Here, we address these problems by sharpening our considerations to a certain subclass of frustration-free 
Hamiltonians, called stabilizer Hamiltonians, which have been actively studied in quantum information science. We 
^f) propose a model of frustration-free Hamiltonians which covers a large class of physically realistic stabilizer Hamilto- 

nians, constrained to only three physical conditions; the locality of interaction terms, translation symmetries and scale 
symmetries, meaning that the number of ground states does not grow with the system size. We show that quantum 
Q, phases arising in two-dimensional models can be classified exactly through certain quantum coding theoretical oper- 

' ators, called logical operators, by proving that two models with topologically distinct shapes of logical operators are 

^ always separated by quantum phase transitions. 
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1. Introduction 

In condensed matter physics, studies on quantum phases in two-dimensional spin systems have been central in 
addressing the underlying mechanisms behind correlated spin systems 1 1 1. There have also been considerable interests 
in quantum phases in two-dimensional systems from quantum information science community [2|, as topological 
phases of spin systems are useful in realizing quantum information theoretical ideas |3-6|. In this paper, we ask two 
fundamental questions concerning quantum phases arising in two-dimensional spin systems. 

(a) What kinds of quantum phases are allowed to exist in two-dimensional spin systems? 

(b) How do we classify these quantum phases? What kinds of properties distinguish different quantum phases? 

Finding exactly solvable models is a key to searching for possible quantum phases. The AKLT state, which is 
the ground state of an exactly solvable spin 1 Hamiltonian, provides useful clues about the ground state properties 
of quantum anti-ferromagnets ||7]|8l. The Toric code, the simplest exactly solvable model with topological order, 
gives insights on how topological phases emerge in correlated spin systems |[3]|5|. These Hamiltonians mentioned 
here are called frustration-free Hamiltonians since ground states can be obtained by minimizing each term locally. 
While these exactly solvable frustration-free models may not be physical, involving more than two spins at one time 
or higher order interactions, it is relatively easy to find their ground states compared with frustrated Hamiltonians 
which we encounter in physically realistic systems. A key idea is to realize that these frustration-free Hamiltonians 



may appear as low energy effective theories for real spin systems ||9l, and thus, may serve as a "bureau of quantum 
phases" with which one may approximately study the ground state properties of actual Hamiltonians. 

So far, a lot of frustration-free spin Hamiltonians with various physical properties have been discovered in a search 
for possible quantum phases arising in two-dimensional spin systems ||5][T0l. However, whether these models could 
exhaust all the possible quantum phases or not is far from obvious. A brute force approach might be to analyze all the 
possible frustration-free Hamiltonians and study their ground state properties. However, being frustration-free does 
not mean that Hamiltonians are exactly solvable since studying their ground state properties may be challenging since 
computations of order parameters to distinguish quantum phases may involve a large number of spins. Moreover, 
finding frustration-free Hamiltonians is much more difficult than finding a ground state of frustration-free Hamiltoni- 
ans since it is computationally difficult to check whether a given Hamiltonian is frustration-free or not. For example, 
even when a Hamiltonian consists only of projectors as in the AKLT Hamiltonian and the Toric code, to determine 
whether the Hamiltonian is frustration-free or not is known to be the hardest problem in the computational complexity 
class QMA (Quantum Marlin- Arthur), a quantum analog of NP ifTTlfTSl . (So, it could be a difficult problem even for 
a quantum computer to solve !). 

As for a classification of quantum phases, two frustration-free Hamiltonians may be classified through the presence 
or the absence of a quantum phase transition (QPT). In quantum many-body systems, two quantum phases with 
different physical properties are separated by a QPT, which is a sudden non-analytic change of the ground state 
properties as a result of parameter changes in the Hamiltonian bridging two quantum phases [jj. Two frustration-free 
Hamiltonians may be considered to belong to different quantum phases when a parameterized Hamiltonian connecting 
them undergoes a QPT as a parameter varies in the Hamiltonian fTT-TSl. 

While studying the ground state properties of parameterized Hamiltonians is generally difficult, QPTs occurring 
in parameterized Hamiltonians have been studied through Renormalization Group (RG) transformations which are 
certain scale transformations acting on parameterized Hamiltonians |T6l or parameterized ansatz of ground states ifTSl 
[T7I420I . The central idea in classifying quantum phases through RG transformations is to capture only the scale 
invariant physical properties of corresponding quantum phases by using fixed point Hamiltonians which are invariant 
under RG transformations. Since fixed point Hamiltonians are obtained by washing out short-range correlations, it is 
known that fixed point Hamiltonians in gapped quantum phases can be represented as frustration-free Hamiltonians 
where each term can be minimized locally without considering neighboring terms. 

However, the challenge in classifying quantum phases in correlated spin systems underlies behind the fact that 
the existence of a QPT depends on how Hamiltonians are varied from the one to the other. In particular, there 
are cases where two frustration-free Hamiltonians are connected through some parameterized Hamiltonian without 
a QPT, but are separated by a QPT for other parameterized Hamiltonian ||2TI4251 . To show that two frustration- 
free Hamiltonians belong to different quantum phases, one needs to analyze all the possible paths of parameterized 
Hamiltonians connecting them and see if they are always separated by QPTs or not. Thus, the need is to find some 
order parameter to classify quantum phases in a way independent of choices of parameterized Hamiltonians. 

Searches for possible quantum phases and their classifications lie at the heart of studies on quantum many-body 
systems. Currently, both problems however face challenges as illustrated in the above. In this paper, we make an 
attempt to overcome these difficulties concerning studies on quantum phases by combining theoretical tools developed 
in quantum information science and the notion of topology. 

Possible quantum phases and stabilizer codes: While analyses on arbitrary frustration-free Hamiltonians are 
indeed challenging, there exists a certain subclass of frustration-free Hamiltonians, called stabilizer Hamiltonians ||2l 
|26| . which play a crucial role in quantum information science. Stabilizer Hamiltonians serve as a natural architecture 
to realize quantum error-correcting codes in correlated many-body spin systems. The basic idea of stabilizer codes 
is to encode a qubit in strongly entangled and correlated states so that no small error can break the encoded qubit. 
A remarkable feature of stabilizer codes is the existence of system Hamiltonians which supports encoded qubits in 
degenerate ground states with a finite energy gap. Such system Hamiltonians, called stabilizer Hamiltonians, are 
frustration-free Hamiltonians which consist only of Pauli operators commuting with each other Since protecting a 
qubit from decoherence is essential in realizing quantum information theoretical ideas such as fault-tolerant quantum 
computation 121, studies on stabilizer Hamiltonians have been central in quantum information science. 

Stabilizer Hamiltonians are certainly limited compared with arbitrary frustration-free Hamiltonians, yet they may 
characterize a large number of quantum phases. In fact, it has been realized that many interesting spin systems may be 
described in the language of stabiUzer codes ll9l l27H32l . as the Toric code is one of the earliest examples of stabilizer 



codes discovered in a search for useful quantum codes |[3]. Also, a good news in analyzing quantum phases arising in 
stabilizer Hamiltonians is that there are many useful theoretical tools to analyze the ground state properties in stabilizer 
codes, which are originally developed in quantum coding theoretical contexts El l27l 128.1 . but are potentially useful for 
studying non-local correlations arising in many-body systems ||3Tl[32l . Thus, the analysis on quantum phases arising 
in stabilizer Hamiltonians and their classification will be the necessary first step in order to address quantum phases 
arising in arbitrary frustration-free Hamiltonians and search for all the possible two-dimensional quantum phases. 

There are several useful theoretical tools in analyzing the ground state properties of stabilizer Hamiltonians. Basic 
analysis tools for the ground state properties of stabilizer Hamiltonians were developed for the cases where a ground 
state is unique Il27l|28l . These tools were further generalized to stabilizer Hamiltonians with the ground state degen- 
eracy with an aim to study global entanglement arising in many-body spin systems such as topological order 021 . 
However, there is still an important gap between physically realistic correlated spin systems and stabilizer codes 
discussed in quantum information science community, as most stabilizer codes encode qubits in highly non-local 
Hamiltonians which are not physically realistic. 

Several authors have initiated studies on physically realistic stabilizer Hamiltonians by analyzing coding properties 
of stabilizer codes supported by local terms 13111331 [34i . However, physically realizable Hamiltonians must have 
some physical symmetries such as translation symmetries too if they are to be realized in some lattice of spins. 
Also, there is another important physical constraint which must be considered in analyzing quantum phases arising in 
stabilizer Hamiltonians. According to the underlying philosophy behind RG transformations, quantum phases must 
be characterized by non-local correlations which survive even at the thermodynamic limit (a limit where the system 
size becomes infinite) and are invariant under scale transformations 1 16|. Since quantum phases must be classified in 
a scale invariant way, stabilizer Hamiltonians need to have scale invariance too if they are to be used as candidates for 
possible quantum phases. 

In this paper, to search for possible quantum phases in two-dimensional spin systems, we begin by proposing a 
model of stabilizer Hamiltonians which have both physical realizability and scale invariance. In particular, the model 
is built on the following physical constraints. 

• Locality of interactions: A system of qubits, defined on a two-dimensional lattice, is governed by a stabilizer 
Hamiltonian which consists only of local interactions. 

• Translation symmetries: The Hamiltonian is invariant under some finite translations of qubits. 

• Scale symmetries: The number of degenerate ground states of the Hamiltonian does not depend on the system 
size. 

We call stabilizer Hamiltonians which satisfy the above three physical constraints Stabilizer codes with Translation 
and Scale symmetries (STS models). 

Classification of quantum phases and topology: The goal of this paper is to analyze and classify quantum phases 
arising in STS models completely. In the analysis and classification of quantum phases arising in STS models, quan- 
tum phases must be distinguished by some quantities or objects which are scale invariant. Here, it is useful to realize 
the similarity between the notion of topology and the classification of quantum phases by observing that geometric 
shapes of objects are scale invariant. In topology, geometric shapes of objects are classified in terms of smoothness 
and non-analyticity. Two objects are considered to be the same when they can be transformed each other continuously, 
while they are considered to be different when they can be transformed each other through only non-analytic changes 
of geometric shapes. In a similar fashion, quantum phases are classified in terms of continuous changes of ground 
states induced by changes of parameterized Hamiltonians. Two frustration-free Hamiltonians are considered to belong 
to the same quantum phase when their ground states can be connected continuously. Two frustration-free Hamilto- 
nians are considered to belong to different quantum phases when they can be connected through only parameterized 
Hamiltonians which undergo QPTs with non-analytic changes of ground states. Thus, the notion of topology and a 
classification of quantum phases are fundamentally akin to each other 

The notion of topology has been playing crucial roles in classifying physical properties of quantum many-body 
systems in various fields of physics |6, 35 1. Quantum field theory, equipped with the notion of topology, has been 
proposed in high energy physics where physical properties of systems depend on topological characteristics of the 
geometric manifolds where the systems are defined |36|. As a physical realization of such topological theories in 



low-dimensional many-body systems, quantum hall systems defined on a geometric manifold were studied and the 
notion of topological order has been introduced in order to discuss and characterize topology-dependent non-local 
correlations 1.37.1 . The notion of topology appears not only in real physical space, but also in the momentum space 
too where QPTs are triggered by topological changes in the spectrum |[38l . Finally, when Hamiltonians have several 
parameters, the parameter space may have non-trivial topological structures which may result in interesting behaviors 
of geometric (Berry) phases ||39ll40l . 

With this intimate connection between topology and quantum phases, we wish to classify quantum phases arising 
in STS models by introducing the notion of topology into quantum coding theoretical tools developed in studies of 
stabilizer Hamiltonians. What plays a central role in analyses on the ground state properties in stabilizer Hamiltonians 
are certain operators, called logical operators ||2l, which are PauU operators commuting with the system Hamiltonian, 
but act non-trivially inside the ground state space. A useful empirical knowledge commonly shared in quantum 
information science community is that ground states are highly entangled when logical operators are supported by a 
large number of spins for a fixed system size. In particular, it has been pointed out that geometric non-localities of 
logical operators may characterize global entanglement arising in ground states of stabilizer Hamiltonians, such as 
topological order 1.31 . .32 1. 

In this paper, we classify quantum phases arising in STS models by introducing the notion of topology into 
geometric shapes of logical operators. In particular, we achieve the following two programs. 

• Exact solution of the model: We solve the STS model exactly by identifying all the possible geometric shapes 
of logical operators. We establish the connection between the scale invariant ground state properties, such as 
topological order, and geometric shapes of logical operators completely^ 

• QPTs and topology in logical operators: We show that STS models with different geometric shapes of logical 
operators are always separated by a QPT We find that the existence of a QPT originates from changes of 
geometric shapes of topologically distinct logical operators. We show that topological structures of geometric 
shapes of logical operators can be used as order parameters to distinguish quantum phases arising in one and 
two-dimensional STS models. 

Organization and contents: The paper is organized as follows. In Section l2] we begin by introducing basic 
notions used throughout this paper by giving introductions to QPTs and stabilizer codes. A definition of quantum 
phases and a criteria for their classification are given. Also, logical operators, which are central in the analyses on 
stabilizer Hamiltonians, are reviewed. In Sectionl3] we introduce the STS model by imposing physical constraints and 
scale invariance. Then, we develop a basic analysis tool, concerning a certain symmetry on logical operators, which 
arises as a result of physical symmetries and scale invariance we impose on stabilizer Hamiltonians. In Section[4]and 
SectionIS] we analyze quantum phases arising in one and two-dimensional STS model and show that different quantum 
phases are characterized by geometric shapes of logical operators. We find that the existence of a QPT originates from 
changes of geometric shapes of topologically distinct logical operators. We show that topological characteristics of 
geometric shapes of logical operators can be used as order parameters to distinguish quantum phases arising in one 
and two-dimensional STS models. 

Here we make some comments on the contents of the paper First, the paper is written in a way accessible to 
both condensed matter physicists and quantum information scientists, providing introductions to quantum coding 
theory and the notion of QPTs and topological order. In particular, the paper is presented in a self-consistent way. 
Second, while a mathematical rigor is a virtue, but so is the physical intuition. Thus, we have presented all the 
derivations and calculations of logical operators in appendices (from [Appendix A| to [Appendix D| i so that we can 
concentrate on physical discussions concerning quantum phases arising in STS models. Finally, while the locality 
of interaction terms and translation symmetries are commonly used as physical constrains in studies of many-body 



spin systems, the importance of scale symmetries may not be obvious. In Appendix E we expand our analyses to 
stabilizer Hamiltonians without scale symmetries in order to study the role of scale symmetries. We point out the 
connection between translation symmetries of ground states and scale symmetries by analyzing the phenomena called 
weak breaking of translation symmetries |l9l- 



' Here, we claim that the model is exactly solved by finding logical operators since logical operators allow us to determine essential physical 
properties of the model as we shall see in the main discussion of this paper. 



2. Review: classification of quantum pliases and stabilizer Hamiltonians 

We begin with quick reviews of theoretical notions which are central to discussions in this paper. In Section [ZT] 
we review how different quantum phases are classified through frustration-free Hamiltonians with some examples of 
quantum phase transitions (QPTs). Then, we state the criteria for classifying quantum phases clearly for the sake of 
later discussions. The notion of scale invariance and its role in a classification of quantum phases are reviewed and 



discussed. In Section 2.2 we give a review of stabilizer codes which are quantum error-correcting codes realized by 
a certain subclass of frustration-free Hamiltonians. Some theoretical tools which become central in analyzing STS 
models are introduced here. 

Here, we make a comment on the notations used throughout this paper. A state of a qubit is described as a 
superposition of |0) and |1), and Pauli operators are represented by alphabets X, Y, and Z in this paper. 

2.1. Classification of quantum phases through frustration-free Hamiltonians 

In a quantum many-body system at zero temperature, physical properties of ground states may change abruptly 
as a result of parameter changes in the Hamiltonian [IJ. In such a quantum phase transition (QPT), quantum phases 
may be characterized by frustration-free Hamiltonians. Here, we first review two parameterized Hamiltonians which 
undergo QPTs and discuss how quantum phases can be classified. Then, we state a criteria to classify quantum phases 
in this paper clearly. 

The models of QPTs we review in this subsection are known to be exactly solvable, meaning that one can char- 
acterize the ground state properties and obtain the energy spectrum through analytical solutions of the Hamiltonians. 
The solutions to models may be found in the original literature |i2 1 J and some standard textbooks on condensed mat- 
ter physics (for example, see |[T]). Solutions and discussions on these models with quantum information theoretical 
viewpoints may be found in Il22l |4T|| . 

An example of QPTs and frustration-free Hamiltonians: We consider the Ising model in a transverse field in 
one dimension with periodic boundary conditions: 

H(e) = -( 1 - e) ^ Z<-')Z(^+') - ej^ X*^' (1) 

where Z*-'^ acts on j-th qubit (spin 1/2 particle). The system consists of A^ qubits (spin 1/2 particles). At e = 0, the 
system Hamiltonian is 

HiO) ^-Y^Z^J^Z^J^^^ (2) 



and the ground states are |0 ■ ■ ■ 0) and |1 ■ ■ ■ 1). At e = 1, the system Hamiltonian is 

H(l) =-Xi^''' (3) 

j 

and the ground state is | + ■ ■ ■ H-) where |-i-) = ~k(IO) + |1))- It is known that physical properties of ground states 
drastically change around e - 1/2. When < e < 1/2, physical properties of a ground state are close to ones of H(Q). 
When 1/2 < e < 1, physical properties of a ground state are close to ones of //(I). This drastic change of physical 
properties at e = 1/2 is called a QPT. Two Hamiltonians H(Q) and //(I) are separated by a QPT at e = 1/2, and may 
represent different quantum phases (Fig.fTJa)). The existence of a QPT can be detected by considering the energy gap 
between a ground state and excited states. Let us consider how the energy gap changes around the transition point 
in this model. Then, at the transition point e = 1/2, the energy gap A(e) between a ground state and excited states 
becomes zero at the thermodynamic limit: A(l/2) -^ for A^ — > oo. Thus, the vanishing energy gap serves as an 
indicator of a QPT. 

Since a ground state of the original parameterized Hamiltonian H(e) cannot be obtained by minimizing each term 
in the Hamiltonians except at e = 0, 1, //(e) is said to he frustrated. In general, finding a ground state of a frustrated 
Hamiltonian is computationally difficult. However, Hamiltonians //(O) and //(I) are not frustrated, meaning that their 
ground states can be obtained by minimizing each term in the Hamiltonian independently. In this example, physical 
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Figure 1: (a) Two quantum phases separated by a QPT. Black circles represent frustration-free Hamiltonians and a blank circle represents a QPT. 
Different quantum phases are represented by frustration-free Hamiltonians. (b) A path dependence of the presence of a QPT. Two frustration-free 
Hamiltonians are considered to belong to the same quantum phase when there exists a parametrized Hamiltonian which connects them without 
closing the energy gap. 

properties of frustration-free Hamiltonians i/(0) and //(I) are easy to analyze. Then, physical properties of the original 
Hamiltonian //(e) may be studied by adding perturbations to frustration-free Hamiltonians //(O) and //(I). Thus, //(O) 
and //(I) may serve as easily solvable reference models to characterize quantum phases (Fig.[TJa)). 

Path dependence of the presence of a QPT: Then, we are tempted to call two quantum phases different if they 
are separated by a QPT. However, there is a subtlety in this classification of quantum phases since the presence of a 
QPT may depend on how we change the Hamiltonian. Let us see an example: 

//(e) = -(1 - e) J] Z*^'Z<-'+i' - e ^ X<»X<>+'\ (4) 

! J 

This model exhibit a QPT at e = 1/2 where the energy gap A(e) becomes zero. Then, one might hope to consider 
two quantum phases at < e < 1/2 and at 1/2 < e < 1 different. However, there exists another parameterized 
Hamiltonian //'(e) which connects //(O) and //(I) without closing an energy gap between ground states and excited 
states (Fig. [TJb)): 

//'(e) = -^M(e)<^'M(e)*-'+'\ M(e)<^' = (1 - e)Z<^n eX*-'' (5) 

;■ 

with //'(O) - //(O) and //'(I) - //(I). It is easy to see that the energy gap of the above Hamiltonian //'(e) is 
independent of e, and remain finite. Also, since //(O) can be obtained by rotating the axis of each qubit in //(I), we 
may consider two Hamiltonians to be in the same quantum phase. 

Classification of quantum phases through frustration-free Hamiltonians: As we have seen in the above ex- 
amples, the existence of a QPT depends on paths of parameterized Hamiltonians which connect two frustration-free 
Hamiltonians. A standard way to classify quantum phases is to consider two Hamiltonians belonging to the same 
quantum phase when they can be transformed each other without closing the energy gap through some path of a pa- 
rameterized Hamiltonian |fl3HT5l . Here, following this spirit, we shall state the criteria to classify quantum phases 
arising in frustration-free Hamiltonians as follows (Fig.l2]). 

• Two frustration-free Hamiltonians //^ and Hg belong to different quantum phases if and only if there is no 
parameterized Hamiltonian //(e) which connects Ha and Hb without closing an energy gap or changing the 
number of ground states. 

Here, by "parameterized Hamiltonians", we mean the following: 

• Parameterized Hamiltonians //(e) consist only of local terms which change continuously with some external 
parameter e, and amplitudes (norms) of local terms do not depend on the system size. 

By "without closing an energy gap", we mean that 

• The energy gap between degenerate ground states and a first excited state is finite even at the thermodynamic 
limit. 



Thus, two different quantum phases are always separated by a QPT regardless of choices of parameterized Hamilto- 
nians H(e). 

Here, we make some comments on the validity of the above classification of quantum phases. While the vanishing 
energy gap may be a signature of QPTs, the original definition of a QPT is a non-analytic change of the ground state 
properties. In fact, there are examples of parameterized Hamiltonians whose ground state energy changes continu- 
ously while the energy gap vanishes [23]. Thus, the connection between the energy gap and the non- analytic ity in 
ground states is not completely established. However, we use the vanishing energy gap as a criteria for classifying 
quantum phases for simplicity of discussion in this paper. In this sense, the classification presented above relies on the 
assumption that QPTs are triggered by the vanishing energy gap or the change of the number of degenerate ground 
states ID. 



Phase 1 



Phase 2 



.*. 



Phase 3 \ ^ QPT 



Figure 2: A classification of quantum phases. Black dots represent frustration free Hamiltonians in corresponding quantum phases. Frustration-free 
Hamiltonians belonging to dift'erent quantum phases are always separated by a QPT. 

RG transformations and scale invariance: QPTs in parameterized Hamiltonians have been analyzed commonly 
through RG transformations lfT6l . which are certain scale transformations. Here, we review the basic idea of RG 
transformations and discuss the relation between RG transformations and frustration-free Hamiltonians in terms of 
scale invariance. Note that while we shall describe a general idea of RG transformations, we do not give any specific 
procedure of RG transformations since detailed reviews on various procedures of RG transformations are beyond the 
scope of this paper 

RG transformation is a way to capture only the non-local correlations from many-body systems by analyzing 
how the system properties changes under scale transformations. RG transformations usually consist of two elements, 
called coarse-graining and rescaling. Coarse-graining is a process to group several spins into a single particle with 
a larger Hilbert space. In this coarse-graining, a microscopic structure of original spins inside a larger spin (coarse- 
grained spin) is completely lost and one ends up with a system consisting only of larger spins. Rescaling is a process to 
replace a coarse-grained larger spin with a original small spin by eliminating some spin degrees of freedom inside each 
coarse-grained spin, which is interacting weakly with other coarse-grained spins. By repeating scale transformations 
consisting of coarse-graining and rescaling, only the non-local correlations extending all over the lattice may survive. 

When quantum phases are classified through RG transformations, the key idea is the use of Hamiltonians which 
are invariant under scale transformations, called fixed point Hamiltonians. By applying RG transformations to some 
parameterized Hamiltonian H(e), one can systematically obtain Hamiltonians which are invariant under RG transfor- 
mations. Since fixed point Hamiltonians are invariant under RG transformations, they do not have any length scale 
and can capture long-range physical properties which may survive even at the thermodynamic limit where the sys- 
tem size goes to the infinity. Then, fixed-point Hamiltonians help us to capture only the scale invariant properties of 
corresponding quantum phases. 

Now, since fixed point Hamiltonians are obtained by washing out short-range correlations, it is known that fixed 
point Hamiltonians in gapped quantum phases can be represented as frustration-free Hamiltonians where each term 
can be minimized locally without considering neighboring terms. Then, one may take a slight liberty and say that a 
study of all the possible frustration-free Hamiltonians is almost equivalent to a study of all the possible fixed point 
Hamiltonians which may appear as a result of RG transformations. In fact, one of the physical constraints (scale 
symmetries) which we shall impose on frustration-free Hamiltonians comes from the observation that fixed point 
Hamiltonians are free from any length scale (see discussions in Section[J!2)i. 



Frustration-free Hamiltonians: Since our goal is to study quantum phases arising in frustration-free Hamiltoni- 
ans, we now give their formal definition. 

When a Hamiltonian is said to be frustration-free, a ground state can be obtained by minimizing each local term 
in the Hamiltonian simultaneously. A useful feature of frustration-free Hamiltonians is that they can be represented 
as summations of projection operators. Consider a class of Hamiltonians which can be represented as summations of 
projectors: 

H^-J^Pj for all;, (6) 

J 

where Pj are projectors. Such Hamiltonians are called frustration-free when projectors commute with each other 
([Pj,Pj'] - 0) and whose ground states satisfy 

PM) = l-A). (7) 

A frustration-free Hamiltonian needs not to be a summation of projections. However, for any Hamiltonian whose 
ground states minimize each local term independently, one can represent the Hamiltonian as a summation of com- 
muting projectors. For example, one may obtain H(Q) and H{1) appeared in the discussions on the Ising model in a 
transverse field by setting Pj = ^{I + Z'-'-'Z*-'^'-') and Pj = ^(I + X*-'') up to some constant correction. 

2.2. Stabilizer codes and logical operators 

While frustration-free Hamiltonians may be used as representatives of quantum phases, their analyses are generally 
difficult both analytically and computationally [42 1. Fortunately, for stabilizer Hamiltonians which are quantum error- 
correcting codes realized by a certain subclass of frustration-free Hamiltonians, there exist various theoretical tools 
to analyze their ground state properties. Of particular importance are certain operators, called logical operators, 
which are central in studying coding properties of stabilizer codes. While properties of logical operators are discussed 
mainly in coding theoretical contexts, they become particularly useful in characterizing quantum phases arising in 
stabilizer Hamiltonians as we shall see in the main discussion of this paper. Here, we give a quick review of stabilizer 
codes and logical operators [26 ,| . We also review basic theoretical tools to analyze properties of entanglement through 
logical operators |32|. Some notations concerning stabilizer Hamiltonians are also introduced here. At the end of this 
subsection, an example of stabilizer Hamiltonians is presented. 

Stabilizer group: A stabilizer code can be characterized by an Abelian subgroup of Pauli operators. Consider an 
A^ qubit system and the Pauli operator group 

r = <//,x"\z(i>,...,x<^>,z*^>> (8) 

which is generated by Pauli operators X*^* and Z*^* acting on single qubits labeled by j. A stabilizer code is defined 
with the stabilizer group 

S = <5<",5^2\--- ,5<^-*>> c r, (9) 

which is a self-adjoint Abelian subgroup of the Pauli operator group V without containing -/ or // where k >Q. S *-'' 
represent independent generators for S. Elements in the stabilizer group .S are called stabilizers. Qubits are encoded 
in a subspace Vs spanned by states \ijj) which satisfy 

5W|i/r) = W) for all; (10) 

where Vs is called the codeword space of the stabilizer code. In total, k qubits can be encoded in V^, and encoded 
qubits are called logical qubits. 

System Hamiltonian: In the stabilizer formalism, there exists a system Hamiltonian which can support the en- 



coded states as degenerate ground states with a finite energy gap (Fig.lSta)). If a Hamiltonian 



H 



Z' 



(j) 



(11) 



satisfies {{S^-'\Vj}} = S, the ground state space of the Hamihonian is the same as the codeword space Vs since the 
energy of the system can be minimized for states satisfying 5''^'|(/') - \i//}. There are 2* degenerate ground states where 



k logical qubits can be encoded. We call Hamiltonians in Eq. (Ill stabilizer Hamiltonians . 

One may notice that stabilizer Hamiltonians are frustration-free up to some constant correction by setting Pj - 
5(7 + 5'*-''). Then, a ground state of a stabilizer Hamiltonian satisfy -Pj|i/') - \4') since S'^^^ij/) = |i/r), and thus, 
stabilizer Hamiltonians are frustration-free. Note that frustration-free Hamiltonians appeared in examples of QPTs in 
SectionlZTIare stabilizer Hamiltonians. 
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Figure 3: (a) The gapped energy spectrum of a stabilizer Hamiltonian witli degenerate ground states, (b) A bi-partition of a system into two 
complementary subsets of qubits A and B = A. Small dots represent qubits, and rectangles represent logical operators. 



Logical operators: Logical operators are Pauli operators ^ € !P which satisfy 

[A5'^>] = for all ; and t i {US). 



(12) 



Two logical operators I and I' are said to be equivalent if and only if ( and (' act in a similar way inside the ground 
state space: 



« (\ilj) = e'\ip) for all |(A) e Vs 

^ W e (S) 



(13) 
(14) 



where "~" represents the equivalence between logical operators { and /". Given a set of logical operators L, all the 
logical operators inside L are said to be independent when 



{{il 6 L)) n .S = 0. 



(15) 



Here, "0" represents a null set. 

Centralizer group: Logical operators can be found in the centralizer group which consists of all the Pauh opera- 
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tors commuting with stabilizers: 

C ^ {{Uer : [t/,5<^'] = 0, for all j ] ). (16) 

The centralizer group can be represented in a canonical form 



If, ... f, c(i) ... c(w-« \ 
C ^ I '' ' ^*' ^ ' ' ^ (17) 

\ ru ■■■ , rt, I 

where each operator represents independent generators for C. Pairs of generators £p and r^ anti-commute with each 
other while any other pair of generators commute with each other. For example, {{p, rp} = while [£p, r^] - for 
p + q. Pairs of anti-commuting generators {p and rp {p - 1 , ■ ■ ■ ,k) are independent logical operators in a stabilizer 
code. 

Canonical set of logical operators: A set of logical operators with the following commutation relations is called 
a canonical set of logical operators: 

Here, the commutation relations are represented in a way similar to a canonical representation where only the operators 
in the same column anti-commute with each other. Note that the choice of a canonical set of logical operators is not 
unique. With a canonical set of logical operators, the subspace Vs can be decomposed as a direct product of k 
subsystems: 

k 

l-A) = (^Kl0)p+j6p|i)p) (19) 

where fp and rp act non-trivially only on \Q}p and \l}p. One can choose the basis \Q)p and \a)p such that fp and rp act 
like Pauli operators applied to a logical qubit represented by |0)p and 11)^: 

£p\0}p = \l}p, ^p\l}p = |0)p, r^lO)^ = 10)^, rplDp = -|1)^. (20) 

Number of generators: The number of independent generators for a group of Pauli operators (9 e 'P is denoted 
as G{0). Note that G{0) does not count trivial generators such as / and //. Here, we give the numbers of generators 
for groups of Pauli operators which are central in discussions on stabilizer Hamiltonians: 

G{r) = 2N, G(S) ^ N-k, G(C) ^ N + k. (21) 

Code distance: The code distance measures the robustness of the stabilizer code: 

d = mm(w(U)) where U e C and U i {iI,S}. (22) 

Here, w{U) is the number of non-trivial Pauli operators in U. The code distance corresponds to a minimal number of 
single Pauli errors necessary to destroy a encoded qubit. Roughly speaking, a large code distance d for fixed A^ means 
that the code can encode qubits securely, and ground states of a stabilizer Hamiltonian are highly entangled. 

Projections: Properties of entanglement can be studied by considering a bi-partition of the entire system into two 
complementary subset of qubits A and B - A. In discussing the bi-partite entanglement, projections of Pauli operators 
are often used. The projection of a Pauli operator U e V onto a subset of qubits A is denoted as U\a. This keeps 
only the non-trivial Pauli operators which are inside A and truncates Pauli operators acting outside the subset A. The 
restriction of a group of Pauli operators O into some subset of qubits A is defined as 

Oa ^ {[UeO : U\^^l]). (23) 

Here, A represents a complement of A. Oa contains all the operators in O which are defined inside A. In this paper, 
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we use the restrictions of the stabiHzer group S, the centraHzer group C and the PauH operator group 'P, which are 
denoted as Sa, Ca and Pa respectively. 

Bi-partition of the system: For a stabiHzer code in a bi-partition into two complementary subsets of qubits 
(Fig.lSlb)), the following theorem holds 1321 . 

Theorem 1 (Bi-partitioning theorem of logical operators). For a stabilizer code with k logical qubits, let the number 
of independent logical operators defined inside a subset of qubits A be gA- Then, for two complementary subsets of 
qubits A and B — A, the following formula holds: 

gA+gB = 2k. (24) 

Note that A and B may support the same logical operators if they have two equivalent representations supported 
inside A and B respectively. Then, gA and gs may count the same logical operators twice. This formula becomes par- 
ticularly useful in determining sizes and geometric shapes of logical operators as we shall see in the main discussions 
of this paper. 

An example of a stabiUzer code: Here, we present an example of stabilizer codes to give some intuitions on 
definitions of logical operators and related groups of Pauli operators. The example we consider is called a five qubit 
code fi3ll44ll . Consider a system of five qubits which is governed by the following Hamiltoruan: 

5 

H = -YjS'-'\ ^W = xWy(./+i)yO+2);^o+3) (25) 

>=i 

where /' € Z5. Note that 5*-'' commute with each other, and this Hamiltonian is a stabilizer Hamiltonian. The stabilizer 
group is 

S = (5(i>,5(2),5<3\5(4') (26) 

since 5*''' is not independent from other stabilizers. This may be seen from the following equation: 

SW X 5(2> X 5(3> X 5(4' X 5<5' = /. (27) 

This stabilizer Hamiltonian has two degenerate ground states with ^ = 1 since G(S) - 4, and logical operators are 

( = X(')Z(2)X(3\ r = Z(')F(2)2(3) (28) 

where {I, r] - 0. Note that these are examples of logical operators, and there are many equivalent representations 
for logical operators. One may see that logical operators I and r commute with all the stabilizers 5*-'' through direct 
computations. The centrahze group can be represented in the following way: 

/ f r(l) c(2) c(3) c(4) \ 

Here, we denote the subset of qubits which consists of the first, second and third qubits as A. Then, we have gA - 2, 
gB - Q and gA+gB = 2 since logical operators £ and r can be defined inside A while there is no logical operator 
defined inside B. Since both I and r have non-trivial supports on three qubits, the code distance is t/ = 3. 

Quantum phases and logical operators: Having introduced stabilizer Hamiltonians and logical operators, let us 
discuss the connection between quantum phases and logical operators in stabilizer Hamiltonian. 

Quantum phases arising in correlated spin systems have been traditionally characterized by physical quantities, 
called order parameters, where different quantum phases are distinguished by different values of order parameters |[T| . 
Recently, it has been realized that entanglement, a non-local manifestation of the indistinguishability of quantum 
states, plays an essential role in characterizing quantum phases. In particular, various entanglement measures, such as 
entanglement entropies, have been used as order parameters and shown to be remarkably powerful in distinguishing 
quantum phases ll22l I45H 49 1 . With this connection in hand, we hope to find some quantities or objects which may 
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characterize entanglement in stabilizer Hamiltonians and serve as order parameters to distinguish quantum phases. 

Now, a basic idea of stabilizer codes is to encode qubits in highly entangled ground states so that small errors 
cannot break encoded qubits. Since the robustness of the code can be measured by the number of supports in logical 
operators, the existence of large logical operators for a fixed system size implies the presence of strong entanglement 
in ground states of stabilizer Hamiltonians f3T\ f32l . Then, geometric sizes of logical operators may be used as 
indicators of entanglement arising in ground states of stabilizer Hamiltonians. In the main discussion of this paper, 
we shall use geometric non-localities of logical operators as order parameters of quantum phases arising in stabilizer 
Hamiltonians. 

3. The model: Stabilizer codes with Translation and Scale symmetries (STS model) 

Now let us move on to the main problem of this paper, concerning possible quantum phases and their classifications 
in two-dimensional spin systems on a lattice. We wish to approach this problem by analyzing and classifying quantum 
phases arising in stabilizer Hamiltonians. 

In quantum information science, stabilizer codes have been particularly important since protecting or encoding 
qubits is essential in realizing various quantum information theoretical ideas such as fault-tolerant quantum computa- 
tion and quantum communication |2|. Recently, stabilizer codes have been increasing their importance in condensed 
matter physics too, for it has been realized that many interesting correlated spin systems may be described in the 
language of stabilizer codes Il3ll5l 12714321 . The connection between correlated spin systems and stabilizer codes is 
potentially useful and powerful since stabilizer codes are described through Pauli operators which obey an opera- 
tor algebra based on a finite group, and there is a possibility to analyze their properties analytically, or at least in a 
computationally efficient way |28 32|. 

However, there still remains a huge gap between realistic correlated spin systems commonly discussed in con- 
densed matter physics community and stabilizer codes commonly discussed in quantum information science com- 
munity. For example, most stabilizer codes, discovered from quantum coding theoretical motivations, encode qubits 
in ground states of highly non-local Hamiltonians. However, in real physical systems, Hamiltonians must have only 
local interaction terms. While there are a few pioneering works which have analyzed stabilizer Hamiltonians with 
local interaction terms 131113311341 . physically realizable Hamiltonians must have some physical symmetries such as 
translation symmetries too. 

Stabilizer codes supported by local and translation symmetric Hamiltonians may seem to have sufficient physical 
realizability. However, there also exists another important physical constraint which must be satisfied by stabilizer 
Hamiltonians in order for them to be viewed as representatives of quantum phases. Let us recall that fixed point 
Hamiltonians, which are frustration-free Hamiltonians obtained after RG transformations, must have some scale in- 
variant properties [16|. Then, we wish to analyze stabilizer Hamiltonians whose physical properties do not depend 
on the system size. 

Here, we propose a model of frustration-free Hamiltonians which are supported by not only local but also phys- 
ically symmetric Hamiltonians with scale invariance. In particular, our model is constrained to only the following 
three physical conditions: 

• Locality of interactions: A system of qubits, defined on a two-dimensional lattice, is governed by a stabilizer 
Hamiltonian which consists only of local interactions. 

• Translation symmetries: The Hamiltonian is invariant under some finite translations of qubits. 

• Scale symmetries: The number of degenerate ground states of the Hamiltonian does not depend on the system 
size. 

We call stabilizer Hamiltonians which satisfy the above three physical constraints Stabilizer codes with Translation 
and Scale symmetries (STS models). 

The main goal of this paper is to analyze quantum phases arising in STS models. In this section, we give the 
definition of STS models precisely, and develop a basic analysis tool of STS models. In Section [TT] we construct 



stabilizer codes realized by local and translation symmetric Hamiltonians. In Section 3.2 we discuss how the number 



of ground states changes as the system size varies, and impose scale symmetries on stabilizer Hamiltonians. In 
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Section 3.3 we give a formal definition of STS models. In Section 3.4 we introduce a basic analysis tool for STS 
models, which we shall call the translation equivalence of logical operators. Readers who want to start with specific 
examples of STS models may jump to Section|4]and Section|5] and return to this section later. 

3.1. Translation symmetries and locality of interactions 

Here, we construct stabilizer codes realized by local and translation symmetric Hamiltonians. Consider a system 
of qubits defined on a D-dimensional square lattice (hypercubic lattice) which consists of N = Li x ■ ■ ■ x L^, qubits 
where L,„ is the total number of qubits in the m direction for m - 1 , ■ ■ ■ ,D. We particularly consider stabilizer codes 
realized by system Hamiltonians with locally defined interaction terms S y. 



H^-Y,Sj (30) 



where each of S j is defined inside some geometrically localized regions. The meaning of "locally defined" will be 
clarified soon. The stabilizer group S is generated from each of the interaction terms S f S - {{S j, Vy)). 
In addition, we assume that the system Hamiltonian possesses translation symmetries: 

T„(H) ^ H (m= I,--- ,D) (31) 

where r„, represent translation operators which shift the positions of qubits by some positive integer v,„ in the m 
direction. For simplicity of discussion, we define: 

n,n = L,„/v„ (m=l,---,D) (32) 

which are integer values so that translations by T,,, can cover the entire lattice in a commensurate way. One may 
consider a hypercube of v'l x ■ ■ ■ x vo qubits as a unit cell of the system where v„, are periodicities of the lattice in the 
m directions. We also assume that the system has periodic boundary conditions for simplicity of discussion. Then, 
the entire space may be viewed as a square lattice (hypercubic lattice) defined on a D-dimensional torus: 

T^ = S'x---xS' (33) 

where SMs a circle. We may represent translation symmetries of the Hamiltonian in terms of the stabilizer group as 
follows: 

T„,{S) ^ S (m=l,---,D). (34) 

WTien a system is invariant under some finite translations, it is more convenient to treat a hypercube of 

V = vi X ■ ■ ■ X vb (35) 

qubits (a unit cell) as a single entity. One may view these v qubits as a single composite particle (Fig. |4]). Then, 
Tm is a translation operator with respect to composite particles, and the system Hamiltonian is invariant under unit 
translations of composite particles. The integer n,„ = L,„/v„j is viewed as the number of composite particles in the m 
direction. 

We may view this picture through composite particles as a coarse-graining of a system of qubits. When a system 
of qubits is coarse-grained through composite particles, a microscopic structure of qubits comprising each composite 
particle is lost, and one can neglect the effect of unitary transformations acting on qubits inside each composite 
particle. To make the difference between a picture through qubits and a picture through composite particles clear, we 
call a square lattice with composite particles a coarse-grained lattice. From now on, the entire system is considered 
as a square lattice o/«i X ■ ■ ■ xno composite particles. Thus, we treat each of composite particles as a single particle 
whose inner state corresponds to a subspace (C^)®^'. 

Having introduced a coarse-grained picture through composite particles, let us formally define the meaning of 
"locally defined" interaction terms. Each of interaction terms S j is defined inside some geometrically localized region 
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Figure 4: Stabilizer codes supported by local and translation symmetric Hamiltonians. A two-dimensional example is shown. Dotted lines represent 
the periodicities of the Hamiltonian. In the above example, the Hamiltonian is invaiiant under translations of three qubits in the 1 direction and 
translations of two qubits in the 2 direction. A system of qubits is coarse-grained by considering a unit cell of 3 X 2 qubits as a composite particle 
so that the Hamiltonian is invariant under unit translations of composite particles. Interaction terms Sj are defined locally inside a region with 2x2 
composite particles. 

of composite particles. Then, we can find a finite integer c <K n„, such that all the interaction terms S j are defined 
inside translations of a hypercubic region with cx- ■ -xc composite particles. For simplicity of discussion, we assume 
that all the interaction terms S j can be defined inside translations of a hypercubic region with 2 x ■ ■ ■ x 2 composite 
particles by assuming c -2 (Fig.|4]i|^ 

Translation symmetries and locality of interactions: Here, we define stabilizer codes supported by local and 
translation symmetric Hamiltonian in the following way: 

• The entire system is a square lattice (hypercubic lattice) of composite particles which consist of v qubits. 

• The system Hamiltonian H - - 2 5; is invariant under unit translations of composite particles: T„{H) = H for 
m = I,-- ,D. 

• Interaction terms S j are defined locally inside a hypercube of 2 x ■ ■ ■ x 2 composite particles. 

Though we have constructed stabilizer Hamiltonians on square lattices, our discussion in this paper can be readily 
applied to stabilizer Hamiltonians defined on any general lattices. 

3.2. Scale symmetries 

Now, we impose another physical symmetry, which we shall call scale symmetries. An important feature of fixed 
point Hamiltonians is that they do not have any length scale since they are invariant under RG transformations. In order 
for stabilizer Hamiltonians to be used as reference models for quantum phases, they must not have scale dependence 
too. So far, we have discussed the cases where the system size n = («i, ■ ■ ■ , no) is fixed. Here, we consider changes 
of the number of composite particles n^ while keeping interaction terms S y the same. In particular, we impose scale 
symmetries on translation symmetric stabilizer Hamiltonians by requiring that the number of degenerate ground states, 
or the number of logical qubits, does not depend on the system size n. 

Let us begin by representing the stabilizer group S for different system sizes ft. For this purpose, we label each 
composite particle according to its position on a coarse-grained lattice. Each composite particle is denoted as Pf 



^One might wonder if our constraint on the locality of interaction terms S j is too strict since it rules out stabilizer codes where S j are defined 
inside geometiically localized regions, but not inside regions with 2 X ■ ■ ■ X 2 composite particles. However, by increasing the size of composite 
particles, one can make all the interaction terms S / defined inside some regions with 2 X ■ ■ ■ X 2 composite particles. Of course, such newly 
defined composite particles may leave n,„ = L,„/v„ a non-integer value. However, it may be possible to analyze the original stabilizer code with a 
non-integer n,„ through the analysis of the newly defined stabilizer code with an integer n,„ . 
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where r — (ri , ■ ■ ■ , ro) represents the relative position of a composite particle with 1 < r,„ < n„, (see Fig. Bl. Now, 
let us denote a hypercubic region with 2 x ■ ■ ■ x 2 composite particles as Piocai which includes composite particles 
P, I ... ^ro with r,„ = 1 or 2 (Fig. l5]). We define the group of stabilizers which is generated from all the interaction terms 



S i defined inside Piocai as n: 



^ = {{Sj--Sj\p,„^.„, = /}). 



(36) 



Here, Piocai denotes a complement of Piocai- We call this group of stabilizers the local stabilizer group. Then, the 
stabilizer group S for the system size ft is generated from translations of the local stabilizer group n as follows: 



S = Sf, = {[t^' ■■■ T^^"in) : jc,„ = 1, ■ ■ ■ , «,„ for all m e Zo] ). 



(37) 



Here, T^" represents a translation of x„ composite particles in the m direction, and T„(7r) represents translations of 
elements in the local stabilizer group n. 
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Figure 5: Labeling of composite particles and the local stabilizer group n defined inside P/^ 



We may choose interaction terms S j defined inside P/ocai so that they generate the local stabilizer group ti. For 
example, let us pick up generators for n = {S '""'', S 2"'', ■ ■ ■)■ Then, the stabilizer code can be supported by a Hamil- 
tonian which consists of translations of 5'f™': 



« = -Z 






^s'; 



local \ 



(38) 



We note that some of interaction terms 5'°"'' may be defined inside regions smaller than a region Pi„cai- 

Now, we denote the number of logical qubits k for a stabilizer code for the system size n as kfj. More precisely, 
if we denote the number of independent generators for Sa as GiSa), we have kf, = N - G(Sfi) where N - v Hm «m- 
Here, v is the number of qubits inside each composite particle. Then, a stabilizer code has scale symmetries when the 
number of logical qubits kjj does not depend on the system size n: 






for all ft. 



(39) 



Thus, the number of logical qubits is always k regardless of the system size ft. Here, we emphasize that in a system 
with scale symmetries, the number of logical qubits k remains constant under not only global scale transformations: 
ft ^ eft where c is some positive integer, but also any changes of n,„. 

However, one may notice that the above constraint on scale symmetries is not well-defined for the case with n,„ - 1 
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Figure 6: (a) Effectively lowered dimensions. An example of a two-dimensional system with ni X H2 composite particles is shown. For the cases 
where n[ = I 01112 = I, the spatial dimension is effectively lowered to one. Interaction terms are defined inside some regions with 1x2 and 2x1 
composite particles respectively, (b)-(e) Examples of folded stabilizers in two dimensions, (b) ni,n2 i^ \- (a original stabilizer), (c) «i ^ 1 and 
n2 = 1. (d)ni = 1 and n2 ^ 1 ■ (e) " 1 = "2 = 1 ■ 

for some m since the local stabilizer group n is defined inside Piocai which is a hypercubic region with 2^ = 2 x ■ ■ ■ x 2 
composite particles (Fig. ^a)). This difficulty could prevent us from studying stabilizer codes at small ft to help us 
to study the properties of stabilizer codes at large ft. However, due to the commuting properties of interaction terms 
S '°™' in stabilizer Hamiltonians, one can define the local stabilizer group n properly even for the cases where n,„ - 1 
for some m. In order to show how this problem can be fixed, we consider a one-dimensional case (D - 1) first. Let us 
recall that Pi„cai consists of two composite particles Pi and P2 in one dimension. We can represent interaction terms 
^^"'^'''eTras 



g local ^ UjTi(Vj) 



(«i ^ 1) 



(40) 



where Uj and Vj are Pauli operators acting on a composite particle Pi. Then, Uj acts on a composite particle Pi 
while Ti(Vj) acts on a composite particle P2. Here, since the translation of 5'°™' must commute with 5'?™', we have 
[Uj,Vf] = for all j and /. Also, we have either [Uj,Uf] = [Vj,Vf] = or {Uj,Uj'} = {Vj,Vf} = since 



[S 



local c- local -\ 



\S , ] - 0. Now, for the case with rii - 1, we define S 



local 



in the following way: 



■t local 



= UjVj («i = 1) 



(41) 



by folding s'""'' at the boundary. Then, one can see that these folded ^'f™' commute with each other through some 
calculations. Thus, by this definition of S''"^"', the local stabilizer group n is well-defined for the case with ni - I. 

In a way similar to this, we can define interaction terms s'"™' e tt in higher-dimensions so that n is well-defined 
for any ft by folding S ^°"'' at the boundary in the m direction where n„, = 1 . For example, in two dimensions, let us 
consider the following interaction term (see Fig.l6tb)) 



^ local 



= Ui,iTi{U2,\)T2{Ui,2)TiT2{U2,2) «i > land«2 > 1 



(42) 



where t/i_i, f/i_2, t/2,1 and f/2,2 are Pauli operators acting on a composite particle Pij . Then, for the cases with n„, - 1 
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for some m, the interaction term is defined as follows (see Fig.l6tc)(d)(e)): 

S'"'"' ^ UuU2,iT2(UiaU2,2) «i = landn2>l (43) 

5'"™' = UijUi,2Ti(U2,iU2,2) m > landnz = 1 (44) 

giocai ^ Ui,iU2,iUi,2U2,2 rix = 1 and «2 = 1. (45) 

With this definition of interaction terms, one can discuss stabilizer codes for any n including the cases where 
n,„ - 1 for some m. As shown in Fig.l6|a), by setting «,„ = 1 for some m, the spatial dimension of the system can be 
effectively lowered. For example, a two-dimensional stabilizer code with n\ - I can be treated as a one-dimensional 
systemjj In such a case, the local stabilizer group, which is originally defined inside a region with 2x2 composite par- 
ticles, is now defined inside a region with 1x2 composite particles since interaction terms are folded in the 1 direction. 

Scale symmetries: Here, we summarize constraints resulting from scale symmetries: 

• The number of logical qubits k does not depend on the system size n and remains constant: 

kfj - k for all n. (46) 

Compared with the importance of translation symmetries, the importance of scale symmetries may be less obvious . 
In fact, there exist many examples of stabilizer Hamiltonians without scale symmetries. However, most examples 
without scale symmetries are trivial, as seen in the following example: 

H = -^z<'J)z('J+i> (47) 

where Z*''^' acts on a qubit labeled by (/, /) in a two-dimensional square lattice with periodic boundary conditions. The 
model is an array of one-dimensional classical ferromagnets which have interactions only in the vertical direction (the 
2 direction). The model has 2^' ground states where L\ is the number of qubits in the 1 direction. This model should 
be discussed as a one-dimensional system rather than a two-dimensional system since each chain of ferromagnets 
is completely decoupled. Thus, while we have imposed scale symmetries to limit our considerations to stabilizer 
Hamiltonians without length scales, scale symmetries help us to focus on physically interesting examples by excluding 
trivial stabilizer Hamiltonians such as the above example. 

There are also many stabilizer codes whose A;,^ remain finite {kf, < k) for all n, but do not remain constant. 
However, one can reduce such stabilizer Hamiltonians to STS models by increasing the size of composite particles 
so that stabilizer Hamiltonians possess scale symmetries with respect to newly defined composite particles. This 
restoration of scale symmetries is further discussed in [Appendix E| 

3.3. Formal definition of STS model 

Based on these discussions, let us formally define STS models. 

Definition 1. A stabilizer code defined with the stabilizer group Sa where n - {n\,- ■ ■ , no) is called a Stabilizer code 
with Translation and Scale symmetries (STS model) if and only if there exists a local stabilizer group n which satisfies 
the following conditions. 

• All the elements in n are defined inside a hypercubic region Piocai with 2^ = 2 x ■ ■ ■ x 2 composite particles. 

• The stabilizer group Sfj is generated from translations of the local stabilizer group n: 

S^ ^ {[ rf ■ ■ ■ T^^in) : jt;„ = 1, ■ ■ ■ , n„, for all m G Zd, } ) (48) 



When the system possesses symmetries, it is sometimes possible to lower spatial dimensions of the system effectively as a result of symmetries. 
See 1501 for general treatments on dimensional reductions in topologically ordered systems. 
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The number of logical qubits kfj is independent of the system size n: 

kjf - k for all n. 



(49) 



3.4. Translation equivalence of logical operators 

We have introduced a model of frustration-free Hamiltonians, an STS model, which may cover a large class of 
physically realizable stabilizer Hamiltonians. Now, let us develop a basic analysis tool for analyzing quantum phases 
in STS models, which we shall call the translation equivalence of logical operators. 

The analysis on logical operators is central in studying properties of entanglement arising in ground states of 
stabilizer Hamiltonians ll26l [311 [32l . Here, we wish to know the effect of translation and scale symmetries on logical 
operators. Translation symmetries of Hamiltonians are particularly useful in analyzing properties of logical operators. 
For example, a translation of a logical operator i is also a logical operator due to translation symmetries: 



^ e Ls 



T„,{() e U 



(50) 



where L^ is a set of all the logical operators. This observation may help us to determine possible forms of logical 
operators in STS models. 

However, the relation between the original logical operator ( and the translated logical operator T^it) is not 
immediately clear In particular, while ( and T^{€) have similar forms which can be transformed each other just 
by translations, { and T„,{€) might be not equivalent in general. For example, consider an array of one-dimensional 



classical ferromagnets discussed in Section 3.2 Then, we notice that unit translations of logical operators might not 
be equivalent to the original logical operators when logical operators are translated in the 1 direction. Thus, the need 
is to establish the relation between ( and its translation T„{€). 

In this subsection, we show that £ is always equivalent to its own translation T„ii€) as a result of scale symmetries 
in STS models. In particular, the following theorem holds. 

Theorem 2 (Translation equivalence of logical operators). For each and every logical operator i in an STS model, a 
unit translation of { with respect to composite particles in any direction is always equivalent to the original logical 
operator t: 

T„{{) ~ {, "it e Lft (m = I,--- ,D) (51) 

where L^ is a set of all the logical operators for an STS model defined with the stabilizer group Sj{. 

We call this property of logical operators under translations the translation equivalence of logical operators. 
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Figure 7: The translation equivalence of logical operators. Logical operators remain equivalent under any translations. Rectangles surrounding 
composite particles represent logical operators defined inside corresponding regions. Two-sided arrows mean that two connected logical operators 
are equivalent. 



The proof of Theorem|2]is presented in |Appendix A| Here, we only show the existence of finite translations which 
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keep logical operators equivalent for some sufficiently large n. In particular, we show the existence of a finite integer 
a„, which satisfies T"l"{{) ~ t for all the logical operators { for some fixed n. 

Let us consider translations of logical operators in the m direction. Due to the translation symmetries of the system 
Hamiltonian, translations of a given logical operator ( are also logical operators. Let us label each of translated logical 
operators as £{j) = T /„({): 

((0) -^ {(I) -^ ■■■ -^ E(n,„ - 1) ^ ^(0) (52) 

where "^" represents translations in the m direction. Here, due to the periodic boundary conditions, we have the same 
logical operator l{nm) = {(0) after translating n„, times. In principle, {(j) can be independent each other in the absence 
of scale symmetries. However, since there are at most 2^* diff'erent logical operators in this stabilizer Hamiltonian, 
there always exists a finite integers b„, < 2^'^ such that t{b„j) ~ {(0) where «„, is a multiple of b„,. By repeating this 
argument for all the 2k independent logical operators, one can find a finite integer a,„ such that f ~ T^'it) for all the 
logical operators (. 

The above argument shows that there exists some finite integer a,„ where T^"'{{) ~ I for STS models with suf- 
ficiently large n„,. However, one can obtain a much more general and stronger result than this, as summarized in 
Theorem l2] In other words, a,„ - 1 for any n„, in the above argument. The proof of Theorem l2] relies on the fact 
that the number of logical qubits k, or the number of logical operators, is not only small, but also independent of the 
system size «. 

The translation equivalence of logical operators is a key in analyzing quantum phases arising in STS models. In 
particular, since any translations of logical operators are equivalent to the original logical operators, in characterizing 
the ground state properties, the positions of logical operators are not important and only the geometric shapes of 
logical operators become essential as we shall see in the next section. The notion of topology arises also due to the 
translation equivalence of logical operators. 

4. One-dimensional STS model: quantum phases and logical operators 

Now, we start the analysis on quantum phases in STS models. In this section, we discuss quantum phases arising 
in one-dimensional STS models. 

Quantum phases arising in one-dimensional spin systems are relatively easy to analyze compared with higher- 
dimensional spin systems. In particular, there are several powerful numerical algorithms to analyze the ground state 
properties of parameterized Hamiltonians which connect frustration-free Hamiltonians. For example, the density ma- 
trix renormalization group (DMRG) approach, whose basic idea comes from RG transformations, provides efficient 
numerical algorithms to compute various physical quantities in one-dimensional systems 115111521 . Also, recent de- 
velopments on the matrix product state formalism solidified the usefulness of DMRG approaches and expanded its 
applicabilities |17, 18 1. In fact, it has been proven that any gapped spin systems on one-dimensional lattices can be 
efficiently simulated [53 1. Then, one might hope that a problem of finding and classifying quantum phases arising in 
one-dimensional STS models is not a difficult one. 

However, a classification of quantum phases is much more challenging than an analysis on a single parameterized 
Hamiltonian through a numerical simulation. In particular, since the existence of a QPT depends on paths of parame- 
terized Hamiltonians, one needs to analyze all the possible parameterized Hamiltonians connecting two STS models 
to see if they belong to different quantum phases or not. Thus, the need is to find order parameters to classify quantum 
phases in a path-independent way. 

A key idea behind RG transformations in characterizing quantum phases is the use of fixed point Hamiltonians 
which are invariant under RG transformations. These fixed point Hamiltonians capture only the scale invariant ground 
state properties of the corresponding quantum phases. Then, following the spirit of the analyses on quantum phases 
through RG transformations, we wish to classify quantum phases arising in STS models through some quantity or 
object which is scale invariant. 

A useful observation in addressing quantum phases arising in STS models is to realize that geometric shapes 
of logical operators do not change under scale transformations. Then, we hope that geometric shapes of logical 
operators may be used as "order parameters" to distinguish different quantum phases. In this section, we show that 
different quantum phases in one-dimensional STS models are completely characterized by geometric shapes of logical 
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operators. In particular, we show that two STS models belong to the same quantum phase if and only if geometric 
shapes of logical operators are the same. Thus, it is shown that two STS models with the same geometric shapes of 
logical operators can be connected without closing the energy gap, while two STS models with different geometric 
shapes of logical operators are always separated by a QPT. 



In Section 4. 1 we start our discussion by analyzing three specific examples of one-dimensional STSs. We discuss 



the role of geometric shapes of logical operators and their relation to the scale invariant ground state properties such 
as global entanglement. In Section 4.2 we extend our analysis to arbitrary one-dimensional STS models and give 



possible forms of logical operators. The ground state properties of STS models are also analyzed through geometric 



shapes of logical operators. Finally, in Section 4.3 we discuss the relation between logical operators and quantum 
phases, and show that quantum phases arising in one-dimensional STS models can be completely classified through 
geometric shapes of logical operators. 



4.1. Role of logical operators: concrete examples 

In this subsection, we analyze three specific examples of one-dimensional STS models in order to discuss the 
role of logical operators in classifying quantum phases. In particular, we discuss how the scale invariant ground state 
properties can be studied through geometric shapes of logical operators in these examples. 



In Section 4.1.1 we begin by analyzing a classical ferromagnet, which is the simplest example of a one-dimensional 
STS. In Section 4.1.2 we discuss an example without degenerate ground states or logical operators {k - 0) to compare 



physical properties in the presence and absence of logical operators. In Section 4.1.3 we analyze another example of 



a non-trivial one-dimensional STS which can be obtained by extending the five qubit code to a one-dimensional spin 
chain. 



4.1.1. Classical ferromagnet as a quantum code 

A classical ferromagnet, the simplest model of interacting spins, can be seen as a stabilizer code. Though the model 
has been completely analyzed a century ago, we utilize the simplicity of the model to give a concise, but insightful 
demonstration of the analysis on entanglement in ground states through geometric shapes of logical operators. 

Let us consider the following Hamiltonian (Fig.lHli: 



H 



.^zO-)zO-+i) 



(53) 



where Z*-'* represents the Pauli operator Z acts on a j-\h qubit. The total number of qubits is A^ and the system has 
periodic boundary conditions. The ground states of a classical ferromagnet are |0 ■ ■ ■ 0) and |1 ■ ■ ■ 1). 
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Figure 8: A classical feiTomagnet. 



The classical ferromagnet can be viewed as a stabilizer code with a translation symmetry since interaction terms 
Z*-''Z*-'^'' commute with each other The stabilizer group is 



^2(i)z(2),z(2)zW,...,Z<^>Z<i>). 



(54) 
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This stabilizer code is an STS model since it satisfies a scale symmetry. While there are A^ qubits and A^ stabilizers in 
this system, one stabilizer is redundant since 

Z(1)Z(2) X 2(2)^(3) X ... X Z^^^Z'^'^ = /. (55) 

Thus, the stabilizer group has only A^ - 1 independent generators, and A; = 1 for any A^ with G(S^) - N - I. Here, 
G{Sn) is the number of independent generators for the stabilizer group Sn- 
Logical operators of a classical ferromagnet are 

( = Z(i>, r = X(1'X(2)...X<^), K,r) = 0. (56) 

One can see that both of logical operators satisfy the translation equivalence of logical operators since T\ (r) = r and 
{Ti({) - Z^'-'Z*^' e Sn- According to geometric shapes of logical operators, we may call { a zero-dimensional logical 
operator and r a one -dimensional logical operator (Fig. [8]). Logical operators characterize transformations between 
degenerate ground states: 

^|0---0) = I0---0), f|l---l) = -|1---1> (57) 

r|0---0> = I1---1), r|l---l) = I0---0). (58) 

With the above classification of logical operators, we naturally think that a one-dimensional logical operator, 
non-locally defined all over the lattice, may characterize the existence of some global entanglement in ground states. 
Here, the ground state space is spanned by two orthogonal basis |0- ■ -0) and |1 ■ ■ ■ 1). Then, in principle, at zero 
temperature, a classical ferromagnet can support a GHZ state: |GHZ) = -)=(|0 ■ ■ ■ 0) + |1 ■ ■ ■ 1)). This GHZ state is 
"stabilized" by the one-dimensional logical operator r since r|GHZ) = |GHZ). Now, we discuss entanglement in a 
GHZ state in a relation with the one-dimensional logical operator r. A GHZ state is a globally entangled state since 
it is a superposition of two globally separated ground states |0 ■ ■ ■ 0) and 1 1 ■ ■ ■ 1) = X*'' ■ ■ ■ X*'^'|0 ■ ■ ■ 0). One can also 
quantify the global entanglement of a GHZ state by computing a mutual information E(A : B) lIU, which is a measure 
of entanglement between two subsets of qubits A and B: 

E{A : B) = E(A) + E{B) - E{A U B) (59) 

where E{A), E(B) and E(A U B) are entanglement entropies for regions A, B and AU B. Then, for a GHZ state, which 
is stabilized by the one-dimensional logical operator r, we always have E(A : B) = 1 for any pairs of disjoint regions 
A and B (Fig.lSll. Thus, global entanglement arising in a GHZ state is scale invariant since E{A : B) does not depend 
on the distance between A and B. 

4.1.2. Cluster state: a model without logical operators 

Next, let us discuss an example which does not have degenerate ground states or logical operators {k - 0). The 
example we discuss is called a cluster state, possessing short-range entanglement between neighboring qubits. Its two- 
dimensional generalization is particularly useful as a resource to realize quantum information theoretical ideas such as 
measurement based quantum computation ll54l . Though a cluster state has strong entanglement between neighboring 
qubits or composite particles, this short-range entanglement does not survive over the entire lattice, and is not a scale 
invariant property. In fact, such short-range entanglement can be removed by applying local unitary transformations 
on neighboring composite particles as we shall see below. 

The Hamiltonian for a one-dimensional cluster state is the following (Fig.l9]i: 

H = -2z<^'-'>X(»Z(>+'>. (60) 

i 

This model is a stabilizer Hamiltonian since interaction terms Z*^"''X*^'Z*^"^'^ commute with each other The model 
does not have degenerate ground states since all the interaction terms z'^-'^^'X^-'^Z*-'^'* are independent. The model has 
a unique ground state, called a cluster state, which satisfy the following conditions: Z*-'"'^X*-''Z^-'^'^|t/r) - \il/). 

A cluster state has short-range entanglement between neighboring qubits. While it is possible to see the existence 
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of such entanglement by writing down a cluster state as a superposition of product states explicitly, a representation 
would look ugly and give us few physical insights. 

A useful property of stabilizer Hamiltonians is that entanglement entropies of ground states can be computed 
easily through the stabilizer group S. Here, we begin by reviewing computations of entanglement entropies for 
ground states of stabilizer Hamiltonians. When a stabilizer Hamiltonian has only a single ground state, the ground 
state can be represented in terms of the stabilizer group as follows: 



I^WI - ^U^^ + Sj) 



(61) 



SjeS 



where A^ is the total number of qubits, and S j are independent generators for the stabilizer group S. Then, the density 
matrix of a ground state of a stabilizer Hamiltonian can be represented in the following way: 



^« = ^n^^^^^) 



(62) 



SjESr 



where S j are independent generators for a restriction of the stabilizer group into R which is denoted as Sr. G(Sr) is 
the number of independent generators for Sr. Since the entanglement entropy is defined as 



Er = Ti\pR\ogpR], 
it can be represented in terms of the restriction of stabilizer group in the following way fSS^: 

Er = Vr-G(Sr), 



(63) 



(64) 



where Vr is the number of qubits inside R. 

Now that we have a formula to compute entanglement entropies, let us compute the entanglement entropy for a 
region Ai which consists of only one qubit. Then, we have E/^^ - 1 since V^j = 1 and G(Sa,) = 0. Next, let us 
compute the entanglement entropy for a region Aj which consists of j consecutive qubits. Then, we have £4, = 2 for 
1 < j < N - I since Va - j and G{Sa) - j -2. This indicates the existence of short-range entanglement between 
neighboring qubits. 

However, a cluster state does not have global entanglement. This can be seen from the fact that E{A : B) - for 
any disjoint regions A and B. Thus, entanglement arising in a cluster state is short-range and is not scale invariant. As 
we shall soon see, this is a direct consequence of the absence of one-dimensional logical operators. 
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Figure 9: A cluster state. A mapping to an STS model is shown. 



Reduction to an STS model: This model can discussed in the framework of an STS model by considering two 

even. Let zf and xf 
, Nil) particle with the 



consecutive qubits as a single composite particle (v — 2) when the total number of qubits A^ is even. Let Zp and xf 



be Pauli operators acting on a p-\h qubit (p — 1,2) inside a j-Xh composite particle {j - 1, 
following commutation relation: 



X 



7(fl 
(J) 






(65) 
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Then, after applying some appropriate unitary transformations on qubits inside each composite particle, one can 
represent the Hamiltonian for a cluster state in the following way: 



H 



■Y/^7^^^. 



yt^'^t''. 



(66) 



Disentangling short-range entanglement: Even after a coarse-graining, the model still possesses short-range 
entanglement between composite particles. However, there is no global entanglement since E(A : B) - Q for any 
disjoint regions A and B of composite particles. One may also see why there is no global entanglement by considering 



how stabilizers z!-'''z'-''^" and xi^'x! 
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stabilizers 



act on qubits inside each composite particle. As shown in Fig. 
connect only the first qubits in j-th composite particles and the second qubits in j + 1-th composite particles. Thus, 
there is no entanglement which survives over the lattice, and entanglement is established only between neighboring 
composite particles. 

As long as global entanglement is concerned, a cluster state is similar to a product state since both states have 
E(A : B) -0 for any disjoint regions A and B. In fact, one can disentangle neighboring entanglement between qubits, 
or between composite particles by applying local unitary transformations and reduce a cluster state to a product state. 
Here, we first look at how short-range entanglement can be disentangled in a system of qubits. Then, let us discuss 
how neighboring entanglement can be disentangled in a coarse-grained lattice. 
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Figure 10: A disentangling operation in a cluster state described through composite particles. Short-range entanglement between composite 
particles can be removed by some local unitary transformations acting on neighboring composite particles. 

The Hamiltonian for a cluster state can be obtained by applying control-Z operations on the Hamiltonian for a 
product state H — - YijX'^^\ Here, the control-Z operation acting on two qubits is defined as follows: 



CZ : 100) ^ 100), 101) ^ 101), |10) ^ |10), |11) ^ -|11) 



(67) 



where CZ represents the control-Z operation. The effect of the control-Z operation on two qubits can be represented 
in terms of Pauli operators in the following way: 



CZ 



Zi, Z2 

^1. ^2 



CZ 



Zi, 

^lZ2, 



Z2 

Z1X2 



(68) 



where the operation transforms Xi into X1Z2 and X2 into Z1X2. Here, Z\ and X\ act on the first qubit, and Z2 and X2 
act on the second qubit. CZ ' represents the inverse of CZ. 

Now, we describe the reduction from a cluster state to a product state by representing control-Z operations in 
terms of composite particles. In particular, through some observations, we notice that the following local unitary 
transformations U^J^ acting on j-\h and 7 -I- 1-th composite particles can disentangle neighboring entanglement between 
composite particles in Eq. ([66|: 
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where Pauli operators are transformed each other through U'-^K Here, we note that {/'^' commute with each other: 
[U^^\ U^^ '] - 0. Then, by applying a unitary transformation U - Ylj U^^\ one can transform the original Hamiltonian 
in Eq. ([66]l to 



UHU- 



•Z^l"' 



Z<' 



(70) 



and short-range entanglement between neighboring composite particles can be washed out (Fig. 10 1. Thus, the model 
can be reduced to a Hamiltonian which supports a "product state of composite particles". 

These observations imply that short-range entanglement arising in a cluster state is not scale invariant. According 
to the basic principle of a classification of quantum phases through RG transformations, such short-range entanglement 
is not relevant for characterizations of quantum phases. At the end of this section, we shall confirm this expectation 
rigorously for quantum phases arising in STS models. 



4.1.3. Extended five qubit code: reduction to a classical ferromagnet 

We discuss our final example of one-dimensional STS models. We consider an example which can be obtained 
by generalizing the five qubit code to a one-dimensional spin chain. By introducing composite particles, one can 
discuss the model in the framework of STS models. We see that geometric shapes of logical operators are the same 
as those in a classical ferromagnet, and physical properties are similar. Then, we show that the model can be reduced 
to a "classical ferromagnet of composite particles" by disentangling short-range entanglement between neighboring 
composite particles. 



As we have seen in Section 2.2 the five qubit code is constructed on five qubits, which is defined with the following 
stabilizer generator: x'^^^Y'^^^Y^^X'-'*^ and its translations. Inspired by this five qubit code, we consider the following 
Hamiltonian defined with A^ qubits: 






(71) 



One can easily see that this extended five qubit code is also a stabilizer code since all the interaction terms x^-''F^'^'^y'-'^^^X'^^^^' 
commute with each other. The code satisfies a scale symmetry since 



pr T^U)yO+l)yO+2)j^O+3) _ j 



(72) 



and ^ = 1 for any A^. 
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Figure 1 1 : The extended five qubit code. 
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This model can be reduced to an STS model by considering three consecutive qubits as a single composite particle 
- 3). Let Zp and Xp be Pauli operators acting on a 
(j - 1, ■ ■ ■ , n) with the following commutation relations: 



(v - 3). Let Zp and Xp be Pauli operators acting on a p-th qubits (p - 1, ■ ■ ■ ,3) inside a j-th composite particle 



^U) 2^7) ^U) •\ 

x'j\ ^/, 4'' r ^^^^ 

Then, after applying some appropriate local unitary transformations on qubits inside each composite particle, one can 
represent each of stabilizers in the following way: 

Z;^'>zf", x;-'''4-''+", Z'^^ZI^^'K (74) 

With this notation through composite particles, one can easily write down logical operators: 

i = Zf\ r = Xf^xf^---Xf\ {£,r} = 0. (75) 

We notice that geometric shapes of logical operators are the same as those of a classical ferromagnet. Then, we expect 
that physical properties of the model are similar to those of a classical ferromagnet. In fact, one can discuss properties 
of entanglement in a way similar to a classical ferromagnet. The existence of a one-dimensional logical operator r 
implies the presence of a global entanglement. Indeed, for a ground state satisfying r\tp} = |i/'), we have E{A : B) = 1 
for any disjoint regions A and B, which can be computed easily by using a formula in Eq. (|64]). Here, we choose two 
degenerate ground states of the model in the following way: 

mo) = liAo), m) = -I^Ai), r\t/,o} = liffi), r\if,i} = \ifro). (76) 

Then, we notice that li/r) has a GHZ-like entanglement since \4f) - -L(|j^q) + |(/ri)) and |i/') is a superposition of two 
globally separated ground states. 

As long as the scale invariant ground state properties are concerned, an extended five qubit code is similar to a 
classical ferromagnet. In fact, this model can be reduced to a classical ferromagnet by applying unitary transformations 
acting on neighboring composite particles in a way similar to the reduction of a cluster state: 

f/W : Z^^z!:i^'^ ^ Z'f, X</'xf' > ^ Xf ". (77) 

Note that {U^-'\ U^-' ^] = 0. Now, after this unitary transformations, only the third qubits in each composite particle are 
correlated through stabilizers Z^ Z^'^ while the first and second qubits in each composite particle are decoupled due 
to the existence of Z^ and Xj . Thus, the extended five-qubit code is equivalent to a classical ferromagnet, when 
decoupled qubits are removed. Later, we shall show that this model and a classical ferromagnet belong to the same 
quantum phase. 

Connection with RG algorithms based on the tensor product state formalism: We have seen that an extended 
five qubit code can be reduced to a classical ferromagnet by applying disentangling operations on neighboring com- 
posite particles. Here, it may be worth noting that these disentangling operations are closely related to a "disentangler" 
which is a key element in a novel RG algorithm based on the tensor product state formalism fT5l [T9l . 

In the previous section, we have coarse-grained the system of qubits by introducing composite particles where 
microscopic properties of each qubit inside a newly defined composite particle are completely lost. In a coarse-grained 
lattice, one may consider two STS models as the same when they can be transformed each other through unitary 
transformations acting on qubits inside each composite particle. For example, Hamiltonians H - - 2jZ'-'*Z*-'^'* and 



2.1 



H - - J^yX'^'X*^^'' may be considered to be the same, as discussed in Section : 

However, coarse-graining is not sufficient to characterize the scale invariant ground state properties since there 
still remain some short-range correlations which cannot be washed out by coarse-graining, as we have seen in anal- 
yses on a cluster state and an extended five qubit code. Such short-range entanglement must be also removed in 
classifying quantum phases, according to the basic philosophy of the classification of quantum phases through RG 
transformations. 

Recently, a remarkable idea of removing short-range correlations has been proposed in a search for efficient RG 
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algorithms based on the tensor product state formaUsm. The key idea is the use of a disentanghng operation, called 
a disentangler, through local unitary transformations to remove these short-range entanglement between neighboring 
particles |,19J . It has been demonstrated that RG transformations which combine coarse-graining and disentangling 
operations work in a remarkably efficient way in analyzing quantum phases arising in two-dimensional strongly cor- 
related spin systems. 

A disentangling operation used in the analysis on a cluster state is essentially similar to a disentangler used in this 
RG algorithm. In particular, a central idea behind these disentangling operation is that short-range entanglement is 
irrelevant to characterizations of quantum phases. This observation will be rigorously confirmed for quantum phases 
arising in STS model at the end of this section. 

Based on observations obtained in the analyses on a cluster state and an extended five qubit code, we wish to 
consider two STS models as the same when they can be transformed each other by local unitary transformations. By 
local unitary transformations, we include the following two elements: 



Unitary operations on composite particles: 

particle. 



Unitary transformations acting on qubits inside each composite 



• Disentangling operations: Unitary transformations acting on neighboring composite particles. 

Here, we note that this is consistent with our original approach to distinguish quantum phases through geometric 
shapes of logical operators since geometric shapes of logical operators are invariant under local unitary transforma- 
tions. Later, we shall see that one-dimensional STS models connected through local unitary transformations belong 
to the same quantum phase jj 

4.2. Logical operators in one-dimensional STS models 

We have seen that geometric shapes of logical operators are central in analyzing the scale invariant ground state 
properties of STS models through several examples. In this subsection, the analysis is extended to arbitrary one- 
dimensional STS models. We obtain a canonical set of logical operators (all the independent logical operators) of 
arbitrary STS models, and see that logical operators are either zero-dimensional or one-dimensional, as in a classical 
ferromagnet. We discuss how the ground state properties, such as global entanglement, in STS models can be studied 
by geometric shapes of logical operators, and show that any one-dimensional STS models can be reduced to multiple 
copies of a classical ferromagnet, or a product state, by disentangling neighboring entanglement between composite 
particles through local unitary transformations. 
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Figure 12: Stabilizers and logical operators in arbitrary STS models. 

Let us consider an STS model defined with composite particles which consist of v qubits. For simplicity of 
discussion, we neglect stabilizers acting on single composite particles since decoupled qubits are to be removed. Then, 

'^In general, local unitary transformations are used in much broader sense. In particular, any unitary evolutions induced by time evolution of 
local Hamiltonians for a finite time may be called local unitary transformations 1151 . 
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stabilizers in one-dimensional STS models can be represented in the following way (Fig. 12 1, as shown in Appendix 
11 

Ferromagnetic pait: z|^''z|^'""\ ■ ■ ■ , Z^^'^Z^^'""* (78) 

Short-range entanglement: Z^f jZ'^^j"' ^uXkti^ ^ ' ' ' ' z|;i\z[,-''"", x'^\x\!^'^'^ (79) 

where v - fc is an even integer. Here, stabilizers zfz\^^ , ■ ■ ■ ,Z'^ Z^^ create ferromagnet-like correlations, while 

'^kl\^k+2 ' ^1+1-^1+2 ' ■ ■ ■ ' -^y i-^i' ' ^,1 i^v' ^ create short-range entanglement between neighboring composite par- 
ticles as in a cluster state. Then, logical operators are 

n(^„) = <! ^1' ;;;' ^* ^ (so) 



where 



= Z(», rj^X^pxf---X^p, forp = l,...,L (81) 



"p - ^p 



Thus, geometric shapes of logical operators in one-dimensional STS models are either zero-dimensional or one- 
dimensional. Zero-dimensional logical operators and one-dimensional logical operators always form anti-commuting 



pairs (Fig. 13 i 
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Figure 13: Reduction of one-dimensional STS models to classical fen'omagnets by disentangling short-range entanglement between neighboring 
composite particles. Geometric shapes of logical operators characterize scale invariant propeilies of STS models. 

We notice that geometric shapes of logical operator are the same as those in a classical ferromagnet. Then, we 
expect that the scale invariant ground state properties in one-dimensional STS models can be discussed in a way 
similar to discussions on a classical ferromagnet. This expectation turns out to be true. In fact, due to the existence 
of one-dimensional logical operators, ground states of STS models can have GHZ-like entanglement, as in a classical 
ferromagnet. For example, when the number of logical qubits is k, there exists a ground state which has a mutual 
information E{A : B) - k for any disjoint regions of composite particles A and B. Note that this can be confirmed 



easily by using a formula in Eq. (64 1. 



Even more, one can show that all the one-dimensional STS models can be reduced to classical ferromagnets by 
applying disentangling operations between neighboring composite particles. Let us focus on a pair of stabilizers 
^1+1^1+2 ^^^ ^k+\^k+2 ■ Th^ii' these two stabilizers can be transformed into Z^-^j and Xj^^^ by removing short-range 
entanglement between neighboring composite particles, as in the discussion of a cluster state. One can repeat the 
similar arguments for all the other pairs of stabilizers. Thus, one-dimensional STS models with k logical qubits can 
be reduced to the following STS model with v - k through local unitary transformations: 
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This model consists of "multiple copies" of a classical ferromagnet embedded in a non-interacting way. The 
/7-th classical ferromagnet connects p-th qubits inside each composite particle with Zp Zp , and different classical 
ferromagnets embedded in the model are decoupled from each other. Therefore, when two STS models have logical 
operators with the same geometric shapes, they have similar global entanglement in ground states and can be reduced 
to an STS model in the above form. 

Here, we mention the importance of geometric shapes of logical operators as indicators of global entanglement in 
ground states. Since geometric shapes of logical operators are invariant under local unitary transformations, they can 
capture the scale invariant ground state properties such as global entanglement as seen in a GHZ state (Fig [T3|). We 
shall soon see that geometric shapes of logical operators distinguish quantum phases in one-dimensional STS models 
completely, serving as order parameters. 

4.3. Quantum phases and geometric shapes of logical operators 

We have seen that the ground state properties, such as global entanglement, in STS models are similar when 
geometric shapes of their logical operators are the same. Here, our hope is to use geometric shapes of logical operators 
as order parameters to distinguish quantum phases. However, there is still an important gap between geometric shapes 
of logical operators and the notion of quantum phases. In particular, it is not clear if different geometric shapes of 
logical operators lead to different quantum phases separated by a QPT Also, it is not clear if the same geometric 
shapes of logical operator imply that two STS models belong to the same quantum phase or not. 

Here, we establish the relation between quantum phases and geometric shapes of logical operators. In this sub- 
section, we show that two STS models belong to different quantum phases if and only if geometric shapes of logical 
operators are different. In particular, we show that two STS models with the same geometric shapes of logical op- 
erators can be connected without closing the energy gap, while two STS models with different geometric shapes of 
logical operators are always separated by a QPT. 

Cases with the same number of logical operators: Let us begin by showing that two STS models belong to the 
same quantum phase when they have the same geometric shapes of logical operators. 

Consider two STS models Ha and Hg with the same number of logical qubits k, or the same number of anti- 
commuting pairs of logical operators. For simplicity of discussion, let us assume that the number of qubits inside 
each composite particle is v for both H^ and H^ iy > k). We also assume that both H^ and H^ have the same system 
size. We denote the projection operators onto the ground state spaces of Ha and Hg as Pa and Pb- From the discussion 
in the previous subsection, there always exist a local unitary transformation U which transform Pa into Pb'- 

UPaU-' = Pb. (83) 

Then, one can show that there exists a parameterized Hamiltonian which connects Ha and Hb without closing the 
energy gap. 

Here, we give a sketch of a proof for the existence of such a parameterized Hamiltonian H{e). From the discussion 
in the previous section, unitary transformations U which connects Pa and Pb can be decomposed into two parts: 

^ ^composite ^ ^disentangle WnCrC V^ composites ^disentangle! ~ ^- v^'^J 

The first part U composite represents unitary transformations acting on qubits inside each composite particle, while the 
second part Udisentangie represents disentangling operations acting on neighboring composite particles. Let us begin 
with the case where U disentamiie - I- Then, by gradually inducing unitary transformations on each composite particle, 
one can transform Ha into UHaU^^ without closing the energy gap or changing the number of ground states. Now, 
since the projection into the ground state space of UHaU^^ is the same as Pb, one can change UHaU^ to Hb without 
closing the energy gap. Thus, Ha and Hb belong to the same quantum phase. Next, let us consider the case where 
U composite - I- Let US recall a transformation from a product state Hamiltonian Ha - -2;^*''' to a cluster state 
Hamiltonian Hb - - J]jZ*''"''X*^'Z*''^''. Since the control-Z operation acting on y-th and j + 1-th qubits can be 
written as 



CZ = exp 



/-(Z*-'') - /)(Z('+'> - /) 
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(85) 



by gradually inducing the control-Z operation such that 



CZ(e) - exp 



/^(ZW-/)(Z(^'+1>-/) 



(86) 



from e = to e = 1, one can transform a product state Hamiltonian to a cluster state Hamiltonian without changing 
the energy gap by setting //(e) - CZ(e)///i(CZ(e)) '. By generalizing this idea, one can show that Ha and Hg can 
be transformed each other without closing the energy gap when the number of logical operators is the same. It is 
straightforward to extend these discussions to the cases where Ucomposiie + I and Udisemangie + /■ Also, when the 
numbers of qubits inside each composite particles are different for //^ and //^ (say, v^ and vg with v^ + vg), we 
coarse-grain by grouping v qubits into a composite particle where v is the least common multiple of va and vb- Then, 
one can repeat the similar argument and show that Ha and Hb belong to the same quantum phase. This completes the 
sketch of the proof. 

Cases with the different numbers of logical operators: Next, let us show that quantum phases represented by 
two STS models are different when geometric shapes of logical operators are different. Since zero-dimensional and 
one-dimensional logical operators always form anti-commuting pairs, different geometric shapes of logical operators 
imply different numbers of degenerate ground states. Then, when connecting two STS models through a parameterized 
Hamiltonian, the energy gap must close at some point as excited states need to be ground states in the course of a 
parameter change. Thus, two STS models with different k belong to different quantum phases which are separated by 
aQPT. 

Here, we would like to make some comments on phase transitions between two STS models Ha and Hb which 
commute with each other, but have different numbers of logical operators. In such cases, phase transitions may occur 
exactly at e = or e = 1, instead of some intermediate point between e = and e - I. Let us look at an example. 
Consider the Ising model in a parallel field: 



AJ) 



H(e) = -(1 - e) ^ zWz(-'+i> - e Ji ^' 
j J 

Ha =//(0) = -^Z<-''Z<-'+i\ Hb = //(I) ^-^i^' 



(87) 
(88) 



Since Ha and Hb belong to different quantum phases with different numbers of logical operators, we expect that the 
model undergoes some phase transition. This parameterized Hamiltonian is exactly solvable for any e since all the 
terms in //(e) commute with each other for any e. Though this model connects Ha and Hb from e = to e = 1, the 
phase transition occurs at e = 0. To see this more clearly, let us extend our analysis to the cases where e < too. At 
e = 0, the model has degenerate ground states |0 ■ ■ ■ 0) and |1 ■ ■ ■ 1). At e > 0, the model has a single ground state 
|0 ■ ■ ■ 0). At e < 0, the model has a single ground state 1 1 ■ ■ ■ 1). Since the ground state properties change drastically at 
e = 0, the model undergoes a phase transition at e = 0. In this model Ha and Hb may not be "separated" by a phase 
transition since Ha lies at the transition point. However, we consider that Ha and Hb belong to different quantum 
phases since the ground state degeneracy is lifted for e ^ 0. 
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Figure 14: A "phase diagram" of one-dimensional STS models. Different quantum phases can be characterized by geometric shapes of logical 
operators (the number of logical operators). Different quantum phases are separated by QPTs. 
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Summary and applications: We summarize the main result of this section (Fig. 14 1. 



• Quantum phases in one-dimensional STS models can be characterized by geometric shapes of logical operators. 
Quantum phases represented by two STS models are different if and only if the numbers of logical operators, 
or the numbers of logical qubits k, are different. 

Here, we mention the importance of the translation equivalence (Theorem l2]i in this classification of quantum 
phases in one-dimensional STS models. The underlying reason why we can specify quantum phases only through 
geometric shapes of logical operators is due to the translation equivalence of logical operators. In particular, as a 
result of the translation equivalence of logical operators, any translations of a logical operators are equivalent to 
the original logical operators. Then, one can distinguish quantum phases only through geometric shapes without 
considering the positions where logical operators are defined. 

Now, let us look at some examples. Within the framework of STS models, one can classify frustration-free 
Hamiltonians appeared in Sectionl2]in the following way: 

Ha = -2z(^>Z(-'+'> ~ //; = -^xWx(-''+'^ oo Hb = -Z^^'' (89) 

; J J 

where Ha and Ha' belong to the same quantum phase while Hb belong to a different quantum phase. This classification 
is consistent with the fact that the Ising model in a transverse field, 

H(e) = -(1 - e) 2 Z(^>z(-'+'> - e^ ^<^', (90) 

; j 

undergoes a QPT. Also, Ha - - 2;Z'-''z'^^'' and H'^ - - ^^X'-'^X'-'^'' belong to the same quantum phase since they 
can be transformed each other through single qubit rotations without closing the energy gap. In combining the models 
discussed in Section |4T[ we have the following classification: 



Ha = -^Z^Z'^+i^ - H'^ - -^X^»X^>+'> ~ H2 = -^xWy(>+i)FO+2)^0+3) 
j j J 

^ Hb ^ -2^*'' ~ ^B = -2z('-'>X(»Z(>+'>. (91) 

j j 

While our discussions show only the existence of QPTs between STS models, analyses on QPTs in specific pa- 
rameterized Hamiltonians are also important problems. However, detailed theoretical analyses on each parameterized 
Hamiltonian are beyond the scope of this paper. We note that QPTs occurring in some parameterized Hamiltonians 
connecting the above STS models are studied in previous works Il23ll25l . 

Connection with tlie symmetry-breaking theory: We have seen that changes of geometric shapes of logical op- 
erators may characterize different quantum phases. Here, we make comments on the similarity between our approach 
through geometric shapes of logical operators and the Landau's symmetry breaking theory on QPTs. 

The key idea of the Landau's symmetry breaking theory is to explain QPTs through changes of symmetries. Let 



us consider an example, the Ising model in a transverse field, described in Eq. (90i. The parameterized Hamiltonian 
has a symmetry with respect to the following global operator r. 

r = X"'---X*^\ rHr^ = H. (92) 

Let us see how the symmetry with respect to r changes as e varies. For e = 1, the ground state is \if/o} - | + ■ ■ ■ +). 
Then, the symmetry with respect to r is said to be spontaneously broken since the expectation value of r is one and 
fixed to a single value. On the other hand, for e = 0, the ground states are li/^o) - |0 ■ ■ ■ 0) and \iffi) = |1 ■ ■ ■ 1). Then, 
the symmetry with respect to r is not broken since the expectation value of r inside the ground state space are not 
fixed to a single value. Thus, the symmetry with respect to a global operator r changes during a QPT, which implies 
the existence of a QPT in this model. 

Now, let us recall that logical operators are Pauli operators which commute with the system Hamiltonian, and 
(H(^^ - H where // is a stabilizer Hamiltonian. Then, logical operators ( characterize symmetries of stabilizer 
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Hamiltonians. Here, one may notice that r is a one-dimensional logical operator of a classical ferromagnet Ha = H{Q). 
Since r is a logical operator, r has an anti-commuting pair t - Z*'\ Then, the expectation value of r is not fixed inside 
the ground state space of Ha, and the symmetry with respect to r is not broken for Ha- However, this symmetry is 
broken for Hg since Hs does not have any logical operator. Thus, our approach in classifying quantum phases through 
logical operators is consistent with the Landau's symmetry breaking theory on QPTs. 

5. Two-dimensional STS model: topological quantum phases and geometric shapes of logical operators 

In one-dimensional STS models, different quantum phases are completely characterized by geometric shapes of 
logical operators. However, the only possible geometric shapes are anti-commuting pairs of zero-dimensional and 
one-dimensional logical operators. As a result, one-dimensional quantum phases are classified through only the 
number of logical operators without actually considering geometric shapes of logical operators. 

Unlike one-dimensional systems, two-dimensional systems have rich varieties of quantum phases, which may 
be seen from rich varieties of possible geometric shapes of logical operators in two dimensions. Here, we wish to 
search for possible quantum phases in two-dimensional STS models by finding logical operators and classify different 
quantum phases through geometric shapes of logical operators. 

Topological order and logical operators: To begin with, in order to discuss quantum phases through geometric 
shapes of logical operators, we need to study the ground state properties of STS models and their relations to logical 
operators. What makes two-dimensional systems strikingly distinct from one-dimensional systems is the possibility 
of the existence of topological order Topological order is a highly non-local quantum correlation which is stable 
against any local and small perturbations. This peculiar correlation is known to be "topological", meaning that it is 
a scale invariant physical property which depends on topological structures of geometric manifolds where systems 
are defined [37|. As we shall soon see, some of two-dimensional STS models, including the Toric code and the 
topological color code l'55'l, have topological order. Since topological order is a scale invariant characterization 
of non-local correlations arising in quantum many-body systems, in classifying quantum phases with topological 
order (topological phases), systems with "different kinds of topological order" must be appropriately distinguished. 
Here, we wish to distinguish topological phases in two-dimensional STS models through geometric shapes of logical 
operators. 

The main puzzle in analyzing topological order is the fact that topological order in correlated spin systems is 
currently characterized in various methods, such as the existence of anyonic excitations |i3]|5]|6l|56l|57l, topological 
entanglement entropy P6ll47l and the absence of local order parameters lfT4l l37, 42 1. While these characterizations 
manifestly captures topological properties of non-local correlations, they must be some aspects of what is called 
"topological order". Though the connection between anyonic excitations and topological entanglement entropy has 
been established through the framework of topological quantum field theory B6l . topological entanglement entropy 
cannot characterize anyonic excitations completely lf58l . Therefore, if we hope to use geometric shapes of logical 
operators as order parameters for systems with topological order, all these existing characterizations of topological 
order must be explained through geometric shapes of logical operators in an unified way. 

We start our analyses on quantum phases in two-dimensional STS models by characterizing topological order 
arising in STS models completely through geometric shapes of logical operators. In Section [5T| we study three 
examples of two-dimensional STS models, and analyze the connection between topological order and geometric 
shapes of logical operators. In particular, we approach three different characterizations of topological order, the 
existence of anyonic excitations, topological entanglement entropy and the absence of local order parameters, through 
geometric shapes of logical operators. In Section |5.2| we extend our analysis to arbitrary two-dimensional STS 
models and describe possible geometric shapes of logical operators. Based on geometric shapes of logical operators, 
topological order arising in two-dimensional STS models are discussed. Then, we show that geometric shapes of 
logical operators can distinguish topological properties of two-dimensional STS models appropriately. 

Quantum phases and logical operators: Then, we shall move on to the analyses on quantum phases in two- 
dimensional STS models and their classifications. Since a topological STS model and a non-topological STS model 
have totally different physical properties, one may naturally think that they belong to different quantum phases. This 
expectation is also consistent with the basic spirit of RG transformations since topological order is a characterization 
of scale invariant non-local correlations. In fact, there are some analytical and numerical evidences for the existence of 
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QPTs between topological phases and non-topological phases on specific models of parameterized Hamiltonians ll59l - 
l63l . However, the connection between a QPT and geometric shapes of logical operators has not been completely 
established yet. While a change in the number of ground states certainly leads to a QPT as we have seen in the 
discussion on one-dimensional STS models, it is not clear whether the change of geometric shapes of logical operators 
implies the existence of a QPT or not. 

A useful insight in analyzing quantum phases in two-dimensional STS models is to realize the resemblance be- 
tween a mathematical notion of topology and a classification of quantum phases. In classifying geometric shapes of 
objects by using the notion of topology, two objects are considered to be the same when they can be transformed each 
other through continuous deformation, while they are considered to be different when they can be transformed each 
other through only non-analytic changes of geometric shapes. In a similar fashion, when quantum phases are classi- 
fied, two frustration-free Hamiltonians are considered to belong to the same quantum phase when their ground states 
can be connected continuously, while two Hamiltonians are considered to be diff'erent when they can be connected 
only through parameterized Hamiltonians which undergo QPTs. 

The underlying ideas behind topology and a classification of quantum phases are fundamentally akin to each other, 
and the notion of topology has been adopted to classifications of quantum phases in various systems Il35l438l . In STS 
models, geometric shapes of logical operators are central in characterizing the scale invariant physical properties such 
as topological order In particular, it will be shown that the ground state properties drastically changes as a result 



of different geometric shapes of logical operators in Section 5.1 and Section 5.2 Therefore, in classifying quantum 



phases arising in STS models, we introduce the notion of topology into geometric shapes of logical operators. 



In Section 5.3 we show that parameterized Hamiltonians connecting two STS models with different geomet- 
ric shapes of logical operators are always separated by a QPT by proving that the energy gap must close at some 
point. Then, we show that different quantum phases in two-dimensional STS models are characterized by geometric 
shapes of logical operators, and thus, topological characteristics of logical operators can serve as order parameters in 
distinguishing quantum phases arising in two-dimensional STS models. 

5.7. Role of logical operators: concrete examples 

In this subsection, we analyze geometric shapes of logical operators in three specific examples of two-dimensional 
STS models. Then, we analyze topological order arising in these models through geometric shapes of logical opera- 
tors. 



In Section 5.1.1 we begin by analyzing a two-dimensional classical ferromagnet, which is a model without topo- 
logical order, and show that it has an anti-commuting pair of zero-dimensional and two-dimensional logical operators. 
In Section [5.1.2| we discuss the Toric code, an STS model with topological order, and see that it has anti-commuting 
pairs of one-dimensional logical operators. We demonstrate that various characterizations of topological order can be 
explained through geometric shapes of logical operators by reviewing the notion of topological order through anal- 



yses on the Toric code. In Section 5.1.3 we give another example of an STS model with anti-commuting pairs of 
one-dimensional logical operators and compare topological order arising in the model with topological order of the 
Toric code. 

5.7.7. Two-dimensional classical ferromagnet 

We begin by presenting geometric shapes of logical operators in a two-dimensional classical ferromagnet. A 
two-dimensional classical ferromagnet is described by the following Hamiltonian: 

H = _2]z<''-'''Z<''^+"-^Z<'-^''z('+i-''^ (93) 

'J ij 

where Z*'-'^ represents the Pauli operator Z acts on a qubit labeled by a vector (/, j). The total number of qubits is 
N - Li XL2, and the system has periodic boundary conditions. The classical ferromagnet can be viewed as a stabilizer 
code with translation symmetries. The stabilizer group is >Sl,,l2 - ({Z*''-'*Z'^''-'^'\Z*'''*Z'^'^''-'-')v/j). This stabilizer code 
is an STS model with v = 1, satisfying scale symmetries since k = I for any Li and L2 with G(Slj,l^) - L1L2 - 1. 
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Logical operators are 



^ = Z' 



d,i) 



/, /, ■ ■ ■ , / 

/, /, ■ ■ ■ , / 

\Z, I, ■■■ , 1 1 



Y]x^'" 



X, X, ■■■ , X 
X, X, ■ ■ ■ , X 



X, X, ■■■ , X 



(94) 



One can see that both of logical operators satisfy the translation equivalence of logical operators. According to 
geometric shapes of logical operators, we may call I a zero-dimensional logical operator and r a two-dimensional 
logical operator. 

As in a one-dimensional classical ferromagnet, a two-dimensional classical ferromagnet can support a GHZ state 



at zero temperature: |GHZ) 



V2 



(|0 ■ ■ ■ 0) -t- 11 ■ ■ ■ 1)). A GHZ state is stabilized by a two-dimensional logical operator 



r since r|GHZ) - |GHZ). A GHZ state is a globally entangled state since it is a superposition of two globally 
separated ground states |0 ■ ■ ■ 0) and |1 ■ ■ ■ 1). Properties of entanglement may be characterized in a way similar to a 
one-dimensional classical ferromagnet by using a mutual information E(A : E) as we have discussed in Section 4.1 
Note that the model does not have topological order. 

5.7.2. The Toric code as an STS model and topological order 

Next, let us discuss the Toric code El 13, which is a stabilizer code with topological order. We first show that 
the Toric code can be reduced to an STS model by introducing composite particles. Then, we describe geometric 
shapes of logical operators and show that one-dimensional logical operators form anti-commuting pairs. Finally, 
we demonstrate that topological order arising in the Toric code can be characterized by geometric shapes of logical 
operators. In particular, we discuss topological entanglement entropy 1291 . anyonic excitations 111 and the absence of 
local order parameters fT4','42l through logical operators. 

In the Toric code, qubits are defined on edges of a square lattice. The Hamiltonian is (Fig. [T5| ) 



^ = -Z^^-Z- 



.^s 



= n^*^^' ^p = n^'*'- 



(95) 



jes 



jep 



Here, s represent "stars" and p represent "plaquettes" (Fig. 15 i, and j represents the position of qubits. One may see 
that these interaction terms commute with each other, and the Hamiltonian is a stabilizer Hamiltonian. The Toric code 
has A; = 2, as seen from the following equations: 



n^- 



n% 



(96) 



With some observation, one may see that there does not exist any other set of stabilizers whose product becomes an 
identity I. 
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Figure 15: Reduction of the Toric code to an STS model. 
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Reduction to an STS model: The Toric code can be reduced to an STS model by grouping two qubits into a 
composite particle (v - 2) as shown in Fig. 15 By applying some appropriate unitary transformations on qubits 



inside each composite particle, interaction terms yij and Sp can be represented in the following way (Fig. 15 i 
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(97) 
(98) 



A'-fl . 



Here, Zj '^ represents a Pauli operator Zi acting on a composite particle labeled by (/, j). Z\, Z2, X\ and X^ are single 
Pauli operators acting on a single composite particle. 2x2 matrices represent stabilizers graphically. One can see that 
this model satisfies scale symmetries by noticing H/j'^'''''' = I ^nd riy^*'''* = I- 

Now, the two pairs of anti-commuting logical operators in the Toric code can be described in the following way: 
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Here, rii x 112 matrices represent logical operators graphically where ni and 112 are the numbers of composite particles 
in the 1 and 2 directions. According to geometric shapes of these logical operators, we may call them one-dimensional 
logical operators. 

One can easily see that the translation equivalence of logical operators holds from the following equations: 

]^^*i'^) = (iTdn, n^"''^ " '■2T2(r2), []«('•» = ^2^1(^2), n^^''" = '-iTii'-i). 

j '■ j ' 

One can also see that only the pairs of mutually orthogonal logical operators, whose intersections are zero-dimensional, 
may anti-commute with each other. Thus, pairs of anti-commuting logical operators always have zero-dimensional 
intersections. 

Topological order and geometric shapes of logical operators: Having reduced the Toric code to an STS model, 
let us discuss topological order arising in the Toric code through one-dimensional logical operators. Currently, charac- 
terizations of topological order are ambiguous since there are several methods to identify the existence of topological 
order in ground states of correlated spin systems. In order to distinguish quantum phases with different kinds of topo- 
logical order appropriately through geometric shapes of logical operators, it is a necessary first step to demonstrate 
that all the existing characterizations of topological order may be explained through logical operators. Here, we re- 
view these characterizations by analyzing topological entanglement entropy, anyonic excitations and the absence of 
local order parameters in the Toric code through geometric shapes of logical operators. 

(1) Topological entanglement entropy: Ground states of the Toric code have some global entanglement. This 
global entanglement can be quantified by an entanglement measure, called topological entanglement entropy Il46ll47l . 
Here, we follow the definition presented in ll47l . Consider four regions A, B, C and D described in Fig. [16] Then, 
topological entanglement entropy is defined as follows: 



' topo 



— Ee + Er — Ea — Er 



(99) 



at the limit where regions A, B, C and D become infinitely large, and at the thermodynamic limit. Here, Er represents 
the entanglement entropy defined for a region of qubits R. It is known that the topological entanglement entropy Stopo 
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becomes non-zero when a ground state has topological order, while S topo becomes zero when a ground state does not 



have topological order For simplicity of discussion, we set the width of A, B, C and D to be unity (see Fig 16 1 
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Figure 16: Regions used in computing topological entanglement entropy. Each circle represents composite particles. The width of these regions is 
taken to be unity for simplicity of discussion. 



Usually, the computation of topological entanglement entropy is challenging since we need to estimate the quantity 
at the thermodynamic limit. However, this quantity can be easily computed for ground states of stabilizer Hamiltonians 
by considering the numbers of generators for restrictions of the stabilizer groups into A, B, C and D. Let us recall 
that the entanglement entropy can be written as follows: Er -Vr- G{Sr) when there is only a single ground state in 
the stabilizer Hamiltonian. Here, Vr is the total number of qubits inside R and G{Sr) is the number of independent 
generators for the restriction of the stabilizer group, denoted as Sr. This formula also holds for the cases where the 
Hamiltonian has degenerate ground states if there is no logical operator defined inside R. In the Toric code, there is 
no logical operator defined inside A, B, C, or D as we shall show later in the discussion on the absence of local order 
parameters. Thus, we obtain a concise formula for the topological entanglement entropy: 



• topo 



G(Sa) + GiSo) - G(Sb) - GiSc) 
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since Va + Vd -Vb + Vc- 

Now, one can compute the topological entanglement entropy just by counting the number of stabilizers defined 
inside each region. However, most generators appear multiple times and do not contribute to the summation. For 
example, if there is a stabilizer S j defined inside D, S j can be also defined inside A, B and C. Then, such S j does not 
contribute to the topological entanglement entropy. If there is a stabilizer S j is defined inside B (but not inside D), 
S i can be also defined inside A. Such S y does not contribute to the topological entanglement entropy. Then, through 
some observations, one may notice that only the stabilizers defined inside A, but not defined either inside B, C or D, 



contribute to the topological entanglement entropy, 
defined only inside A 



Thus, S topo is equal to the number of stabilizers which can be 
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Such stabilizers can be constructed from stabilizers ^^'•■'^ and S^'-'* in the following way: 






S(x,y) = P]P]S*''^' = 
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Here, jc + 1 x y + 1 matrices represent stabilizers graphically. We have left the entries for identity operators blank 
for brevity of notation. Then, we notice that S topo - 2 for the Toric code as proven in f29l. Though we have 
considered cases where A, B, C and D have the width of unity, it is immediate to generalize the above argument for 
the cases where the width is larger than one. Thus, the Toric code has non-zero topological entanglement entropy, 
which implies the existence of topological order. The connection between geometric shapes of logical operators and 
non-zero topological entanglement entropy will be clarified in the discussion of anyonic excitations, presented soon 
in the below. 

An important observation is that the topological entanglement entropy does not depend on the system size, and 
represents non-local correlations which are scale invariant. Thus, in classifying quantum phases, systems with differ- 
ent topological entanglement entropies must be distinguished appropriately. 

(2) Anyonic excitations: Topologically ordered systems are known to have a finite energy gap between degenerate 
ground states and excited states, which does not vanish even at the thermodynamic limit. Due to this finite energy gap, 
Hamiltonians may support quasiparticle excitations, which are called anyonic excitations. These anyonic excitations 
are localized (involving only a finite number of qubits inside some localized regions), and can be considered as 
particles, called anyons Il64ll65l . When a system has a topological order, these anyons obey unusual quantum statistics. 
Here, we characterize properties of anyons under braiding in terms of geometric shapes of logical operators. 

Anyonic excitations arising in the Toric code can be completely characterized by one-dimensional logical op- 



(1,1) 



7(1,1) 



erators. A pair of anyons can be created by applying Zj ' ' to a ground state of the Toric code li/r). Since Z 
anti-commutes with S("i>"-* and S*"'"'), ZJ ' creates a pair of excitations by flipping eigenvalues of S*"''"-^ and S*"'') 
from 1 to -1 (Fig. [T7|a)). These excitations can be viewed as quasiparticles, in other words, anyons. Anyons can 
propagate around the torus in the 2 direction, and their propagations can be characterized by a logical operator {^ and 



its segment (Fig. 17 a)): 



Zi, /, 

Zi, /, 



£i(yuy2) = Y\z\ 



i,j) 



(103) 



j=yi 



where 1 < yi < ^2 ^ "2- Then, one notices that ^i(3'i,3'2) creates excitations by flipping eigenvalues of S^"'''"" and 
(gini.yi)^ as shown in Fig. 17 a). A pair of anyons can be also annihilated by applying £i{l,n2) = £i. Thus, £1 can 
characterize the creation of anyons and their propagations around the torus in the 2 direction, and their annihilation. 

Here, we call these anyons associated with stabilizers S^''^' "b-type" anyons. Though a logical operators £1 can 
characterize propagations of b-type anyons only in the 2 direction, their propagations in the 1 direction can be charac- 
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terized by a logical operator r2- A pair of b-type anyons can be also created by applying Xj to a ground state of the 
Toric code \ij/) since X^ ' ' anti-commutes with S*^"'"-' and S^-''^' (Fig. 17 a)). Their propagations in the 1 direction 



can be characterized by a logical operator ro and its segment (Fig. 17 a)): 



ri = 



^2, X2, 



I 

X2 



r2{xuX2) = f^^: 



(U) 

2 



(104) 



A similar discussion holds for other one-dimensional logical operators ^2 and ri (Fig. [T7|c)(d)). These logical 
operators characterize anyons associated with stabilizers J\}-''^\ We call these anyons "a-type" anyons. Thus, we 
notice that there are two different kinds of anyons which can travel around the torus because of one-dimensional 
logical operators. 
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Figure 17: (a) Propagations of anyons characterized by a segment of a one-dimensional logical operator f i and r2. Shaded circles represent b-type 
anyons which are associated with stabilizers S'''^'. (b) A segment of .?((a:, y) and propagations of b-type anyons. 

Though we have described propagations of anyons only in the vertical or horizontal direction in the torus, anyons 
can propagate around the torus freely and their propagations are not restricted to the vertical or horizontal direction. 
In order to discuss propagations which do not wind around the torus, stabilizers ^(x,y) and S{x,y) become essential. 
Let us recall that ^(x,y) and S(x,y) have geometric shapes like a loop. We consider the following segment of^(x,y): 



X2, X2, X2, 



X2, X2, 



Zi 
Zi 
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Then, we notice that this segment also creates a pair of b-type anyons at the endpoints of the segment (Fig.[T7|b)). 
With some observations, we also notice that any segment of ^(jc,}') can create a pair of b-type anyons at the endpoints. 
Thus, J[(x,y) can characterize the propagation of b-type anyons around the loop. A similar discussion also holds for 
S(x, y) where S(x, y) characterizes the propagation of a-type anyons around the loop. 

By combining J[(x, y), !B(x, y) and one-dimensional logical operators, one can characterize arbitrary propagations 
of anyons in the Torus. Here, let us analyze the connection between loop-like stabilizers ^{x,y), !B{x,y) and one- 
dimensional logical operators £1, £2, f-i and r2. Here, let us recall the form of ^(x,y) described in Eq. ( 101 1. Then, 



we notice that J{{x,y) includes a segment of Ei inside a vertex in the 2 direction and a segment of r2 inside a vertex in 
the 1 direction. Thus, one can smoothly connect propagations of b-type anyons by attaching J[{x,y) to €1 and r2. A 
similar observation holds for a-type anyons too. 
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Finally, let us show that braiding of a-type and b-type anyons create a non-trivial phase. One can braid an "a-type" 



anyon around a "b-type" anyons by using propagations characterized by ^{x,y) and !B(x,y) (Fig. 18 1. Through the 
braiding, one obtains an extra phase factor "-1" as a result of anti -commutations at the intersection of propagation 
paths characterized by ^(x,y) and !B(x,y). This extra phase factor is a direct consequence of anti-commutations 
between one-dimensional logical operators. Thus, one-dimensional logical operators and loop-like stabilizers ^(x,y) 
and Six, y) characterize the existence of anyonic excitations with non-trivial braiding property. 



b-type 



{ 



A(x,y) 



B{x,yP 
\_/a-type 



Figure 18: Braiding between a-type and 6-type anyons. Their propagations are cliaracterized by ^(x,y) andS(jc, >')■ 

Having analyzed anyons through one-dimensional logical operators, one can see the connection between anyonic 
excitations and topological entanglement entropy clearly. In the computation of topological entanglement entropy, we 
have seen that only the loop-like stabilizers ^{x,y) and S{x,y) contribute to S ,„po- This implies that the topological 
entanglement entropy counts the number of different types of anyons which may propagate around the loop A. This 
observation is consistent with the fact that topological entanglement entropy is related to the total quantum dimension 
of anyons in the following way; 



' topo 



21ogD 



(106) 



where D is the total quantum dimension of anyonic excitations. For the Toric code, we have D - 2 since there are two 
different anyons r] 

(3) Local indistinguishability: An important feature of topologicaUy ordered ground states is that there is no 
local order parameter which can distinguish topological phases. We have already seen that topological entanglement 
entropy is a global quantity which is computed at the limit where regions A, B, C and D become infinitely large and 



at the thermodynamic limit (Fig. 16 1. In general, globally defined quantities or objects are necessary in characterizing 



global entanglement arising in topologicaUy ordered ground states. 

The characterization of topological order we discuss here is closely related to the fact that there is no local order 
parameter to characterize topological phases . In a topologicaUy ordered system, it is known that, for mutually 
orthogonal ground states Itfr^} and |i^i) in a topologicaUy ordered system, the following equations hold: 



<iAo|0|^o) = <(AilO|«Ai>, <^^olO|^i) = 



(107) 



for any locally defined physical observable O at the thermodynamic limit lfT4ll42l . Therefore, one cannot distinguish 
two different topologicaUy ordered ground states through locally defined physical observables. We may call these 
conditions local indistinguishability conditions. 

While this characterization of topological order is rather abstract, it is mathematically convenient in formulating 
topological order and is used commonly in rigorous treatments of topologicaUy ordered systems. Here, we confirm 
that ground states of the Toric code satisfy the above indistinguishability conditions. A useful observation is that there 



^One may say that anyonic excitations are local objects and topological order could be characterized locally. However, topological order refers 
to non-local correlations arising in ground states of topologicaUy ordered systems at zero temperature, and anyonic excitations are signatures of the 
non-local ground state properties. While anyonic excitations are localized objects, their braiding properties still retain scale invariance since the 
non-trivial phase after the braiding does not depend on specific paths where anyons propagated. 
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is no local physical observable O which can change a ground state li^o) to li/ri ). In the language of quantum codes, this 
implies that there is no logical operator I which can be defined locally. Then, we naturally expect that it is sufficient 
to show that there is no locally defined logical operators. 

The Toric code has two pairs of one-dimensional logical operators which are globally defined. Then, one might 
think that there is no local logical operator However, these four logical operators are some particular representations, 
and there exist many equivalent representations of the same logical operators since applications of stabihzers keep 
logical operators equivalent. Thus, one needs to show that one-dimensional logical operators in the Toric code do not 
have any equivalent representations which are defined locally. 

This can be easily proven by a theoretical tool developed for studies of a bi-partite entanglement in stabilizer 
codes (Theorem [T] introduced in Section 2.2 1. In particular, for a given bi-partition of a system of qubits into two 
complementary subsets A and B = A, we have 



gA+KB = 2A: 



(108) 



where §a and gB are the numbers of independent logical operators which can be defined inside A and B respectively. 
Here, we take A as the largest zero-dimensional region which consists of ni - 1 x «2 - 1 composite particles (Fig. 19 1. 
Since the entire system consists of ni xn2 composite particles, B is a union of one-dimensional regions which extend in 
the 1 and 2 directions. Now, since the Toric code has ^ = 2, we have gA+ gB- 4. Then, since all the one-dimensional 
logical operators can be defined inside B, we have gB - 4, and g^ - 0. Thus, there is no logical operator defined 
inside a zero-dimensional region A. 




Figure 19: A bi-partition into a zero-dimensional region A and its complement B. Numbers shown indicate ; 
independent logical operators defined inside A and B respectively. 



and gs which are the numbers of 



Now, consider a physical observable O defined inside some localized region R. Since there is no logical operators 
inside R, the density matrices of all the ground states of the Toric code are the same: pr - -zp^ Yls eSRi^ + S j) where 
S j are independent generators for Sr. Then, one can compute an expectation value of a local physical observable O 
defined inside R in the following way: 



<iA()|0|(Ao> = <<Ai|0|^i> = TrfOpfi]. 



(109) 



One can also easily show that {\j/{)\0\4'\) = by using the fact that there is no logical operator defined inside R. Thus, 
the Toric code satisfy the local indistinguishability conditions. 

Here, we emphasize that one-dimensional logical operators are indeed "one-dimensional" since they cannot be 
defined inside zero-dimensional regions. Therefore, the dimensions we have assigned to logical operators have topo- 
logically invariant meanings which are commonly shared by sets of all the equivalent logical operators |32|. 

A GHZ state, revisited: Here, we make some comments on the difference between a GHZ state and the Toric 
code. While both a GHZ state and a ground state of the Toric code have global entanglement, a GHZ state is not 
topologically ordered. This may be seen by the fact that a GHZ state does not satisfy the conditions (1) and (3). 
While a classical ferromagnet can support localized excitations, these excitations cannot propagate around the lattice 
since the excitation energy increases as excitations move. As for the indistinguishability condition, a zero-dimensional 

0) + 1 1 ■ ■ ■ 1)) and ^(|0 ■ ■ ■ 0) - 1 1 ■ ■ ■ 1». Though a classical 



logical operator Z*^'-'* connects two ground states -i=(|0 
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ferromagnet does not satisfy conditions (1) and (3), a GHZ state has a non-zero topological entanglement entropy since 
S topo - 1 ■ This is because topological entanglement entropy characterizes properties of a single ground state while 
characterizations of topological order need considerations on all the ground states. 

5.7.5. Another model with topological order 

In the Toric code, the existence of one-dimensional logical operators is responsible for topological order. Then, 
we expect that one-dimensional logical operators are central in characterizing topological phases. In order to confirm 
this expectation, let us give another STS model which has one-dimensional logical operators and discuss its physical 
properties. 

We consider a system of Li x L2 qubits governed by the following Hamiltonian: 



where 



H = -^.<' 



J) 



(110) 
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Here, blank entries represent identity operators I. One can see that interaction terms S ^'•■'"' commute with each other. 
In a strict sense, this system is not an STS model since the number of logical qubits change according to Li and L2- 
In particular, with some calculations, we have 
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However, the model can be treated as an STS model when we view 3x2 qubits as a composite particle. Let us define 
«! - L\IT> and n2 - I^il'^ by choosing L^ to be a multiple of 3 and L2 to be a multiple of 2. Then, the system possesses 
scale symmetries since A: = 4 for any «i and n2- 

Now, let us represent stabilizers in this model in terms of composite particles. By applying appropriate unitary 
transformations on qubits inside each composite particle, we can represent stabilizers in the following way: 
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Then, logical operators can be easily found by looking at stabilizers S ^ , ^j ^ S ^ and S^ since they have forms 
similar to y[*''^' and S*'-') in the Toric code. Logical operators are 



^K'-^ni.ni) — 
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(116) 
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Thus, the model has four pairs of anti-commuting one-dimensional logical operators. 

Since the model has logical operators whose geometric shapes are similar to those in the Toric code, we naturally 
expect that the model has similar physical properties as those in the Toric code. Let us confirm this expectation by 
checking the conditions (1), (2) and (3) one by one. First, topological entanglement entropy can be easily computed 
since stabilizers which can be defined only inside a loop-like region A are: 



X y 



si{x,y) = nn^i''' ^2(x,y) = nn^^ 



(118) 



i=\ j=\ 



'=1 i=i 



As a result of the existence of these loop-like stabilizers, this model has S ,opo - 4. Second, anyonic excitations arising 
in this model can be characterized by one-dimensional logical operators. For example, {i and r2 characterize prop- 
agations of anyons associated with stabilizers S^ ir* the 2 and 1 directions respectively. Also, S i{x,y) characterize 
propagations of these anyons around a loop-like region. In total, there are four kinds of different anyons, which may 
be called a-type b-type c-type and d-type anyons, which are characterized by "^i and r2", "^2 and r\\ "^3 and r4" and 
"^4 and r-i". A braiding between a-type and b-type anyons and a braiding betwee c-type and d-type anyons create non- 
trivial phases -1. Other braidings do not generate any phase. Finally, from a similar reasoning used in the discussion 
on the Toric code, it can be shown that ground states of this model satisfy local indistinguishability conditions. 
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Figure 20: Reduction to two copies of tlie Toric codes. Botli models have one-dimensional logical operators with k = 4. 

Now, since topological properties of ground states are similar to those of ground states in the Toric code, we 
naturally expect that this model and the Toric code may belong to the same quantum phase. In fact, the model can 
be reduced to two copies of the Toric code by applying local unitary transformations. From the forms of stabilizers 
5j , ■ ■ ■ ,5g , one may expect that stabilizers S ^ ,■ ■ ■ ,S^ are responsible for the existence of topological order. 
In particular, it can be seen that 5th and 6th qubits inside each composite particle are decoupled from other qubits, 
and are not relevant to topological order. In fact, one can remove S^ and S ^ by applying disentangling operations 
between neighboring composite particles in a way similar to the reduction of a cluster state to a product state used in 
the discussion of one-dimensional STS models. S I , ■ ■ ■ ,S^' look the same as interaction terms yi*''^' and S*''^' in 
the Toric code. Also, it can be seen that 1st and 2nd qubits are decoupled from 3rd and 4th qubits. Thus, one may see 
that two copies of the Toric code are embedded in this model in a non-interacting way. 

Here, in order to see the reduction of the model to the toric code in a more physically intuitive way, we give an 
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actual model which includes two Toric codes: 
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(119) 



with N - Li X L2 qubits where Li and L2 are even integers. Here, qubits at (/, j) with i + j - (mod 2) is decoupled 
from qubits at (/, j) with i + j - 1 (mod 2). If we pick up qubits only at (/, j) with i + j - Q (mod 2) and rotate the 
lattice by 45 degree, we have the following Hamiltonian: 
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Through some observations, we may notice that this model is equivalent to the Toric code. Thus, the Hamiltonian in 



Eq. ( 1 19 1 includes two copies of the Toric code in a non-interacting way. 

In summarizing observations obtained so far, the following two Hamiltonians are equivalent; 
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Since these two Hamiltonians can be transformed each other through local unitary transformations, they belong to the 
same quantum phase. This can be shown through a similar discussion used in the classification of quantum phases in 
one-dimensional STS models. Here, we note that local unitary transformations do not change the geometric shapes of 
logical operators, and thus, logical operators indeed characterize the scale invariant physical properties. 

5.2. Logical operators in two-dimensional STS models and topological order 

We have analyzed possible geometric shapes of logical operators in two-dimensional STS models through several 
examples, and found that logical operators form anti-commuting pairs in the following way: 



Classical Ferromagnet (k - \) 
The Toric code (k - 2) 



0-dim - 2-dim 
1-dim - 1-dim 



We have also seen that one-dimensional logical operators may be responsible for the existence of topological order. 

Here, we wish to establish the relation between physical properties of arbitrary two-dimensional STS models and 
geometric shapes of logical operators. However, it is generally difficult to characterize physical properties in higher- 
dimensional systems (D > 1) both analytically and computationally. In particular, since there are many possibilities 
for geometric shapes of logical operators, one might think that finding logical operators are much more difficult in 
two-dimensional STS models than those in one-dimensional STS models. 

Fortunately, it is possible to determine geometric shapes of logical operators in two-dimensional STS models 
completely. Moreover, we can characterize topological order arising in two-dimensional STS models completely 
from geometric shapes of logical operators. In this subsection, we describe possible geometric shapes of logical oper- 
ators. In particular, we show that possible shapes are "anti-commuting pairs of zero-dimensional and two-dimensional 
logical operators" and "anti-commuting pairs of one-dimensional logical operators". Based on these logical operators, 
we discuss physical properties of arbitrary STS models. While we shall concentrate on presenting geometric shapes 
of logical operators and discussing topological order arising in two-dimensional STS models, all the derivations of 



logical operators are given in Appendix C 



In Section 5.2.1 we present possible geometric shapes of logical operators. In Section 5.2.2 we discuss topolog- 
ical order arising in two-dimensional STS models through geometric shapes of logical operators. 

5.2.7. Geometric shapes of logical operators in two-dimensional STS models 

We present possible geometric shapes of logical operators. In two-dimensional STS models, possible geometric 
shapes are "anti-commuting pairs of zero-dimensional and two-dimensional logical operators" and "anti-commuting 
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pairs of one-dimensional logical operators" (Fig.pTb. Here, we describe the above statement more precisely. Consider 
a two-dimensional STS model with rti x n2 composite particles, k logical qubits and v qubits inside each composite 
particle, which is defined with the stabilizer group >S„|,,„. A region with x xy composite particles is denoted as Ux,y, 
which includes composite particles Pn.ra with 1 < r; < x and I < r2 < y. Then, for wi, n2 > 2v, there exist a canonical 
set of logical operators 



n(^„,,„,) 



A, 



n 



(122) 



which satisfy the following conditions after some appropriate local unitary transformations on qubits inside each 
composite particle: 

• Zero-dimensional and two-dimensional logical operators: For p - 1, ■ ■ ■ , ^o (^o ^ k), (p are logical operators 
defined inside a region U2v,\, and rp are logical operators defined all over the lattice in a periodic way: 



n^?- 



^p^ ^p; 



\_Xp, Xp, 



X, 
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pj 



One-dimensional logical operators: For p -k^ + X,- ■ ■ ,k,lpWQ logical operators defined inside a region with 
ni X 1 composite particles and Vp are logical operators defined inside a region with 1 x «2 composite particles: 



f, = W"' = 



= n 
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Xp, Xp, 



, I 

, Xp 



Computations of this canonical set of logical operators are presented in Appendix C Logical operators are graph 



ically shown in Fig. 21 Here, we have assigned dimensions to logical operators according to their geometric shapes. 
Note that zero-dimensional logical operators are "not completely zero-dimensional" since they are defined inside a 
region with 2v x 1 composite particles which is denoted as U2v,i- For ease of graphical representations, we draw 
logical operators £p (p - 1, ■ ■ ■ , ko) inside a region C/2,.,1 by assuming that they can be actually supported by a single 



composite particle Pj j (or a region Uii) in Fig. 21 



Here, we notice that the sum of dimensions of logical operators which form anti-commuting pairs is always two. 
Also, intersections of anti-commuting logical operators are always zero-dimensional. We call this constraint on the 
dimensions of logical operators the dimensional duality of logical operators in STS models. This duality also holds 
for one-dimensional STS models. We summarize this observation in the following way|j 



Observation 1 (Dimensional duality of logical operators). In D-dimensional STS models (D 
and (D — m)-dimensional logical operators form anti-commuting pairs. 



1,2), m-dimensional 



5.2.2. Topological order in STS models 

From the discussions in the previous subsection, one may naturally expect that an STS model has topological order 
when there exist one-dimensional logical operators. For example, one-dimensional logical operators may characterize 
propagations of anyons in the vertical and horizontal directions. However, whether there exist anyons which can 



We note that there are several models of stabilizer Hamiltonians constructed in higher-dimensional systems which satisfy this dimensional 
duality of logical operators. For example, generalizations of the Toric code defined on a D-dimensional torus (D > 2) have anti-commuting pairs 
of ;?i-dimensional and (D - m)-dimensional logical operators 1 4 66 1 . However, there exists a three-dimensional stabilizer code without scale 
symmetries which may have a pair of anti-commuting two-dimensional logical operators 1671 . 
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Figure 21: A canonical set of logical operators in a two-dimensional STS model. Two-sided arrows represent anti-commutations. Note that 
zero-dimensional logical operators fp (p = 1, ■ ■ ■ , ko) are represented inside a single composite particle instead of a region with 2v X 1 composite 
particles (t/2v,i ) for ease of graphical representations. Dimensions are assigned to logical operators according to their geometric shapes. 



propagate around the torus is not immediately clear. Also, computations of the topological entanglement entropy are 
not possible through only logical operators. 

What are central in characterizing topological order in the Toric code are loop-like stabilizers ^(x, y) and !B(x, y) 
which are generated from stabilizers ^^'■■i'' and S*^''^'. Then, we hope to find such stabilizers and loop-like stabilizers 
generated from them. In fact, there exists such pairs of stabilizers which can characterize propagation of anyons 



around a loop. See Appendix D for derivations and detailed discussions. Here, based on these stabilizers, we discuss 
topological order arising in two-dimensional STS models. 

Let us begin by the case with ki - 2. Then, after some appropriate unitary transformations on qubits inside each 
composite particle, there exist the following pairs of stabilizers: 



-('■,;■) 



P1Z1X2, PiZi 
P1X2, Pi 



(ij) 



:(ij) 



Pi, 
P2Z2, 



P2X1 
P2Z2X1 



id.j) 



(123) 



where Pi and P2 are some Pauli operators acting on a single composite particle. Here, one may notice that S j 

^('•^■> and 5^''-'* = S*'--'' if Pi ^Pi^I. 

Logical operators can be represented in the following way: 



ii.j) 



nzi''^> i2 = n^"' ^' - n^i' ^- - u^ 



(124) 



From 5 J and S^ , one can construct the following loop-like stabilizers: 

Si{x,y) ^ nfl^f^', S2{x,y) = nfl^'-'^l 
i=\ j=i ;=i j=\ 



(125) 



Note that contributions from Pi and P2 cancel out inside the loop. 

Next, let us consider the case with arbitrary ki where ki is the number of pairs of one-dimensional logical operators. 
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For simplicity, let us consider the case where k\ - k and A;o = 0. One can prove that ^i must be even, and there exist 
following stabilizers: 



Ai,j) 



^(ij) 



P1Z1X2, 
P1X2, 

Pi, 
P2Z2, 



PiZi 
Pi 

PiXi 

P2Z2XI 



(ij) 



(ij) 



' k-l 



:(!J) 



Pk-\Zk-\Xk, Pk-\Zk-\ 
Pk-iXk, Pk-i 

Pk, PkXk-i ^*'''' 
PkZk, PkZkXk-\ 



(iJ) 



where Pi, ■ ■ ■ ,Pk are some Pauli operators acting on a single composite particle. One may see that stabilizers form 
pairs just like a pair of ^*''^' and S*''^' in the Toric code. Based on these stabilizers, one can find loop-like stabilizers 
too. 

With these observations, we notice that there are ki different types of anyonic excitations which are characterized 
by pairs of "5 j'' and 82'", ''Sf and S ^ " and so on. These anyons described by each pair of stabilizers are not 
interacting with each other Thus, the braiding rule between anyons is described as a direct product of ki /2 braiding 
rules of the Toric codes. Topological entanglement entropy is ki , which can be directly computed from the number of 
loop-like stabilizers PI 

Finally, let us see whether STS models satisfy local indistinguishability conditions or not. Here, we take the largest 
zero-dimensional region R - t/„,-i,„2-i which consists of ni - I x n2 - I composite particles. Then, we notice that 
gR = ko since g^ - ko + 2k\ where ^0 + ^1 = ^- Thus, when ko = 0, the model satisfies local indistinguishability 
conditions. 

From these discussions, we may conclude that topological properties of two-dimensional STS models can be 
completely characterized from geometric shapes of logical operators. Thus, geometric shapes of logical operators may 
work as order parameters to distinguish topological phases with different topological order arising in STS models. 

5.3. Quantum phases in two-dimensional STS models 

We have seen that different geometric shapes of logical operators lead to completely different ground state prop- 
erties in STS models. Finally, let us classify quantum phases arising in two-dimensional STS models. 

Since a topologically ordered system and a non-topologically ordered system have ground states with completely 
different physical properties, we naturally expect that topological and non- topological STS models belong to different 
quantum phases. However, it is not immediately clear whether the emergence and loss of topological order lead to 
the existence of a QPT or not. In particular, while changes in the number of ground states (the number of logical 
operators) must lead to a QPT, it is still not clear whether there is a QPT or not between two systems with the same 
number of ground states. 

Here, we show that geometric shapes of logical operators can be used as "order parameters" to distinguish quantum 
phases, including topological phases, arising in two-dimensional STS models by proving that any parameterized 
Hamiltonians connecting two STS models with different geometric shapes of logical operators are always separated 
by QPTs. We also show that the existence of a QPT originates from the non-analyticity of transforming logical 
operators with topologically distinct shapes each other since there is no continuous deformation (diffeomorphism) 
between them at the thermodynamic limit. 

In Section 5.3.1 we show that two STS models with topologically different logical operators are always sepa- 



rated by a QPT. In Section [5. 3. 2| we expand our discussions on topological structures of logical operators to sets of 
equivalent logical operators. 

5.5.7. The vanishing energy gap and topological QPT 

Let us begin by clarifying the problem we address. Consider two STS models Ha and Hb with and without 
topological order When the number of logical operators in Ha and Hb are different, there must always be a QPT 
between Ha and Hb- Therefore, we consider the cases where Ha and Hb have the same number of logical operators, 
but different geometric shapes of logical operators. In particular, for simplicity of discussion, we consider the case 



'in order to prove this rigorously, one needs to show that there is no other stabilizer which can be defined inside a loop-like region. This proof 
is not trivial, but can be obtained by using discussions similar to proofs of lemmafTlwhich appear inlAppendix C.ll 

46 



where Ha has two pairs of anti-commuting zero-dimensional and two-dimensional logical operators {k^ - 2, ki =0 
and k - 2) while Hb has two pairs of anti-commuting one-dimensional logical operators (kd - Q,k\ -2 and k - 2). 
One can generalize this discussion to the cases with arbitrary numbers of logical operators easily. 

Our goal is to show that H^ and Hg belong to different quantum phases. For this purpose, we suppose the opposite; 
there exists a parameterized Hamiltonian Hie) which connects Ha and Hb without closing the energy gap or changing 
the number of ground states: 

HiQ) = Ha, H(l) = Hb, A(e) > c for all e (126) 

where A(e) is the energy gap and c is some constant which does not depend on the system size. 

The most striking difference between Ha and Hb is the presence and the absence of topological order which can 
be seen by different geometric shapes of their logical operators. This drastic change of geometric shapes may underUe 
a non-analytic change of the ground state properties and a vanishing energy gap. Then, one might hope to show 
the existence of a QPT by looking at how geometric shapes of logical operators change with a parameter change. 
However, unlike Ha and Hb, one cannot define logical operators inside the Pauli group for //(e) at < e < 1 since the 
Hamiltonian H(e) is frustrated in general. 

Despite //(e) is frustrated and not a stabilizer Hamiltonian in general, it is possible to define operators which 
act like logical operators. Let us denote projection operators onto the ground state space of //(e) as P(e). Since the 
number of degenerate ground states does not change with e, one can find some unitary transformation which satisfy 
the following condition: 

P(e) = U(e)pAU(er\ Pa = P(0), Pb = P(l) (127) 

where Pa and Pb represent projections onto the ground state spaces of Ha and Hb respectively. Here, we note that such 
a unitary transformation is not uniquely determined and has many degrees of freedom. The unitary transformation 
U(e) may characterize the evolution of ground states with respect to e. Let us pick up some ground state |i^(0)) of Ha 
at e = 0. Then, the following state |i/'(e)) - U{e)\if){Q)) is a ground state of //(e). Here, we denote an anti-commuting 
pair of zero-dimensional and two-dimensional logical operators in Ha as {a and rA- Then, the following operators act 
like logical operators inside the ground state space of //(e), transforming degenerate ground states among them: 

m = U{e)€AU(e)-\ r(e) = U{e)rAU{e)-' (128) 

since they act non-trivially inside the ground state space: 

K(e),r(e)) = 0, {{ef = r(ef = /. (129) 

Note that ^(0) = ( and r(0) = r. We call €{e) and r(e) analytically continued logical operators. Note that tie) and r(e) 
may not commute with the parameterized Hamiltonian //(e). 

Now, let us discuss how analytically continued logical operators l{e) change with respect to e. Here, we first give 
an intuitive explanation why the assumption of no QPT leads to a contradiction. Then, we add some mathematical 
rigor to our intuition. 

Intuitive approach based on topology in logical operators: At the beginning of the discussion, we have assumed 
that there is no QPT between e = and e = 1 . Then, we naturally hope that the geometric shapes of {{e) change 
smoothly with respect to e. A unitary transformation U{e) must not have any non- analytic ity with respect to e even 
at the thermodynamic limit since a ground state |iA(e)) must change continuously with e. Then, £{e) also changes 
continuously with respect to e without non-analyticity. At e = 0, a zero-dimensional logical operator ^(0) = (a is 



defined inside some localized region (Fig. 22 1. For e sufficiently close to zero, ^(0) and {{e) are similar, and £{e) may 
be approximated as: 

({e) ^ {{Q) + e{{Q)' (130) 

where ^(0)' represents the derivative of {(e) with respect to e at e = 0. Then, the geometric shape of l{e) should be 
similar to the one of /'(O) (Fig.[22]l. 

47 



Now, let us analyze geometric shapes of ^(e) at e = 1. At e = 1, the analytically continued logical operators (!(l) 
and r(l) act non-trivially inside the ground state space of Hg- Let us recall that there are no logical operators inside 
any zero-dimensional region R in an STS model Hg'- gu - where R is a region consisting of (ni - 1) x (n2 - 1) 
composite particles. As a result, ground states of Hg satisfy the local indistinguishability conditions, and there is no 
local operator which can transform ground states each other. As a result, neither ^(1) or r(l) can be defined locally. 



Energy 





logical 
operators 



® 



Topologically distinct 
no continuous defornnationl 



non-analyticity 

i 






Figure 22: Geometric sliapes of an analytically continued logical operator, the vanishing energy gap, and the existence of a QPT. Lightly shaded 
regions represent the "shapes" of analytically continued logical operators. 

While an analytically continued logical operator £{l) acts like a logical operator of Hb, it may not be an actual 
logical operators of Hg since it is not a Pauli operator in general. Here, for simplicity of discussion, we assume 
that ^(1) happens to be a one-dimensional logical operator fg in Hb which is a Pauli operator. Then, geometric 
shapes of {(e) changes from a zero-dimensional object {^ to a one-dimensional object ^^ as we vary e from e = to 
e - 1. However, it is impossible to change geometric shapes of a zero-dimensional object to a one-dimensional object 
continuously at the thermodynamic limit since a continuous deformation between them is topologically prohibited. 
This contradicts with our original assumption that there is no QPT in H(e). Thus, there must be a QPT between Ha 
and Hb, and we may conclude that Ha and Hb belong to different quantum phases. Though we have assumed that 
^(1) happens to be a one-dimensional logical operator ^b in Hb, the observation that ^(1) cannot be defined locally is 
sufficient to reach the same conclusion. 

A more rigorous approach based on adiabatic continuation: While the above analysis implies that the topo- 
logical distinction between geometric shapes of logical operators in Ha and Hb lead to the existence of a QPT, the 
discussion may lack in a mathematical rigor. In particular, we have assumed that geometric shapes of ^(0) and {(e) for 
small 6 are close. However, it is not clear if this assumption is correct or not. Also, while we have used the expression 
"geometric shape" naively to describe geometric properties of analytically continued logical operators {(e), we have 
not stated clearly the definition of geometric shapes of analytically continued logical operators. 

Below, we shall clear these ambiguities by borrowing theoretical techniques developed in studies of adiabatic con- 
tinuation lfT3l [141 l42l l68ll ■ Adiabatic continuation is an idea of studying the ground state properties of Hamiltonians 
belonging to the same quantum phase by finding a gapped parameterized Hamiltonian connecting them. Here, for 
simplicity of discussion, we begin with the cases where there is always a single ground state in parameterized Hamil- 
tonians. If two Hamiltonians H and H' are in the same quantum phase, one can always find gapped parameterized 
Hamiltonian H(e) connecting them due to the definition of quantum phases. Then, according to the adiabatic theo- 
rem [69 70], by varying the parameterized Hamiltonian H(e) sufficiently slowly, one can transfer a ground state of H 
(denoted as \tp}) to a ground state of H' (denoted as liA')). A key idea is to realize that a unitary operator U required to 
transform \t(r} to \ifr'} di/r') - U\iff}) can be extracted from the parameterized Hamiltonians H(e). In particular, theoret- 
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ical techniques on adiabatic continuation provide systematic formulas to obtain such a unitary transformation U from 
a gapped parameterized Hamihonian //(e). While we have started our discussions with the cases where there is only 
a single ground state, adiabatic continuation also works when there is the ground state degeneracy and the number of 
degenerate ground states does not change. Even more, it works when the ground state degeneracy is slightly broken. 
(See references for applicabilities of adiabatic continuation lT3 HT4ll42ll68]| ). 



Now, we interpret our previous analyses on the parameterized Hamiltonians in Eq. ( 126 1 through adiabatic contin- 
uation. Let us recall that we have supposed that there exists a gapped parameterized Hamihonian Hie) which connect 
Ha and Hb, and the number of ground stated does not change as we vary e. Then, one can represent the unitary 



transformation [/(e) defined in Eq. ( 127 i through the parameterized Hamiltonian //(e). Now, we discuss properties 
of analytically continued logical operator tie) - U{e)iAU{ey'^ . A remarkable discovery from recent developments 
on studies of adiabatic continuation is the fact that one can approximate analytically continued logical operators {{e) 
with some operator {(e) whose geometric shape is close to the one of (a as long as a parameterized Hamiltonian //(e) 
remains gapped at the thermodynamic Umit lfT3l [141 . In particular, there exists an approximation €{e) defined inside 
some localized region with ^(e) x ^(e) composite particles: 

tie) ~ m (131) 

where ^(e) depends on the energy gap A(e), but remains finite even at the thermodynamic limit. Here, by approxima- 
tion, we mean that lie) and lie) act in a similar way inside the ground state space. In particular, if lie) transforms two 
degenerate ground states |i/'(e)) and |i/''(e)) of Hie) in the following way: 

lie)me)} = me)'}, <f (e)K(e)|«A(e)) ~ 0(1). (132) 

Here, the inner product is some constant which does not depend on the system size. 

Since lie) approximates lie), one may consider a localized region with ^(e) x ^(e) composite particles as the 
geometric shape of lie). Now, at e = 1, we have an approximation ?(1) which is defined inside some finite region with 
^(1) X ^(1) composite particles. Then, we have 

{<P'il)\lil)mi)} ~ Oil) (133) 

for two orthogonal ground states |i/'(l)) and |i^'(l)) of //(I) = Hb- However, from the indistinguishability condition, 
there is no local physical observable which can transform ground states each other in Hb- Thus, ^(1) must be an 
infinite number at the thermodynamic limit, which leads to a contradiction. Therefore, there must be a QPT between 
Ha and Hb since the energy gap vanishes or the number of degenerate ground states changes at some point. 

QPTs and topology in logical operators: We have seen that non-analytic changes of geometric shapes of logical 
operators lead to the existence of the vanishing energy gap. This observation may be better understood by the use of a 
mathematical language developed in studies of shapes of objects. Topology is a study of classifying objects in terms 
of smoothness and non-analyticity. If two objects can be transformed each other through continuous deformations 
(dififeomorphism), they are considered to be the same. Roughly speaking, diffeomorphism is a one-to-one mapping 
between two geometric manifolds where both the map itself and its inverse are differentiable. The notion of topology 
has been particularly useful in characterizing physical properties which may survive even at the thermodynamic limit, 
mainly in the context of field theory where physical properties are to be discussed at the continuum limit. 

In a strict sense, the notion of topology cannot be introduced in discussions of quantum phases arising in lattice 
systems since spins on lattices are discretely distributed and geometric shapes of logical operators are not smoothly de- 
termined. However, by considering a system of qubits at the thermodynamic limit, one can smoothen geometric shapes 
of logical operators effectively. Then, the fact that there is no continuous unitary transformation C/(e) transforming 
^(0) into ^(1) is a direct consequence of the fact that there is no diffeomorphism which map a topologically trivial 
object (a zero-dimensional point) to a topologically non-trivial object (a one-dimensional loop winding around the 
torus). One cannot cut the winding of a one-dimensional logical operator without introducing some non-analyticity. 

Summary and application: We have shown that Ha and Hb belong to different quantum phases when geometric 
shapes of their logical operators are different. Though discussions have been limited to the cases where Ha has two 
pairs of anti-commuting logical operators and Hb has two pairs of anti-commuting logical operators, our analysis can 
be readily generalized to arbitrary two-dimensional STS models. 
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Let us consider the case when Ha and Hb have various numbers of pairs of anti-commuting logical operators. In 
particular, let us denote the numbers of logical qubits as fe^'** and k^^\ and the numbers of zero-dimensional logical 
operators kf^ and fejf * with kf^ = fe*^' - kf^ and kf = k^"^ - fe<^': 



( 1 dim)-( 1 dim) (0 dim)-(2 dim) 
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Hb 
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kf> 
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Then, two STS models belong to different quantum phases when k\ 

Here, we make some comments on the cases where two STS models have the same geometric shapes of logical 
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operators: k^ - fe^ and k\ - k\ . Since topological properties of STS models, including topological entanglement 
entropy and anyonic excitations, can be completely characterized by one-dimensional logical operators, it is natural 
to expect that Ha and Hb belong to the same quantum phases. In fact, in Section 5.1.3 we introduced the model 



which has four pairs of anti-commuting logical operators and showed that the model is equivalent to a model where 
two copies of the Toric code are embedded in a non-interacting way. In this case, two STS models can be connected 
by some parameterized Hamiltonian without closing an energy gap since two models can be transformed each other 
only through local unitary transformations acting on neighboring composite particles. 

At this moment, we do not have a mathematical proof to the statement that Ha and Hb belong to the same 
quantum phases when geometric shapes of their logical operators are the same. However, we believe this statement 
is reasonable since any existing characterization of topological order cannot distinguish two STS models with logical 
operators whose geometric shapes are the same. Therefore, we conjecture that two STS models belong to the same 
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quantum phases if and only if geometric shapes of their logical operators are the same: k\ — k^ 
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Figure 23: A "phase diagram" of two-dimensional STS models. 



Here, we summarize the main result of this section (Fig. 23 i 



Different quantum phases in two-dimensional STS models can be characterized by not only the number of logi- 
cal operators, but also geometric shapes of logical operators. There must be a QPT when there is a topologically 
prohibited change in geometric shapes of logical operators between two STS models. Geometric shapes of log- 
ical operators can be used as order parameters to distinguish quantum phases arising in two-dimensional STS 
models. 



Note that a QPT with a change of the numbers of one-dimensional logical operators involves the emergence or 



the loss of topological order Such QPTs are called topological QPT (TQPT). In Fig. 23 a QPT which crosses the 
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boundary in the lateral direction is a topological QPT, while a QPT which crosses the boundary in the vertical direction 
is a non-topological QPT. 

Now, let us look at some examples. As a result of our analysis on geometric shapes of logical operators, one can 
classify quantum phases arising in STS models in the following way: 
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Here, Ha is the Toric code, Hb is a Hamiltonian supporting a single product state, H'^ is a Hamiltonian supporting a 
two-dimensional cluster state, H^ is a classical ferromagnet and Hd and //^ are Hamiltonians discussed in Section 5. 1 



Ha and H'^ have two pairs of anti-commuting one-dimensional logical operators, Hg and H'^ have no logical operators. 
He has a pair of zero-dimensional and two-dimensional logical operators and Ho and //^ have four pairs of anti- 
commuting one-dimensional logical operators. 



5.3.2. Topological deformation of logical operators 

While we have seen that topological structures of logical operators play crucial roles in the classification of quan- 
tum phases, our discussions have been hmited to some specific representations of logical operators. However, logical 
operators have many equivalent representations and their geometric shapes may not be uniquely determined. We 
conclude the discussion of this paper by expanding our analyses on topological structures of logical operators to a 
set of equivalent logical operators. In particular, we show that the notion of continuous deformation naturally arises 
in geometric shapes of a set of equivalent logical operators in two-dimensional STS models. We note that a similar 
analysis on geometric shapes of a set of equivalent logical operators in the Toric code was presented in |32|. 

We begin by analyzing geometric shapes of two-dimensional logical operators r^ (1 < p < ko). Let us recall that 
rp anti-commutes with a zero-dimensional logical operator fp. Here, the translations of £p are equivalent to £p due 
to the translation equivalence of logical operators (Theorem l2]|. Then, we notice that all the logical operators which 
are equivalent to rp must have supports on every composite particle since rp needs to have some overlaps with all the 
translations of fp in order for rp to anti-commute with l^p. With this observation, we conclude that two-dimensional 
logical operator rp can be defined only inside a two-dimensional region, with supports on all the composite particles. 

Next, let us analyze one-dimensional logical operators £p and rp (ko + I < p < k). In particular, we show 
that one-dimensional logical operators cannot be defined inside any zero-dimensional region Ua,b with a < «i and 
b < n2. In order to show this, it is sufficient to prove that one-dimensional logical operators cannot be defined inside 
a region with n^ - I x n2 - I composite particles. Suppose that £p can be defined inside a region with n^ - I x n2 - I 
composite particles. Then, since rp can be defined inside a region with 1 x n2 composite particles, there exists a 
translation of rp which does not have an overlap with fp. Then, due to the translation equivalence of logical operators. 



Cp and rp commute with each other However, this contradicts with the fact that 



'p and rp anti-commute with each 



other Through similar discussions, we can show that any one-dimensional logical operators cannot be defined inside 
zero-dimensional regions. 

So far, we have seen that dimensions of logical operators have topologically invariant meanings where w-dimensional 
logical operators cannot be defined inside (m - l)-dimensional regions (m - 1,2). Now, let us discuss the changes 
of geometric shapes of logical operators. For this purpose, let us fix some notations. We call a region with ni x 1 
composite particles 2(1) and a region with 1 x n2 composite particles 2(2)- 

A useful observation regarding the topological structures of logical operators can be obtained by considering the 
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number of independent logical operators gA which can be defined inside a region A. Let us consider the case where 
we have two regions A and A' where A is larger than A', meaning that A includes all the composite particles inside A' . 
Then, if ^^ = gA', one can deform all the logical operators defined inside A into A', since logical operators defined 
inside A must have equivalent representations supported inside A'. With this observation in mind, the following 
theorem is central in characterizing topological structures of a set of equivalent logical operators. 



Theorem 3 (Topological deformation of logical operators). 
two-dimensional STSs, we have the following equations: 



For the numbers of independent logical operators in 
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Figure 24: (a) A deformation of logical operators from Pi.i to 2(1) U 2(2). (b) Q(l) U 2(2) and Pij. (c) 2(1) and 2(1)- (d) 2(2) and 2(2). 



Here, let us interpret the meaning of the theorem. Let us begin by analyzing gp^ ^ - gQ(i)uQ(2) (Fig. 24 a)). Pii and 
Q(1)UQ{2) are the regions with windings in the 1 and 2 directions, and one can deform a logical operator as long as we 
do not break the windings in geometric shapes of a logical operator. Next, let us analyze gpnjjQnS - gP\ i (Fig- 24 b)). 



Both Q(\) U 2(2) and Pj i are one-dimensional regions without any winding, and one can deform a logical operator 
defined inside 2(1) U 2(2) until it becomes a point Pi j. Finally, let us analyze gQ{i) - gQ(\) (Fig. 24 c)). Both 2(1) 



and 2( 1 ) have a winding in the 1 direction, but do not have a winding in the 2 direction, and one can deform a logical 
operator as long as we do not break the winding in the 1 direction. A similar discussion holds for gQQ.) - gQ{2) 
(Fig.|24jd)). The discussions so far can be summarized in the following way. 

• One can deform the geometric shapes of logical operators continuously while they remain equivalent in the 
following ways: 



^1, 

where, in "P — > P" 
defined inside /J'. 



2(1) u 2(2) 2(1) u 2(2) ^ Pi, 



2(1) ^ G(i) 2(2) ^ 2(2) 



if a logical operator i defined inside R is given, then there exists a logical operator €' 



Proof of theorem l3 : One can easily prove four equations through a simple counting of the number of logical 
operators with Theorem T] Let us start with the proof of gp^ ^ - gQ(i)uQ(2)- Since two-dimensional logical operators 
cannot be defined inside Pij, we have gp^ ^ - 2ko + ki. Also, since all the one-dimensional and zero-dimensional 
logical operators are defined inside 2(1) U 2(2), we have ^g(i)uQ(2) = 2feo + ^i, and thus gp^ ^ - gg(i)ue(2)- Now, since 
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gA+8B = 2k for B = A for any region A, gp, , = gQ(i)uQ(2) leads to g g(i)ui2(2) = ^^u ■ ^l^o, by setting A = Q(l), we 
have gQ(i) = gQ(i) - k. This completes the proof. D 

As a result of theoreml3] we notice that there exist topologically invariant properties which are commonly shared 
by a set of equivalent logical operators. The theorem essentially states that one can deform a geometric shape of a 
given logical operators freely while keeping it equivalent as long as we do not change its topological structure. We call 
this property of logical operators in two-dimensional STS models topological deformation of logical operatorsF] 
Therefore, the notion of topology arises naturally in geometric shapes of logical operators in STS models through 
topological deformation of logical operators. 

6. Summary, open questions and outlook 

In this paper, we have discussed possible quantum phases and their classifications in two-dimensional spin systems 
on a lattice. Our results may be summarized in the following three points. 

• STS model: We have proposed a model of frustration-free Hamiltonians which covers a large class of phys- 
ically realistic stabilizer Hamiltonians which are constrained to only three physical conditions; the locality of 
interactions, translation symmetries and scale symmetries. 

• Exact solution of the model: As a key to the analyses on physical properties of STS models, we have found 
possible forms of logical operators and their geometric shapes completely. Then, we have characterized topo- 
logical order arising in the model by geometric shapes of logical operators. 

• Quantum phase transitions and logical operators: We have shown that different quantum phases in STS 
models can be characterized by geometric shapes of logical operators. We have shown that the existence of a 
QPT results from non-analytic changes of geometric shapes of logical operators. 

Since possible geometric shapes of logical operators in two-dimensional STS models are pairs of zero-dimensional 
and two-dimensional logical operators as in a classical ferromagnet, or pairs of one-dimensional logical operators as 
in the Toric code, we conjecture that the model can be reduced to a system where classical ferromagnets and the Toric 
code are embedded in a non-interacting way. However, at this moment, we do not have a proof for it, and would like 
to leave it as an open question for the future. 

Our construction of STS models and classification of quantum phases based on geometric shapes of logical oper- 
ators can likely be broadened in many ways. 

Qudit stabilizer codes and non-abelian quantum phases: While we have studied only the systems of qubits 
(spin 1/2 particles), systems of qudits (particles with larger spins) can provide rich varieties of quantum phases. There 
exists an extension of stabilizer codes to systems of qudits by the use of some generalization of Pauli operators iTTl. 
It is possible to construct STS models based on this qudit stabilizer formalism. Then, our analyses on STS models 
may be readily generalized to STS models constructed on systems of qudits. 

Recently, topologically ordered systems with anyonic excitations whose braiding rule is characterized by a non- 
Abelian group have been gathering significant attentions since such systems with non-Abelian topological order may 
be used as a resource for realizing fault-tolerant quantum computation 13|. Since the stabilizer formalism is based 
on a set of Pauli operators which commute with each other, anyonic excitations supported by stabilizer Hamiltonians 
obey a braiding rule characterized by an Abelian group, as in the Toric code. However, the Toric code has natural 
extensions, constructed based on not Abelian groups, but non-Abelian groups |5 1. These extensions of the Toric code, 
called the quantum double model, support non-Abelian anyons. In the quantum double model, there are operators 
which are analogous to logical operators in stabilizer codes. Thus, our analyses and classifications of quantum phases 
through geometric shapes of logical operators may also be applied to non-abelian topological phases too. 

'^While our discussions have been limited to two-dimensional stabilizer Hamiltonians, logical operators in several stabilizer Hamiltonians 
constructed in higher-dimensional systems also have the similar property. For example, generalizations of the Toric code defined on D-dimensional 
torus (D > 2) have m-dimensional and D - m-dimensional logical operators |4 66 1, and they obey the topological deformation of logical operators. 
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Also, while we have limited our considerations to a QPT between a non-topologically ordered system and a 
topologically ordered system, one can consider a QPT between different topological phases. It may be also possible 
to capture such a QPT in terms of geometric shapes of logical operators. 

Higher-dimensional STS models: We have analyzed two-dimensional STS models as a class of stabilizer codes 
which build on physically reasonable systems. However, correlated spin systems in higher-dimensional systems also 
gather significant attentions in quantum information science community, concerning two important open questions in 
quantum coding theory. 

One of the ultimate goals in quantum coding theory is to create a self-correcting quantum memory |31 1. A self- 
correcting quantum memory is an idealistic memory which would correct errors by itself due to the large energy barrier 
separating degenerate ground states. If such a memory could exist, it will be a perfect quantum information storage 
device which may be used commercially in the future. Though there have been significant progresses in hypothetical 
constructions of a self-correcting quantum memory in a four-dimensional space Il4l l72ll . convincing proposals for its 
realization in three-dimensions have not appeared yet. It is shown that a self-correcting stabiUzer codes cannot exist 
in two dimensions f3T|. Now, the feasibility of a self-correcting memory in three-dimensions is one of the most 
important and interesting open questions in quantum coding theory Il30l 17211741 . 

The feasibility of a self-correcting quantum memory is closely related to another important open question con- 
cerning the upper bound on the code distance of local stabilizer codes. It is commonly believed that the code distance 
of local stabilizer codes with A^ qubits is upper bounded by 0{ ^fN) polynomially in the A^ — > oo limit, rl However, 
recently it have been shown that the code distance of local stabilizer codes is upper bounded by 0(L ) in D di- 
mensions where L is a linear length of the stabilizer code f3T|. Though this work has opened possibilities for the 
existence of a local stabilizer code whose code distance may exceed 0( y/N) polynomially, this bound is proven to be 
tight only for D = 1,2. The upper bound on the code distance of local stabilizer codes is also one of the important 
open questions in quantum coding theory. 

Unfortunately, our analyses in this paper cannot give answers to these open questions since our discussions are 
limited to two dimensions. However, STS models may serve as a class of quantum codes which are physically 
realizable, and their analyses may provide partial answers toward these open questions f76l. 

Scale symmetries and weak breaking of translation symmetries: The final problem we address is rather con- 
ceptual, but may have some fundamental importance. While the locality of interactions and translation symmetries 
are important physical constraints, it may not be obvious why scale symmetries are also important. Here, we mention 
the importance of scale symmetries in a relation with translation symmetries of ground states. 

Topologically ordered systems are known to have degenerate ground states, and in analyzing topologically ordered 
systems, we wish to study properties of not only a single ground state, but the entire ground state space. In analyzing 
degenerate ground states of topologically ordered systems, a main challenge is the fact that translation symmetries 
may be broken. For example, consider a translation symmetric Hamiltonian which is invariant under unit translations 
of qubits. One might hope that this Hamiltonian would give rise to ground states which are also invariant under unit 
translations of qubits. However, this naive expectation is generally true only when there is a single ground state. 
There exist examples where degenerate ground states are invariant only under translations of several qubits, while the 
Hamiltonian is invariant under unit translations of qubits. This bizarre breaking of translation symmetries is observed 
commonly in topologically ordered systems, and is sometimes called a weak breaking of translation symmetries |[9l. 
When translation symmetries are weakly broken, there will exist a ground state which diff'ers from its own translation. 
Then, naturally, we wish to further coarse-grain the system so that all the ground states are invariant under unit 
translations of coarse-grained particles. 

Now, a naturally arising question is whether translation symmetries of STS models are broken or not. We have 
coarse-grained the system of qubits by introducing composite particles so that Hamiltonians are invariant under unit 
translations of composite particles. However, it is not clear whether ground states of STS models are also invariant 
under unit translations of composite particles or not. If translation symmetries of ground states are broken weakly, 
STS models must be further coarse-grained so that ground states restore translation symmetries with larger composite 
particles. 



'We note that there exists a local stabilizer code whose code distance scales as 0( yN log N) |75| . Also, there is an example of three-dimensional 
stabilizer Hamiltonians without scale symmetries whose code distance may scale as 0(L?) for some specific choices of the system sizes 1671 . 
However, no example of stabilizer codes has been found whose code distance surpass 0( ^/N) "polynomially" regardless of the system size. 
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Somewhat surprisingly, in STS models, all the degenerate ground states are always invariant under unit translation 
of composite particles, and translation symmetries are not broken weakly while STS models may have topological 
order. We shall give a proof for this statement in [Appendix E| A key for the existence of translation symmetries 



of ground states is the existence of scale symmetries in stabilizer Hamiltonians. In fact, it is shown in Appendix E 
that for stabilizer Hamiltonians with translation symmetries, the existence of scale symmetries is the necessary and 
sufficient condition for translation symmetries of ground states. Then, scale symmetries may be the sufficient and 
necessary condition for the presence of translation symmetries in ground states in arbitrary gapped correlated spin 
systems. 

Currently, the importance of scale symmetries is not well appreciated, while the locality of interaction terms and 
translation symmetries of ground states are well respected as important constraints for physical realizability of system 
Hamiltonians. However, it will be interesting to see the roles of scale symmetries in various quantum many-body 
systems. In particular, the use of scale symmetries may enable us to have some deeper insights on the underlying 
mechanisms behind correlated spin systems, as we have succeeded in introducing the notion of topology into the 
discussions of logical operators thanks to scale symmetries. Possible applications of scale symmetries may include 
studies of computational complexities of Hamiltonian problems [77 J , RG algorithms based on the tensor product 
state representations |fT9l l20l and analyses on coding properties of quantum codes beyond stabilizer codes such as 
subsystem codes Il30ll78ll79l. 

Appendix A. Proof of the translation equivalence of logical operators (Theorem |2j 

In this appendix, we give a proof of the translation equivalence of logical operators (Theorem |2|. In |Appendix| 
we give a proof by using a certain lemma concerning properties of logical operators in STS models (lemma fTb. 
In Appendix A. 2 we give a proof of this lemma to complete the proof of Theoreml2] In Appendix A.3| an extension 



A.l 



of this lemma is presented for the sake of later discussions in finding logical operators in two-dimensional STS models. 

Appendix A.l. Proof of the translation equivalence of logical operators 

We give a proof of the translation equivalence of logical operators. 

Translation equivalence of logical operators (Theorem |2|: For each and every logical operator f in an STS 
model, a unit translation of { with respect to composite particles in any direction is always equivalent to the original 
logical operator {: 

Tje) ~ e, Mi e L,t im = I,--- ,D) (A.l) 

where L,^ is a set of all the logical operators for an STS model defined with the stabilizer group Sfi. 

Proof of Theorem^ It is sufficient to prove the theorem for translations in the 1 direction: Ti{t) ~ (. Here, we 
define the following {D - l)-dimensional regions (hyperplanes) (see Fig. |A.25] i: 

R(x) - Pr|_..._r„, : ri - X and 1 < r,„ < «,„ for m ^ 1 I- (A. 2) 

In total, there are «i of(D- l)-dimensional regions /?(x) for jc = I,--- ,ni. Then, the following lemma holds. 
Lemma 1. In an STS model defined with the stabilizer group Sft, there exists a canonical set of logical operators 

n(^,) = {^;' :::; i] (a.3) 

such that tp are defined inside a (D — l)-dimensional region R(l) (p — 1, ■ ■ ■ ,k) while rp are defined over the entire 
lattice in a periodic way in the 1 direction: 

Tiirp) ^ rp (p ^ !,■■■ ,k). (A.4) 
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The meaning of this lemma becomes clear when all the logical operators are graphically shown as in Fig. A. 25 
Here, rp are represented with Pauli operators r' acting on /?(!) in the following way: 



- WTlir'p). 



(A.5) 



With this canonical set of logical operators, one can prove that all the logical operators remain equivalent under 
unit translations in the 1 direction. In fact, one can immediately see that logical operators r^ satisfy the translation 
equivalence since they are periodic: T\{rp) - rp. One can also prove that ip ~ T\{{p) by seeing that lpT\{(p) 
commutes with all the logical operators listed in a canonical set n(>S,t) given in lemma fl] and thus, ipT\({p) must be 
a stabilizer. By considering the fact that any logical operator can be represented as a product of ^i , ■ ■ ■ ,(k, r\,- ■ ■ , r^ 
and stabilizers in Sj{, one obtains 



Ti{{) ~ (, V£ e U 



This completes the proof of Theoreml2] 






(A.6) 
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Figure A. 25: A canonical set of logical operators. A two-dimensional example is illustrated (D ■■ 
dimensional regions R(x). (Here, R{x) are one-dimensional regions for D = 2). 



2). The entire lattice is separated into (D - 1)- 



Appendix A.2. Proof of lemmaU] 

Now, we give a proof of lemma [T] We begin our proof by proving the following lemma. 

Lemma 2. Consider a stabilizer code defined with the stabilizer group S, the centralizer group C and k logical qubits. 
If there exist centralizer operators £p,rp e C (p — I,- ■ ■ ,a) with a < k which satisfy the following commutation 
relations 






(A.7) 



{p and rp are independent logical operators. 



In other words, if we find pairs of anti-commuting centralizer operators represented in a canonical form, they are 
guaranteed to be independent logical operators in the stabilizer code. The proof of lemma l2] is immediate by seeing 
that any product of operators taken among (p and r^ (p = 1 , ■ ■ ■ , a) is not a stabilizer For example, if {\ is included 
in the product, ri anti-commutes with the product, and the product is not a stabilizer. Thus, these 2a logical operators 
tp and rp are independent. 
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Now, we present a proof of lemmafT] The proof utilizes the fact that the number of logical qubits does not depend 
on the system size, and is equal to k even when rti - I. 
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Figure A. 26: (a) Canonical sets of logical operators for the cases with «i = 1 and n\ = 3. A canonical set for ni = 3 can be constructed from 
a canonical set for ni = 1. (b) A canonical set of logical operators. A canonical set for arbitrary ;ii can be constructed from a canonical set for 
ni = 3. 



Proof of lemma IT} Our proof starts with the cases where n,„ are fixed for m > 2. In the proof, we frequently 
change the value of rii while fixing n,„ for m > 2. We represent the stabilizer group for fixed «„ for m > 2 as 

*-*«! = '-'(ni,n2,-,no)- 

First, let us represent a canonical set of logical operators for the case with ni = 1 as follows: 



Tl(S[) = 



Zi, 

^1, 






(A.8) 



where Zp and Xp are Pauli operators defined inside R(l) for p - 1, ■ ■ ■ ,k. (Note that R(l) represents the entire lattice 
since n\ - 1). Next, let us consider the case with n\ - 3. Then, we notice that following operators form a canonical 



set of logical operators (see Fig. A. 26 a)): 



n(>s') = 



Ziri(Zi)r2(Zi), 
liri(li)r2(li), 



z,ri(z,)rf(Zi) 
Xirrx{Xu)T\{Xu) 



(A.9) 



It is immediate to see that these 2k operators commute with all the stabilizers by considering the folding of stabilizer 



generators described in Section 3.2 Since there are only 2k independent logical operators, a logical operator in this 
STS model is equivalent to some product of these 2k logical operators. Thus, for the case with n\ - 3, any logical 
operator £ has a representation which is periodic in the 1 direction: Ti {{) - (,. (This argument proves the translation 
equivalence of logical operators for the cases with odd n\). 

Let us continue to discuss the cases with «i = 3. We consider a bi-partition of the entire system into R{\) and 
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R{2) U R(3). Then, from Theorem[T] we have 

8R(1) + gR(2)UR(3) = ^k. (A. 10) 

Here, we show that all the logical operators defined inside R(l) commute with each other by supposing that there 
exists a pair of anti-commuting logical operators {' and r' defined inside R(iy. {/",r) - 0. Now, let us increase nj. 
Then, t and r' are also logical operators for the cases with rii > 3. Since translations of £' and r' are also logical 
operators, we have 2ni logical operators T^(£') and T-^{r') with the following commutation relations 



r,i(0, ■■■, T"'{e') 

T\(r'\ ■■■, Tl^ir') 



(A. 11) 



for X - 1, ■ ■ ■ ,n\. From lemmal2] these 2n\ logical operators are independent, and we have kj{ > 2n\. However, this 
contradicts with the fact that the number of logical qubits k^t does not depend on n. Thus, all the logical operators 
defined inside R{\) must commute with each other 

A similar discussion holds for logical operators defined inside R{2) U 7?(3). If there exists a pair of anti-commuting 
logical operators defined inside R{2) U ^(3), one can create a large number of independent logical operators for the 
cases with large n\. Thus, all the logical operators defined inside R{2) U RiJ) commute with each other. 

Since there are at most k independent logical operators which commute with each other, we have gR(i) < k and 



gR(2)uR0) ^ ^- Then, fromEq. (A.lOi, we have 

gR(\) = M2)u«(3) = k- (A. 12) 

Now, let us represent the canonical set of logical operators for iii = 3 as follows: 

n(^^) = {^;' :::; i] (a.i3) 

where (p is a logical operator defined inside R{\) fox p - 1,- ■ ■ ,k. Since all the logical operators for ni - 3 have 
representations which are periodic in the 1 direction, one can choose r^ such that 

Ti(rp) = rp (p ^ l,---,k). (A. 14) 



Thus, we have proven lemmafTlfor the case with ni = 3 (Fig. A.26 'b)). 

Here, we represent r^ as rj, = r' Ti{r')T^(r') for ni = 3 where r' is a Pauli operator acting on ^(1). Then, by 
modifying the definition of r^ a little, one can obtain a canonical set of logical operators for any ni . Here, we redefine 



rp for arbitrary «; as (Fig. A.26 b)) 

"1 

(A. 15) 



nw.)- 



We note that this new definition of rp includes the previous definition of rp originally given only for the case with 
Ml = 3. Then, we obtain the canonical set of logical operators for arbitrary n^: 

n(^,;).{^;; ;:;; ^; }. (A.16) 

This canonical set of logical operators has the form described in lemma[T] Though we started our discussion for fixed 
n,„ for m > 2, one can apply the same discussion for any n. This completes the proof of lemma fl] 

Appendix A. 3. Extension of lemmaU] 

One can also extend lemma [T] and obtain the following lemma, which will be useful in the later discussion in 
computing logical operators in two-dimensional STS models. In the below, Q{1) and Q{2) denote regions of ni x 1 
and 1 X n2 composite particles. 
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Lemma 3. 1. Let the number of independent logical operators defined inside a subset of qubits A be gA- For a 
two-dimensional STS model defined with the stabilizer group S^ ,«2 '^'^d k logical qubits, all the logical operators 
defined inside Q(l) can be defined inside Q(l) when n2 > 2: 

8q(Y) = 8Q(i) = k. (A. 17) 

2. Consider an arbitrary set of independent logical operators {£i, ■ ■ ■ ,£k} which are defined inside 2(1). Then, 
there exists a canonical set of logical operators n(»S„, „,) which includes [ti, ■ ■ ■ ,{i^]: {{i, ■ ■ ■ ,/'j.) C n(»S„,_„,) 
in such a way that 

n(5„,„,) = p|; :::; M (a.i8) 



and r\,- ■ ■ ,ri( are periodic: 



Tiirp) = rp. (A. 19) 



The lemma also holds when we interchange the 1 direction and the 2 direction. Here, let us emphasize the 
difference between the statement of this lemma and the statement of lemma[T] In lemma[T] we have shown that there 
exists "some" canonical set of logical operators where fp is defined inside 2(1) while r^ is periodic: T2{rp) = r^. 
However, this lemma ensures that for any given set of k independent logical operators £p defined inside 2(2), we can 
always find rp such that rp are periodic: Toirp) - rp. We skip the proof of lemmalSlsince the proof can be obtained 
with a little modification to the proof of lemmafl] 

Appendix B. Stabilizers in one-dimensional STS models 

In this appendix, we obtain stabilizers for one-dimensional STS models described in Section |4] Consider a one- 
dimensional STS model with ni composite particles. We list all the independent stabilizers defined inside two con- 
secutive composite particles Pi and P2 as 

Sj = [aj,f5j]. (B.l) 

Here, a^,/?/ represents the form of the stabilizer S j graphically, and should not be confused with a commutator We 
compute the overlapping operator group for a composite particle Pi If32l , defined in the following way: 

O''^ ^ ( { U\p, : U e S,„ } ) = ( { «„y6j,V; } ). (B.2) 

Let us represent the overlapping operator group in a canonical form. Due to translation symmetries of the system, we 
have 

[aj,Pj,] = for ally,/. (B.3) 

Here, by {aj,^^} - 0, we mean that Uj andySy commute with each other. Then, in finding a canonical representation 
of (9^', one can consider the contributions from a j and the contributions from/3y separately. Here, we consider two 
groups of operators 0(a) and Oifi) 

0(a) = {{aj,Vj}}, Oip) = {{/3j,Vj}.) (B.4) 



We have described S j as in Eq. (B.l 1. However, the choice of S j is not unique, and one can choose any set of 



S j as long as Sp^up^ - (S j Vj). Then, one can choose S j such that the canonical representations of <9(q') have the 
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following form: 

0(a) ^ I "'' ■■■' "" ^'' •••' ^^'' ^'^ ] (B.5) 

\ a,+u ■■■ , ait, I 

where 

Sit^j = [ipij] for ;■ = l,---,fc', (B.6) 

and t is some integer For a derivation, see t32|. Here, Sp^ represents all the independent generators for the restriction 
of >S into Pi, and fj are logical operators defined inside Pi. 52/+) are stabilizers due to the translation equivalence of 
logical operators. We note that k' < k since there are at most k independent logical operators which commute with 
each other. Then, one can represent 0(J3) in the following form: 

Oi/3)J^'^ ■■■' f" '" ■■■' '*" ^"^ ]. (B.7) 

One can easily verify thatyS^ satisfies the above commutation relations. Having represented 0(a) and 0(P) in canonical 
forms, it is immediate to represent (9^' in a canonical form: 

->p, _ I o^u ■■■, a,, ySi, •■■, /3,, l\, ■■■ , tk'- 



\ ff^+i, ■■■ , a2,, Pt+i, ■■■ , Pit, I 

Now, let us prove that k' - k. Here, for simplicity of discussion, let us assume that there is no stabilizer in 
Sp^. Then, one can represent the number of independent generators for the stabilizer group S„^ by using the above 
canonical representations. First, stabilizers S j for j - 1, ■ ■ ■ , 2f and their translations are always independent. Second, 
stabilizers S j j -2t + \,- ■ ■ ,k' and their translations are not independent since 

W Tl(Sj) = / (; = 2f + 1, ■ ■ ■ ,2f + k'). (B.9) 

X 

Then, we have G(S„^) = 2f«i + k'(ni - 1). Therefore, k' - k. 

In summary, one can represent the overlapping operator group defined for Pi in the following way: 

^P, ^ / «,, •••, a„ A, ■■■, A, A, ■■■, 4 \ (g^Q) 

\ at+u ■■■ , 0-2,, /3t+u ■■■ , yS2„ / 

Here, we neglected stabilizers inside Sp, . Then, after applying some local unitary transformations, one can represent 



S j as in the forms presented in Section 4.2 



By extending the above discussion to the cases where D > 1, one can obtain the following lemma which will be 



useful in discussions in Appendix E 



Lemma 4. Consider a translation symmetric stabilizer Hamiltonian which is invariant under unit translations of 
composite particles: 

T,„(H) = H (m = 1, ■■■,£)) (B.ll) 

and whose interaction terms are defined inside hypercubic regions with 2 X ■ ■ ■ X 2 composite particles. Let ft' be the 
size of the entire system which is defined with n'j X ■ ■ ■ X n^ composite particles where n^ are arbitrary fixed integers. 
The stabilizer group is denoted as Sf where n' — (n'^, ■ ■ ■ , n^). The number of logical qubits is denoted as k' . Let us 
assume that all the logical operators satisfy the translation equivalence of logical operators in the 1 direction: 

Ti(€) ~ i for alii e L^,, (B.12) 

where Ti represents unit translations of composite particles, and Lfj, is a set of all the logical operators in the stabilizer 
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Hamiltonian. Then, there exists the following canonical set of logical operators: 






4 

rk 



(B.13) 



where (p is defined inside a region with I X ny ■ ■ X n'j^ composite particles, and rp is periodic in the 1 direction: 
T\{rp) = rp. 

The proof can be obtained through a similar discussion given in the above. For example, when D - 2, let us 
denote the region with I xn2 composite particles as Q{2). Then, we list all the independent stabilizers defined inside 
e(2)uri(e(2))as 



Sj = ajTi(/3j) 



(B.14) 



where Oj and /3j are Pauli operators defined inside 2(2). By computing the overlapping operator group O^^^^ through 
Uj and /3j, one can easily prove the above lemma. The generalizations to the cases where D > 2 is immediate. Thus, 
we shall skip the proof of lemma B] 

Appendix C. Logical operators in two-dimensional STS models 

In this appendix, we obtain a canonical set of logical operators for two-dimensional STS models which is presented 
in Section 5.2 We begin by introducing some lemmas which are useful in obtaining a canonical set of logical operators 
in Appendix C. 1 Then, we obtain a canonical set of logical operators for two-dimensional STS models in Appendix 



Here, we remind some notations again (Fig. |C.27| i. An STS model has «i x na composite particles and k is the 
number of logical qubits. Qil) and Q{2) represents regions with n^xl and 1 x n2 composite particles. Ua^i, represents 
a region with ax b composite particles where I < a < n[ and I < b < ni- Note that 2(1) = f^«i,i and 2(2) = t/i,„,. 
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Figure C.27: Some regions of composite particles and their notations. 



Appendix C.l. Some lemmas for obtaining a canonical set of logical operators 

In obtaining a canonical set of logical operators, we frequently use the following lemma. 

Lemma 5. For a stabilizer code S with a canonical set of logical operators 



n(S) 






4 

rk 



(C.l) 



consider some logical operator { where the commutation relations between I and logical operators in Y\iS) are 



tp = {-\f"lpi. 



0,1, 



irp = (-D-^T/, /p - 0,1. 



Then, E is equivalent to the following logical operator: 



Wp'i'- 



(C.2) 



(C.3) 
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Proof. Given a logical operator u, u is equivalent to some product of ^i , ■ ■ ■ ,4 and ri , ■ ■ ■ , r^: 

k 



'~U 



P "P 



^p^fp 



0,1. 



(C.4) 



/'=! 



By considering the commutation relations between t and (\,- ■ ■ ,4 and ri, ■ ■ ■ , r^^, we have e^, = e' and /;, = /'. D 

The following lemma becomes also useful in obtaining a canonical set of logical operators. 

Lemma 6. Consider an STS model defined with the stabilizer group »S„,_„j where n\ is odd. Then, for each and 
every logical operator (,, there exists a logical operator € which is equivalent to I and is periodic in the 1 direction: 

Proof Given a logical operator (, the following logical operator is equivalent to ( due to the translation equivalence 
of logical operators: 






x=\ 



since «i is odd. Note that I' is periodic in the 1 direction: Ti{{') - (' . This completes the proof of the lemma. 



(C.5) 



D 



Now, let us begin finding a canonical set of logical operators. Our goal is to show that there exists a canonical 
set of logical operators in two-dimensional STS models which have the configurations of Pauli operators described in 



Section 5.2 Since the canonical set given in Section 5.2 can be used as a canonical set for any «i and n2 as long as 
«! and n2 are larger than 2v, it is sufficient to show the existence of such a canonical set for specific n^ and «2. Here, 
we consider the case with n\ -2- 2^'!, where "!" represents a factorial: 2^'! = 2^^' x ■ ■ ■ x 1. The reason why we take 
such an artificial value as 2 ■ 2^'! will become clear in the course of computations. 
From lemma [T] one can find a canonical set of logical operators 



n(5„,,„,) = 



r\. 



4 
rk 



(C.6) 



where ip are defined inside Q{\) while r^ are periodic in the 2 direction: T2(rp) - rp (Fig. C.28 i. We first analyze the 
properties of logical operators £p by proving the following lemma. 



Lemma 7 (Extraction of periodicity for £p). Consider a two-dimensional STS model defined with the stabilizer group 
Sni,n2 where n\ —2- 2^^'\. If a logical operator £ is defined inside Q{\), there exist centralizer operators t" , t^ € Cq(\) 
such that 



i~ei 



where 



rf(^^) 



(b < 2^"), 



(C.7) 



(C.8) 



and €" is defined inside U2v,i- Here, (" or (^ may be an identity operator. 

Here, Cq{\) is a restriction of the centralizer group C into Q{\). Cq(\) includes stabilizers and logical operators 
defined inside 2(1). Any logical operator ( defined inside 2(1) can be decomposed as a product of two centralizer 
operators t" and ^^ where t" is defined inside a zero-dimensional region f/2v,i and fi is periodic in the 1 direction 
with the periodicity b (Fig C.28 i. In other words, a logical operator { can be decomposed into a periodic part f^ and a 
non-periodic part t" . 

The proof of this lemma is separated into three steps, summarized in three sublemmas respectively. The proof 
relies on the fact that there are 2v independent generators for the Pauli group acting on a single composite particle. 
We begin by proving the first sublemma. 
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I I ~r- 



Figure C.28: Decompositions of logical operators into a periodic part and a non-periodic part. A logical operator i defined inside Q(l) is decom- 
posed as a product of C" and £^ where i" is defined inside U2v, i while {^ is periodic in the 1 direction. A similar decomposition for r is also shown 
where r, r" and r^ are periodic in the 2 direction. 



Subleinina 1. For any logical operator { defined inside Q(\) with ni - 2 ■ 2^^'!, there exist centralizer operators 
r, eP e Ce(i) with 



€ ~ rf where rf(/) ^ f (b < 2^"), 



(C.9) 



and £" is defined inside f/„,-i,i. 



The sublemma claims that £ can be decomposed into a periodic part {^ and a non-periodic part £" where £" is 
defined inside a region U„^-i^i. In the proof, we shall frequently use the translation equivalence of logical operators. 



Proof. Let us represent { as 



i = ||rf~'(Mx) = [mi,M2, ■■■ ,M„,]. 



(CIO) 



Here, an ni x 1 matrix (vector) represents Pauli operators supporting £ inside 2(1) graphically. We prove the sublemma 
for the case with Ux ^ I for all x since otherwise £ satisfies the condition described in the sublemma. 

Since there are at most 2^'' different Pauli operators which act on a single composite particle, there must exist 
integers 71, 72 e Z„^ {ji < J2) with J2 - j\ = b < 2?'' such that m^-, = Uj^. If we represent I explicitly with Uj^ and Uj^, 
we have 



\u 



1 » ■ ■ ■ •> ^ji » ■ ■ ■ ? ^ji-i ' ^72 ' ' ' ' 



,M„,]. 



(C.ll) 



Since Mjj = Mj^, one can construct a centralizer operator by using a sequence of Pauli operators My, ,■• • ,«;,-! appearing 
in {. In particular, one can construct the following centralizer operator: 



njh-l 






(C.12) 



where J2 = j'l + b. Here, Uj^,- ■ ■ , m^^-i appears periodically in £^. Note that ni/b is an integer since ni - 2 ■ 2^"!. 
{^ is periodic in the 1 direction: T'^({^) = {^. One can show that {^ is a centralizer operator by seeing that £^ 
commutes with all the stabilizers. Here, we emphasize that such a construction of a centralizer operator is possible 
since stabilizers can be defined insider regions with 2x2 composite particles. 
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Since both ^^ and i consist of a sequence of Pauli operators Uj^,- ■ ■ , Uj^^i, one can cancel this sequence of PauH 
operators in £ by applying some translation of {^ to {. A resulting centralizer operator is defined inside Um-h,!, which 
is included by t/„,-i,i since b > I. Thus, due to the translation equivalence of logical operators, the decomposition 
described in the sublemma is possible. This completes the proof of the sublemma. D 

Now, one can see the reason for setting ni =2-2^^'!. This choice of iii makes sure that one can extract a periodic 
centralizer operator ^^ from i by requiring that nilb are always integers for b < 2^''. The reason for the extra factor 
"two" will become clear later. 

Next, we need to prove that a logical operator ( inside t/„,-i,i can be also defined inside U2v,i- We prove the 
following sublemma. 

Sublemma 2. When tii — 2 ■ 2^^'\, a logical operator £ defined inside U„^-i^i can be also defined inside U2^',\: 

^C/„,-U = 8V,2.,y (C.13) 

Proof. Let us consider the case where a logical operator { satisfies 

Ml = / and u„ + l{x + 1). (C.14) 

Then, there must exist integers ji, J2 e Z„2 {j\ < J2) with J2 - ji = b < 2^'' such that My, = Uj^. If we represent £ 
explicitly, we have 

i = [/,M2, ■■■ ,Mj, ,Mj,+l--- ,M;2,--- ,M„,]. (C.15) 

Then, one can construct the following centralizer operator: 

f = [/,M2,--- ,Mj,_l,Mj,,My2 + l,My2+2--- ,««,,/,■■■ ,/]■ (C.16) 

We constructed f by discarding m^ for x = ji + 1, ■ ■ ■ , J2 - 1 and attaching M2, ■ ■ ■ , My, to Uj^,- ■ ■ , m„, . One can 
immediately see that £' is a centralizer operator defined with n^ - b - I composite particles. 

Now, both { and {' consist of a sequence of Pauli operators M2 , ■ ■ ■ , Mj, - 1 , Mj, . Then, if we apply f to {, a resulting 
centralizer operator is defined with ni - ji composite particles. Then, due to the translation equivalence of logical 
operators, one may notice that { can be also defined inside L'^„,-2,i since j\ > 2. One can use the same discussion for 
a logical operator defined inside Unt-2,i and shorten its length by one. By repeating this shrinkage until the length of 
{ becomes 2^^', we complete the proof of the sublemma. D 

One can further shorten the length of £, as summarized in the following sublemma. 

Sublemma 3. When ni — 2 ■ 2^^'\, a logical operator I defined inside f/22>",i can be also defined inside t/2v,i-' 

^£/,2v., = ^£/2v,l- (C.17) 

Proof. Let us suppose that i is defined inside t/22»,i with m^ ^ / for x = 1, ■ ■ ■ , 2^'': 

e = [mi,M2,--- ,M22v,/,--- ,/]. (C.18) 

If there exists an integer j such that Uj - Uj+] (1 < j < 2^^' - 1), one can construct the following centralizer 
operator: 

f = \ui,U2,- ■ ■ ,Uj,Uj+2,' ■ ■ ,U22'',I,- ■ ■ ,l\ (C.19) 

whose length is 2^'' - 1. Then, the length of {£' is 2^'' - j. Thus, { can be defined inside a region J/22'-i,i- 

Next, let us consider the case where there does not exist any integer j such that Uj = Uj+i. Since there are 
2v independent generators for the Pauli operator group defined for a single composite particle, there exist integers 
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1 < b < a < 2v such that {ui,, ■ ■ ■ ,Ua} consists only of independent Pauli operators, and 

Mfl+i e <{«*,■■■ ,Mfl)) and m„+i i ({m^+i, ■ ■ ■ ,m„)). 
Then, there exists some set t of integers which is a subset of aset{b,b + I,- ■ ■ ,a - l,a} where 



(C.20) 



(C.21) 



with let. Note that t must include b due to Eq. (C.20 1. 

By constructing a certain centralizer operator based on t, one can decompose ^ as a product of smaller centralizer 
operators. Here, in order to get some intuition, we consider an example where b - I, a - A and t = {1, 3,4). Then, 
we have 



M5 — U1U2U4, 

Then, consider the following centralizer operator 

f = T*(e)T^(£)Ti(£y 

— [I,I,I,I,Ui,U2,- ■ ■] X [I,I,Ui,U2,U2,U4,- ■ ■] X [I,Ui,U2,U2,U4,U^,- ■ ■] 

X [mi, M2, M3, M4, M5, M6, ■ ■ • ] 

= [mi, U1U2, M1M2M3, M2M3M4, M1M3M4M5, U2U4U^U(j, • • • ] 

= [mi, M1M2, U1U2U2, M2W3M4,/, M2M4M5M6, ■ ■ ■ ] 

since M1M3M4M5 = /. Then, one may notice that 

£" — [mi, M1M2, M1M2M3, M2M3M4, /,/,■■• ] 



(C.22) 



(C.23) 



(C.24) 



is a centralizer operator since {" commutes with all the stabilizers. By applying {" to £, one can shorten the length of 
£ by one. Thus, it is possible to create an "eraser" {" by taking a product of f and its translations. 

Now, let us discuss more general cases. Let us consider the following centralizer operator as an eraser: 



r = n 



^a-{x-l) 

1 



iO- 



(C.25) 



Then, since {'\p^^^ , = /, one can construct a centrahzer operator 



■] 



(C.26) 



where Mj , ■ ■ ■ , m"^j are some Pauli operators which can be computed from Eq. (C.25 1. Thus by using {", one can 



shrink the size of (. 

One can use the same trick to shorten the length of a logical operator £ as long as the length of f is larger than 2v. 
Thus, { can be defined inside t/2v,i- D 

By integrating these three sublemmas, one can prove lemma IT] One can also obtain a similar lemma for logical 
operators r^ in Eq.(C.6 1 (see Fig C.28 1. 



Lemma 8 (Extraction of periodicity for r^). Consider a two-dimensional STS model defined with the stabilizer group 
<S{rn,n2) with «! = 2-2^"!. If a logical operator r is periodic in the 2 direction: T2{r) — r, there exist periodic centralizer 
operators r", r" 



Tiir") 



r2(/) = / 



(C.27) 
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such that 

r ~ r"i^ where T^ii^) = r^ (^ < 2^") (C.28) 

and r" is defined inside U2v,m- 

We skip the proof of lemma [8] since the proof is essentially the same as the proof for ip due to the periodicity of 
rp in the 2 direction. As a result of these two lemmas on th e properties of ( and r, both { and r can be separated into 



non-periodic parts f", r" and periodic parts f^, r^ (Fig C.28 1 



Appendix C.2. Canonical set for a two-dimensional STS model 

Having introduced necessary tools in the previous section, let us now give a canonical set of logical operators for 
two-dimensional STS models. Recall that we have limited our consideration to the cases where n\ - 2 ■ 2^'! in the 
previous subsection. In addition, we consider the case where n2 is odd first. In order to indicate that we are considering 
STS models for some fixed n\ and n2, we use notations n\ -hi = 2 ■ 2^'! and n2 = n2 where h2 is some fixed odd 
integer. Thus, we consider the case with ni -hi and n2 - n2 first. 

From lemma [3] there are k logical operators ^i, ■ ■ ■ ,4 defined inside 2(1). Let the number of logical operators 



defined inside f/2v,i be gy^ , , = ko (for n^ -hi and n2 — h2). Then, we can choose ti,- ■ ■ ,4 such that (Fig. C.29 ^b)) 
• ^1, ■ ■ ■ , 4o are defined inside f/21,1. 



t*:o+l> ■ ■ ■ >tt 



are defined inside 2(1) with T^''({p) = {pfor p - ko + 1,- ■ ■ ,k for some integers bp < 2 . 



Here, we attach the indexes "a" and "/3" to represent whether the corresponding logical operator is periodic in the 1 
direction or not: 

^■■■'4"„ and <,„■■■,< (C.29) 

where £" = (p for p - 1,- ■ ■ ,k() and Cp = ( p for p - ko + I , ■ ■ ■ , k. Then, from lemma 3 there exists a canonical set 
of logical operators which includes {",■■■ ,£'^ and ^^ ^j , ■ ■ ■ , ^: 

{pa pa S f^ \ 

«i' ■■■' **„' Vp ■■■' hV (C.30) 

ri, ■■■ , ?•*„, ?-i„+i, ■■■ , r^ ) 

where rp are periodic in the 2 direction: T2(rp) - rp for p = 1, ■ ■ ■ ,k. Here, Yl{Sn,,m) represents a canonical set of 
logical operators for ni = ni and n2 = n2- We shall see that t",- ■ ■ ,i'^ are independent logical operators and g2v,\ — ko 
for any ni and n2 in the course of discussions. 

For now, we concentrate on the analysis on r,t„+i,-- ,r,t. By using lemma [8] we can decompose rp for p - 
ko + I,- ■ ■ ,kin the following way (see Fig. |C.28| i: 

rp ~ r;4 ( r2(r;) = r; and T2(4) = 4 ) where T^'X^) = 4 (cp < 2^'-). (C.31) 



Here, r" are defined inside U2v,m (Fig- C.29 'b)). Now, we show that [/^, ry] - for any p and q ip,q - ko+ I,- ■ ■ ,k). 
We have T^'i^p) - 4 and T^'C^) - 'q- Take the lowest common multiple ofbp and Cq as iip,j. Then, we have 

Tl'-^^p) = el, t'-'(4) = A^. (C.32) 

In other words, fp and r^ have a common periodicity dp^q. Since hi/dp^q is an even integer for ni = 2 ■ 2^'!, we must 
have [€i,,rq] - 0. (This is the reason for the extra factor "two" in 2 ■ 2^'!). Now, commutation relations between r" 
and ^q are 



r 



*°^'' ..;,[. I ^ Lft„«,. (C.33) 
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(a) 



^1' 






2v 




X„ 



r^ V 



_ t 

— — ►^i 



*- 1 

Figure C.29: (a) Logical operators f" and f^. {" are defined inside [/2v,i- ^ are defined inside 2(1), and periodic in the 1 direction: T^'iffp) = ffp. 
Logical operators Vp for p = Uq + i,--- ,k are anti-commuting pairs of ^, and are periodic in the 2 direction, (b) HJ and their anti-commutation 
with dj, for p = <:o + 1, ■ ■ ■ , ^. A two-sided an'ow represents the anti-commutation: [lfp,rVt = Q. r^ can be represented in a periodic way inside 
2(2): r« = n;'i, 7'2^(^,,). 



Also, r" commutes with /" since I" can be defined so that I" has no overlap with r" due to the translation equivalence 
of logical operators. Thus, we obtain the following set of independent logical operators (see Fig. C.29 a)(b)): 



{1, 






^-ko + l' 
ko+V 



C Lh 



(C.34) 



where Lsj^^ is a set of all the logical operators for n\ - hi and n2 - n2- Since all the logical operators defined inside 
(2(2) can be defined inside Q{2) from lemma 3 logical operators /"^ +1 , ■ ■ ■ ,r^ can be defined inside 2(2). Let us recall 
that we have been considering the case with oad «2- Given 4* +; , ■ ■ ■ ,r1 defined inside 2(2), we can represent r" such 



that they are periodic in the 2 direction from lemma ^ since 



IXmr';) 



(C.35) 



due to the translation equivalence of logical operators. Let us choose r" in the following way (Fig. C.29 b)): 



= fl ^2(^P) - 



x„ 



(p^ko + l,--- ,k) 



(C.36) 



where Xp is some Pauli operator acting on a composite particle Pii . Xp commutes with Xg for p,q - 1 , ■ ■ ■ ,ko. Now, 



ki) + 1, 



k, defined in the forms in Eq. (C.36i, are independent logical operators for 



let us show that r^ for p 

any 112 and ni = hi - 2 ■ 2^'!. r" are centralizer operators for any ht. Since r" anti-commutes with €p, they are not 

stabilizers. Thus, the following logical operators are independent for any 112 and n\ = hi = 2 ■ 2^' !; 



c Lh, 



,V«2 



'/to+1' ' k 

Here, L^, vn^ means that operators on the left hand side are logical operators for any n2 and rti = hi. 



(C.37) 
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and rf , , , 



We continue to consider the case where ni -hi - 2 ■ 2^'! and na - na- Since t'^,--- ,C 
are independent logical operators which can be defined inside 2(2), by using lemma pi we notice that there exists 
following canonical set of logical operators: 



^{Sn,,ny 



pa 









•_k 

Zk 



(C.38) 



where Zp are periodic in the 1 direction: Ti{Zp) - Zp. Here, we denote this new canonical set as n(>S,-,,, ,-,,)' to make 
the difference from n{Spi^fi^) clear. Recall that Zp are periodic in the 1 direction. Since na is odd, we can define Zp 
such that they are periodic both in the 1 and 2 directions: 



Un(Zp). 



(C.39) 



Now, through some calculations of commutation relations, we notice that the following logical operator form a canon- 
ical set: 



n(Sr„,i,r 



Zu 






'h+v 

kn + V 



•k 



(C.40) 



Here, we notice that I" and Z^ (/? = 1, ■ ■ ■ , A:o) are independent logical operators for any n\ and n-i since they are 
centraUzer operators and anti-commute with each other for any n\ and n^. Let us represent Zp as 



^p = nn^>^2(f^.) 



A=l y=l 



Uo, u„ 



Up, Up 

u„, u„ 



u. 



u 



pj 



(p = 1, 



,^o) 



(C.41) 



where Up is a Pauli operator acting on a composite particle Pii. From geometric shapes, Zp (p - I,- ■ ■ ,ko) are 
two-dimensional logical operators, and t"{p- 1, ■ ■ ■ , A;o) are zero-dimensional logical operators. 

While two-dimensional logical operators Zp commute with each other: [Zp,Zp>] = foT p,p' = 1, ■ ■ ■ ,ko when 
«i = «! and «2 = «2> they may anti-commute when both n^ and «2 are odd integers. In order to find two-dimensional 
logical operators which commute with each other, we consider the case where «i and «2 are some odd integers. To 
indicate the difference, we denote rii = hi and n2 = "2 where hi and n2 are odd integers. Then, Zp may not commute 
with each other since Up may not commute with each other for p = 1, ■ ■ ■ ,ko. Let us represent the commutation 
relations between Zp as follows: 



7 7 



(-ir-Z.Zp, e 



'P.q 



±1 for p,q — \, 



{CAT) 



Here, we define the following two-dimensional logical operators: 



/p ^Zp\\{(ir», 



K ^ nn^i^2^^p ~ ^p- 



?=i 



(C.43) 



We note that l" is periodic both in the 1 and 2 directions. Then, we have a following set of independent logical 
operators (see lemma|5]i: 



pa 






C La 



(C.44) 



1' ' ito' ko+V ' 'k 

Here, ^f , • • ■ , ^? are zero-dimensional logical operators and r.,- ■ ■ ,ri are two-dimensional logical operators. In 
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particular, they are logical operators regardless of «i and «2: 



pa 






C Lv 



(C.45) 



Here, Lv„,,vn, means that {" and r^ are independent logical operators for all n\ and 112- 



Let us summarize our discussions so far. In particular, by combining Eq. (C.37i and Eq. (C.45i, we have the 
following independent logical operators for any ni > 2v and n\ -2-2^^'\: 



n(>Sg, Vn,) 



^p 



' ''to' to+I' 
' in' *0 + l' 



(C.46) 



where n(»Sii,,vn2)" means that this canonical set is valid for n\ - hi and any n2- One can represent each logical 
operator in the following way. 

• {" are zero-dimensional logical operators defined inside t/21,1 for p - I,- ■ ■ ,ko. 

• rp are two-dimensional logical operators which are periodic in both directions: Txir],) - T2(rp) - 'p with 



A = nn^i^2(^^^> 



^n, ^n 



^n^ ^n 



X,, 



.x,„ 



x„ 



x^ 



pj 



(C.47) 



for;? - 1, ■ ■ ■ ,kQ. 



r"j are one-dimensional logical operators defined inside 2(2), and periodic: T2{r",) - r" with 



n ^2(^/') 



x„ 



Xp, 
.Xp, 



I, 
I, 



(C.48) 



for p - kQ + \, 



€p are one-dimensional logical operators defined inside 2(1), and periodic with periodicities bp: T^''(€p) - {!], 



where bp < 2^\ 



Here, we notice that logical operators ("{p- 1, ■ ■ ■ , ko), r"{p - ki) + I,- ■ ■ ,k) and rp (p - 1, ■ ■ ■ , Aro) have the forms 



described in a canonical set of logical operators described in Section 5.2 

To complete the proof of the theorem, we need to analyze Cj,. In particular, we need to show that /^ can be defined 
in a periodic way: Ti(€p) - Cp. In fact, it is immediate to show this by considering the case where n2 = 2 ■ 2^'! and «] 
is odd instead of the case where ni = 2 ■ 2^'! and «2 is odd. The proof basically follows the same discussion used in 
the analyses on the case ni = 2 ■ 2^^ ! and ^2 is odd, and thus, we skip it. This completes the derivation of a canonical 
set of logical operators. 



Appendix D. The existence of loop-like stabilizers 

In this appendix, we sketch the proof for the existence of stabilizers which are used to construct loop-like stabi- 



lizers appeared in Section 5.2 For simplicity of discussion, we discuss the cases without zero-dimensional logical 
operators (ko = and ki - k). Generalizations to the cases where ko + Q are immediate. 
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Let us recall that one can represent all the one-dimensional logical operators as follows: 



04' 



Zp, /, 



[Zp, I, 



n^^ 



I, I, 
I, I, 

Xp, Xp, 



X. 



/'J 



(D.l) 



for p — 1 , ■ ■ ■ ,k. These 2k operators are independent logical operators regardless of ni and 112- 

We represent independent stabilizer generators which can be defined inside 2x2 composite particles in the fol- 
lowing way: 



^(U) 



Ca, Da 



n(U) 



(D.2) 



for q — 1, ■ ■ ■ , v' with v' > v where v is the number of qubits inside each composite particle. Stabilizers S g are well 
defined for the cases where n\ - 1 or «2 — 1, as discussed in Section 3 For example, when «2 = I, S q may be 
represented as 



:(/,!) 



= [aq,liq\ 



ao 



(D.3) 



Ai,l) 



where a^ = A^Cq andPq = BqDq. S q is defined inside a region with 2x1 composite particles. 
For convenience of discussion, we introduce the following notations: 

qeR ijeR ^eR 



(D.4) 



where R represents a set of integers in Zy' . We also define the summation of sets of integers R and R' such that 

o-CReR') = aCR)a'(R'). (D.5) 

Thus, "©" satisfies the following properties: 



7 e R ;■ e R' => j e R © R' 

;■ e R ./■ ^ R' => j ^ R © R' 

y ^ R ./■ e R' => j ^ R © R' 

j ^ R ./■ ^R' => ;eR©R'. 



(D.6) 



We begin by proving the following lemma. 
Lemma 9. For a two-dimensional STS model, consider the cases with «i X 1 composite particles fni > 2 and n2 — I). 



Then, there exist sets of integers Rp (p - 1, 



, k) such that 

a{Rp) = PiRpX 



and 



[a{Rp)J,-- ,/] ~ £p ^ [ZpJ,-l\. 



(D.7) 



(D.8) 



Proof. Let us consider the case with n\ -2 and n2 = 1 (2 x 1 composite particles). In this case, logical operators are 
(p - [Zp,I] and rp - [Xp,Xp]. Here, we represent Pauli operators acting on the system of 2 x 1 composite particles 
through two component vectors. These representations should not be confused with commutators. 
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Due to the translation equivalence of logical operators, Ti({p){p is a stabilizer: 

Ti({p){p = [Zp,Zp] eS2,i. 



(D.9) 



Here, S2,i is the stabilizer group. Then, this stabilizer can be decomposed as a product of stabilizers S ^ and their 
translations. In particular, there exist some set of integers R^ and R' such that 

[Zp,Zp] = 5(''i'(R^)ri(5(''i>(R;)) = KR^)/3(R;),a(R;)/?(R^)] 

where 5('''>(Rp) = [a(Rp),/3(Rp)] and 5<''''(R;) = [a(R'p),/3(R'p)]. Here, we have 

[Zp,I] = [a(R^)^(R;),/] = [a(R'pWRp),I] e Cz.i- 

Also, let us consider the following stabilizer: 

5<''i'(R^eR;) = [a(Rp®R'^),/3(Rp®R;)]. 

Then, since 

Zp = aiRpMR'p) = a(R;)/J(Rp), 

we have a(Rp e Rp = /3(Rp e Rp, and 

[a(Rp®R'pXI] = ^(R^eRp,/] e Ci.i. 



(D.IO) 
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(D.13) 
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Then, we notice that [a(Rp)/3(Rp),I] e Cij. Since S^''^\Rp) = [a(Rp),/3(R^)], we also have [a(Rp),Il\fi(Rp),I] e 
C2,i. By using a similar discussion, we eventually have 



[a(Rp,/], KRp,/], b8(RP,/], L0(Rp,/] e Cz.i. 



(D.15) 



Since S'-''^\Rp) = [a;(Rp),yS(R;,)] e >S2,i, we have [a(Rp)f3(Rp),I] e >S2,i due to the translation equivalence of 
logical operators. Thus, we finally have 



and 



5<'''>(R^eRp = [a(R^eRp,a(R^eRp] 



[a(R,,eRp,/] ~ ^(Rpa(Rp,/] = ^^ = [Zp,I]. 



(D.16) 



(D.17) 



Though we have limited our considerations to the cases with 2x1 composite particles, 5*''''(Rp © Rp satisfy the 
properties described in the lemma for the cases with ni x 1 composite particles for arbitrary «i . This can be shown by 
using the fact that {p = [Zp, /, ■ ■ ■ , /] are logical operators regardless of ni. This completes the proof of the lemma. D 

At the beginning of the discussion, we have given a specific set of S g . However, as long as they generate the 
same stabilizer group Sni,m, one can freely choose any stabilizers defined inside a region with 2x2 composite particles 
as S q . Here, we choose S p so that the following operators are independent logical operators: 



\a)p = 



ap, I, 

ap, I, 
a,,, /, 



(D.18) 



for p - 1 , ■ ■ • ,k where ap = ApCp. 
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Next, let us consider the case with ni x 1 composite particles. When «i - 1, S ^ may be represented as 



c(lj) 



yq 



(i,l) 



(D.19) 



where jq = AqBq and 6q = CqDq. Here, we consider the following centralizer operators: 

/, /, ■ ■ ■ , / 

r(r)p = 



yp, 7p, 



I 



(D.20) 



While (,(a)p are independent logical operators for p - 1 , ■ ■ ■ , A:, centralizer operators r(j)p may not be independent 
logical operators. However, through some speculations, one may notice that one can choose S ^ such that both 
t{a)p and r{j)p are logical operators for /? = 1 , ■ ■ ■ ,k. Then, t{a)p and r{j)p are independent logical operators for 
P^ I,--- X 

Though we have obtained 2k independent logical operators, l{a)p and r(y)p do not form a canonical set of logical 
operators. In fact, one will soon see that \t{a)p, r{y)p\ - since \ap,jp\ - 0. In order to represent a canonical set of 
logical operators in terms oil{d)p and r{j)p, let us use the following commutation relations: 

\ap,ap.\ = 0, \jp,yp-\ = 0, \ap,JiA - 0. (D.21) 

Since commutations \ap, ap>] =0 and [y^, jpi] = are obvious, let us show [fia)p, r{y)p] - 0. We use the fact that 



AiJ) 



Cn 



Dp 
B„ 



('■,;■) 



(D.22) 



are stabilizers, and, they commute with their own translations. For p - 1, ■ ■ ■ ,k, we have Dp - ApBpCp. Then, we 
have 



{Ap,Dp] = {Ap,ApBpCp] = {Ap,BpCp] = 0, {Bp,Cp] = 0, 

and, as a result, we have 

[ap,yp] = [ApCp,ApBp] = 0. 
With small extension, one can also show that 

[ff(R),r(R)] = 

for any R where R is a set of integers inside "Lk- Here, ^(R) is defined in a way similar to q'(R). 



(D.23) 



(D.24) 



(D.25) 



Now, let us show that k is even, and the existence of deformers described in Section 5.2 Below, we shall frequently 
change the choices of S q while keeping the group generated from S q . {{S q )'^^^). 

Let us begin by analyzing the case where ^ = 1. Then, we have [Q;i,yi] - 0. However, this contradicts with the 
fact that the following operators are independent logical operators: 

o\, /, ■ ■ ■ , 



a\, I, 



y\, y\. 



I 

n 



(D.26) 



Thus, ki^l. 
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Next, let us consider the case where k -2. Then, the commutation relations between Up and jp must be 



72, n 



Here, we apply local unitary transformations to qubits inside each composite particle in the following way: 



u{ "'' "^ \u-' = 



72, n 

Then, after unitary transformations, we have 

AiCi = Zi, A2C2 = Z2 



Zi, Z2 

'^i, ^2 



A,B 



1£>1 



^2, 



A2B2 = ^1 



(D.27) 
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Thus, with some PauU operators Pi and P2, one can represent S':^ and 5i'^ as follows: 



Al,j) 



P1Z1X2, PiZi 

PIX2, Pi 



(ij) 



sr^ 



Pi, P2X1 
P2Z2, P2Z2X1 



(i,j) 



(D.30) 



which are stabilizers described in Section |S!2l 

Now, we consider the cases where k > 3. Since £(a)p and {ij)p are independent logical operators, there exists 
some integer / ^ 1 such that {o'l,')',) = 0. Without loss of generality, we can set / - 2. Also, one can choose S j',' so 



that [73, O'l] - 0. Then, from Eq. (D.21 1, we have 



O'l, Q'2 

72, 71 



(D.31) 



Also, through some computations of commutation relations, one may notice that the following; if 73 commutes with 
Qf2, then (23 commutes with 72, and if 73 anti-commutes with 02, then aj, anti-commutes with 72. 
When 73 commutes with 02, we have 



When 73 anti-commutes with 02, we have 



Qfi, a2, 03, 73 
72, 71 



O'l, Q'2, a^ai, 7371 
72, 71 



Then, by changing the choice of 5^'^ so that 



:(i,j) 



:(!J) 



AiJ) 



AiJ) 



:(!J) 



o (iJ) c ('.;■) 
o 1 03 , 



we have 



O'l, 02, Q'3, 73 

72, 71 



(D.32) 
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(D.35) 



:(hj) 



Thus, in any cases, by choosing S p appropriately, we have the above commutation relations described in Eq. (D.32 1 
when k> 3. 

Now, if A: = 3, we have a contradiction in a way similar to the case with k - 1. If A: = 4, one can choose S j',' so 



73 



that 



Then, we have the following stabilizers: 



ai, 02, Q'3, a'4 

n, n 74, 73 



5(U) 



^(U) 



P1Z1X2, PiZi 

^1^2, Pi 

Pi, PiXi 

P2Z2, P2Z2XI 



(i,j) 






('../) 



c('j) 



P3Z3X4, P3Z3 
P3X4, Pi 

Pa, PaX^ 
P4Z4, P4Z4X3 



nCj) 



('■jl 



(D.36) 
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One may use a similar discussion for an arbitrary k and show that k must be even, and there exist a choice of S p 
such that 



{i,j) 



a\, ai, 

n, 71, 



ak-i, ak 
7k, Jk-i 



This leads to the existence of stabilizers described in SectionlSSI 



(D.39) 



Appendix E. Scale symmetries and weak breaking of translation symmetries 

In this appendix, we discuss translation symmetries of ground states in STS models in the context of coarse- 
graining. 

A central idea behind RG transformations is coarse-graining of a system of spins by grouping several spins into 
a single particle with a larger Hilbert space. Now, in coarse-graining STS models, the main problem is to find the 
smallest suitable unit cell of qubits to be grouped in coarse-graining. Though one can freely coarse-grain the system 
by choosing an arbitrarily large unit cell of qubits, one may end up with a coarse-grained particle with an arbitrarily 
large Hilbert space, which are hard to analyze. Thus, we naturally hope to choose the smallest possible unit cell of 
qubits. 

The main difficulty in grouping qubits comes from the fact that translation symmetries of ground states may be 
broken. One might hope that a translation symmetric Hamiltonian would give rise to translation symmetric ground 



states. However, this naive expectation is generally true only when there is a single ground state. In Appendix E.2 



we show an example of a stabilizer Hamiltonian whose degenerate ground states are invariant only under translations 
of several qubits, while the Hamiltonian is invariant under unit translations of qubits. Such a bizarre breaking of trans- 
lation symmetries is sometimes called a weak breaking of translation symmetries, which is often seen in topologically 
ordered systems |9|. When a translation symmetry is weakly broken, there exists a ground state whose translation 
is orthogonal to itself. Thus, when coarse-graining the system of qubits, we hope to group qubits so that ground 
states are invariant under unit translations of newly defined particles in order to resolve a weak breaking of translation 
symmetries. 

Here, the whole procedure of coarse-graining may be summarized in the following three steps: 

1. Coarse-grain so that the Hamiltonian is invariant under unit translations of composite particles. 

2. Coarse-grain so that interaction terms are defined inside 2x2 composite particles. 

3. Coarse-grain so that ground states are invariant under unit translations of composite particles. 

We have akeady performed the first and second steps of coarse-graining when the STS model is introduced in Sec- 
tion [3] Here, we hope to further coarse-grain the system of composite particles so that all the ground states are 
invariant under unit translations of coarse-grained particles. 

In this appendix, we show that the composite particle, after the first and second steps of coarse-graining, happens 
to be the smallest suitable unit cell by proving that all the ground states of STS models are invariant under unit 
translations of composite particles. In other words, we do not need to perform the third step of coarse-graining since 
STS models have been already coarse-grained appropriately through composite particles ! 
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In Appendix E. 1 we give an example of a stabilizer Hamiltonian which exhibits a weak breaking of translation 



symmetries and show that this system can be appropriately coarse-grained by introducing composite particles. In Ap- 
Ipendix E.2| we prove that all the ground states of STS models are invariant under unit translations of composite 
particles. 

Having seen that a system of composite particles has been already coarse-grained properly in STS models so that 
all the ground states are invariant under unit translations of composite particles due to scale symmetries, we consider 



the cases where scale symmetries are not satisfied. In Appendix E.3 we extend our discussion to the cases without 



scale symmetries where the number of logical qubits kjj is not constant and depend on the system size n. We show that 
one can always find some finite size of a unit cell of composite particles so that all the ground states are translation 



symmetric when the number of logical qubits k^ is upper bounded by some constant. Finally, in Appendix E.4 
we show that, in translation symmetric stabilizer Hamiltonians, translation symmetries of ground states are weakly 
broken if and only if scale symmetries are broken. 

Appendix E.l. An example of weak breaking of translation symmetries 

A unique ground state of a translation symmetric Hamiltonian always has translation symmetries since the trans- 
lation of the ground state must be the same as the original ground state. However, when ground states are degenerate, 
translation symmetries of ground states may be broken despite the translation symmetries of the supporting Hamil- 
tonian. In particular, there are cases where degenerate ground states are invariant only under translations of several 
qubits, while the Hamiltonian is invariant under unit translations of qubits. This bizarre breaking of translation sym- 
metries, sometimes called a weak breaking of translation symmetries, lies behind the main motivation for further 
coarse-graining the system. 

Here, we give an example of a stabilizer Hamiltonian which exhibits a weak breaking of translation symmetries, 
and show that translation symmetries can be restored by introducing composite particles. 

Breaking and restoration of translation symmetries: Let us consider the following one-dimensional Hamilto- 
nian with periodic boundary conditions: 

// = -^zWz(-'+i>Z(-'+2) (E.l) 

where the total number of qubits is A^, and Z*-'' act on qubits labeled by j. This Hamiltonian is a stabilizer Hamiltonian, 
and has four ground states when A^ is a multiple of three: 

1000000 ■■■ 000), I011011---011), 1101101 ■■■101), I110110---110). (E.2) 

Here, we notice that ground states are not invariant under unit translations of qubits despite that the original Hamilto- 
nian is invariant under unit translations (Fig. |E.30l l. 



Qubit 

io"o' 



N; ' 

Composite particle 



Weakly-broken translation 
symmetry 

Restoration of translation 
symmetry 



Figure E.30: Coarse-graining and restoration of translation symmetries. 

Now, let us coarse-grain the above system by introducing composite particles. Obviously, if one consider three 
consecutive qubits as a single composite particle, one can restore translation symmetries: T^{\il/)) - li/r) where T shifts 
the position of qubits by one. Thus, while translation symmetries of ground states are weakly broken, one may restore 
translation symmetries through coarse-graining. 

Scale symmetries and composite particles: Next, let us see the relation between translation symmetries of 
ground states and scale symmetries. In a strict sense, this model is not an STS model since the number of degenerate 
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ground states changes according to A^. Here, let us return to the original picture through qubits before introducing 
composite particles. When A^ is a multiple of three, there are four ground states with k = 2 since G(S) - N - 2. 
However, when A^ is not a multiple of three, there is only a single ground state with k -Q since G{S) - N. 

In order to treat the model as an STS model, one needs to limit considerations to the cases where A^ is a multiple 
of three. Then, by considering three consecutive qubits as a composite particle (v - 3), one can reduce the system to 
an STS model: ;t„ = 2 for all n = A^/3. 

Now that we have reduced the model to an STS model, one may readily notice that translation symmetries of 
ground states are not broken since all the ground states are invariant under unit translations of composite particles. 
Thus, for this STS model, we do not need to further coarse-grain the system. In particular, we notice that the smallest 
number of qubits to be grouped in order for this model to have scale symmetries is the same as the smallest number 
of qubits to be grouped for restoring translation symmetries in ground states. We shall discuss whether this relation 
between translation symmetries of ground states and scale symmetries holds for arbitrary STS models or not in the 
next subsection. 

Appendix E.2. Translation symmetries of ground states in STS models 

As seen in the above example, even if translation symmetries of ground states are broken, it may be possible 
to enlarge the size of a unit cell to restore translation symmetries of ground states. However, finding the smallest 
suitable unit cell of qubits is generally difficult since ground states are highly entangled in topologically ordered 
systems. In particular, unlike the above example where ground states can be represented simply as direct product 
states, topologically ordered ground states are not direct product states. 

A valuable idea we can import from quantum coding theory is the study of physical symmetries of degenerate 
ground states through the analysis on physical symmetries of logical operators. In the language of quantum coding 
theory, the ground state space of the system Hamiltonian is the codeword space which is stabilized by the stabilizer 
group S. Logical qubits are encoded in the codeword space which is exactly the same as the ground state space. 
Since the encoding of logical qubits can be analyzed by properties of logical operators, it happens that one can reveal 
physical symmetries of ground states by studying physical symmetries of logical operators. 

Here, we show that translation symmetries of ground states can be studied by analyzing how logical operators 
transform under translations. Through the analysis of logical operators, we prove that all the degenerate ground states 
of STS models are invariant under unit translations of composite particles. 

Theorem 4 (Translation symmetries of degenerate ground states). Each and every degenerate ground state li/r) of an 
STS model is invariant under unit translations of composite particles in any direction: 

T„,m) = |<A>, V|.A> ^Vs, (m = I,--- ,D) (E.3) 

up to a trivial phase factor, where Vs^, is the ground state space for an STS model defined with the stabilizer group Sg. 

Below, we sketch the proof of the theorem. We begin by showing that translation symmetries of ground states 
can be studied by analyzing how logical operators transform under translations. In particular, we show that ground 
states are invariant under unit translations of composite particles if and only if the translation equivalence of logical 
operators holds. 

Consider a stabilizer code which is defined on a D-dimensional hypercubic lattice of ni x ■■■ x no composite 
particles. Let us first analyze the necessary condition in order for all the degenerate ground states to possess trans- 
lation symmetries. For this purpose, let us assume that all the degenerate ground states have translation symmetries 
and analyze how logical operators must transform under translations. We start by analyzing the case with two-fold 
degeneracy (k = 1). Then, any degenerate ground state can be represented in the following way: 

\4r} = ck|(Ao>+/3|<Ai) (E.4) 

where 

^l^o) = l«Ao>, mi) = -l^i>, '-I'Ao) = l«Ai), '-I'Ai) = l-Ao) (E.5) 
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with anti-commuting logical operators t and r. Here, from an assumption, li/^o) and |i/ri) have translation symmetries. 
Then one may notice that r,„(^) acts in a way exactly the same as (. inside the ground state space, and one has t ~ T„i({). 
A similar discussion holds for r, and one has r ~ T,„{r). Thus, the translation equivalence of logical operators holds. 
The generalization of the above discussion to the cases with A: > 2 is immediate, and we notice that all the logical 
operators must remain equivalent in order for all the degenerate ground states to have translation symmetries: 

T,„(\ip)) = \iP) for all \4/) e Vs ^ ( ~ T,„({) for alK e L (E.6) 

where V^ is the ground state space and L is a set of all the logical operators. 

In fact, "the translation equivalence of logical operators" is the necessary and sufficient condition for "the existence 
of translation symmetries in degenerate ground states": 

T,„{\i//)) = |(A> for all \4r) e Vs « f ~ T„{{) for alU 6 L (E.7) 

In order to verify the "<=" relation in the above equivalence between ground states and logical operators, we assume 
that all the logical operators remain equivalent under some finite translation: T„ii€) ~ (. For simplicity of discussion, 
we show the existence of translation symmetries in degenerate ground states for the cases with ^ = 1. Let us represent 



ground states of the stabilizer Hamiltonian as eigenstates of t as represented in Eq. ( E.5 i. Then, li^o) can be represented 
in terms of the stabilizer group S and logical operators in the following way 13211 : 

where A^ is the total number of qubits and 5 1, ■ ■ ■ ,5ai_i are independent generators for>S. Then, one can immediately 
notice that T„,(\\I/q)) - |i^o) since T„,{{i,S)) - {(,S). One can also repeat the same argument for |i^i) and obtain 
TmiVfJx)) - I^Ai)- More generally, we have Tn,{\i//}) = \t//} for any ground state \iff}. Thus, all the degenerate ground 
states possess translation symmetries. A generalization to the cases where A; ^ 1 is immediate. 

As a result of the above discussion, we notice that "the translation equivalence of logical operators" is equivalent 
to "the existence of translation symmetries of ground states". Thus, in an STS model, both the system Hamiltonian 
and degenerate ground states are invariant under unit translations of composite particles. 

Appendix E.3. Reduction to STS models: coarse- graining with composite particles 

Originally, we have defined composite particles for convenience of discussion so that the system Hamiltonian 
is invariant under unit translations of composite particles. Fortunately, in the presence of scale symmetries, all the 
ground states are invariant under unit translations of composite particles, and thus, STS models are coarse-grained 
properly. 



However, in Appendix E. 1 we have seen that the number of degenerate ground states may depend on the system 
size in general. In fact, the main challenge is to find a suitable unit cell which restores translations symmetries of 
ground states. From the discussion in the previous section, we know that it is sufficient to restore scale symmetries in 
order to translations symmetries of ground states. However, it is not clear whether it is possible to restore translation 
symmetries of ground states and scale symmetries or not. 

Here, we give a sufficient condition for translation symmetric stabilizer Hamiltonians to be properly coarse- 
grained. We also give a procedure for coarse-graining. In particular, one can show that translation symmetric stabi- 
lizer Hamiltonians can be coarse-grained and discussed in the framework of STS models if the number of degenerate 
ground states are upper bounded by some finite integer, as summarized in the following theorem. 

Theorem 5 (Reduction to an STS model). Consider the case where the number of logical qubits is upper bounded by 
some finite integer k: 

kjf < k for all ft (E.9) 

where the bound is tight for some positive integers «„, = «™" for m — 1 , ■ ■ ■ ,D. Then, there exist some finite positive 
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integers a,,, < 2^'' such that 

^fli«i,fl2«2,-, flo«o =^ forall «!, ■■■,«£, (E.IO) 

where ni, ■ ■ ■ ,no are positive integers. 

Thus, if we consider ai x ■ ■ ■ x ao composite particles as a new composite particle, the stabilizer code can be 
treated as an STS model since the system possesses scale symmetries with newly defined composite particles. 

Here, we describe a procedure to find fli , ■ ■ ■ , ad without giving a proof. It is possible to prove that this procedure 
appropriately coarse-grains stabilizer Hamiltonians through lemma|4] 

. Find ^""■«' = Ur\ ■ ■ ■ . "n"'"^) such that /fc?""' = k. 

• Consider a system with n^'"'"\ and find the smallest integer oi such that 

T"\e) ~( foralU. (E.ll) 

Note that a\ < 2^* and n" /ai is an integer. 

• Consider a system with ft - {a\,n^ , ■ ■ ■ , nl^ ), and find the smallest integer a2 such that 

rfCO -^ for all ^. (E.12) 

• Consider a system with n = (ai,a2,nf ,■ ■ ■ ,n'^ ), and obtain 03. 

• Repeat the same procedure until we obtain oi, ■ ■ ■ , «£>■ 

Appendix E.4. Scale symmetries and weak-breaking of translation symmetries 

In summarizing the discussions so far, we have the following relation between scale symmetries and translation 
symmetries of ground states: 

h ^ k ^ ( ~ T„(£) « T„,im = I.A) (E.13) 

where the above equations represent "scale symmetries", "the translation equivalence of logical operators" and "trans- 
lation symmetries of ground states" respectively. 

A naturally arising question is whether the existence of scale symmetries is the necessary and sufficient condition 
for the existence of translation symmetries of ground states or not. It turns out that scale symmetries and translation 
symmetries of ground states are the equivalent conditions for stabilizer Hamiltonians: 

k,i ^ k o e ~ T„,{{) o T„M)) = |iA). (E.14) 

We summarize the above relation in the following theorem. 

Theorem 6 (Translation symmetries of ground states and scale symmetries). Consider translation symmetric stabi- 
lizer Hamiltonians which are invariant under unit translations of composite particles: 

T,n{H) = H im = I,--- ,D) (E.15) 

and whose interaction terms are defined inside hypercubic regions with 2 X ■ ■ ■ X 2 composite particles. Let kg be the 
number of logical qubits for a system with the size it. Then, the existence of scale symmetries is the necessary and 
sufficient condition for the existence of translation symmetries of ground states: 

kft = k forall it o Tn{\if/)) = |i/r) forall \\jj) e Vs, (E.16) 

where Vs„ represents the ground state space for a system with the size ft. 
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In order to prove this, one needs to show that the translation equivalence of logical operators is a sufficient condi- 
tion for scale symmetries: 

ka ^ k for all n <= { ~ T,„{() for all n (E.17) 

Here, we give only a sketch of the proof since discussions used for the proof are similar to discussions used for the 
proofs of other theorems in this paper The goal is to show that the number of logical qubits is the same for all n. First, 
we consider the cases where n,„ are odd integers for all m. Due to the translation equivalence of logical operators, all 
the logical operators have translation symmetric representations. In particular, after some appropriate local unitary 
transformations on qubits inside each composite particles, we have 

D n,„ D n,„ 

^z' " n n ^""'^^'^' '"" " n n ^™"(^''>- (^-^^^ 

m-l A"m = l m=l j:,„-1 

Then, when n,„ - 1 for all m, Xp and Zp are independent logical operators. This implies that kg must be the same for 
the cases where n„, are odd integers for all m. We denote kjf for the cases where n,„ are odd as k. 

Next, we consider the case where only «i is an even integer and other «„, are odd integers. Let us denote the 
system size as ft' where «'j is an even integer and «J„ are odd integers for m > 1. We also denote the number of logical 
qubits as k' . Then, from lemma|4] one can show that there exists a canonical set of logical operators as follows: 

mS,,) = ( 5' ■■■' 5" I (E.19) 

where £' are defined inside a region with I xn'^x- ■ ■ xn'^ composite particles, and r' are defined in a periodic way in 
the 1 direction: Ti(r') - r' . Then, {\, ■ ■ ■ ,t'^, and r\,- ■ ■ ,r'^, are independent logical operators for the case where the 
system size is ni - 1 and n^ - n'^foT m > 1 since r' are defined in a periodic way. Then, we have k' = k. Thus, for 
any cases where rii is even and all the other «„ are odd, the number of logical qubits is k. 

One can repeat the similar argument for the case where rti and na are even, and other «„ are odd. By using this 
argument, one can show that the number of logical qubits must be k for any ft. 
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